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Introduction to Formal Realizability 
Theory—I 


By BROCKWAY McMILLAN 


Manuscript received October 15, LO51 


This paper offers a general approach to the realizability theory of net- 
works with many accessible terminals. The methods developed are applied 


to give a complete characterization of all finite passive networks. 


I. SUMMARY 

1.0 A principal result of this paper is to characterize those matrices 
Z(p), functions of the frequency parameter p, which can be realized as 
open-circuit impedance matrices of finite passive networks. This char- 


acterization is provided by the following theorem: 


1.1 Theorem:* Let Z(p) be ann X n matrix whose elements are Z,,(p), 
| 


< r,s <n, where 

(i) Each Z,,(p) is a rational function 

(ii) Z,.(p) = Z,.(p) (the bar denotes complex conjugate 
(ill) Z,.(p) = Ze-(p) 


(iv) For each set of real constants hy, --- . ht. , the funetion 


ep) = LY Zra(p)krk 
r,s l 


has a non-negative real part whenever Re(p) > 0. 
Then there exists a finite passive network, a 2n-pole, which has the 
impedance matrix Z(p). 


* Presented to the American Mathematical Society April 17, 1948. Abstract 
260. Bulletin of the A.M.S. No 54, July, 1948 
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Conversely, if a finite passive 2n-pole has an impedance matrix Z(p), 
that matrix has the properties (1), (11), (1), (iv 

\ formally identical dual theorem holds for open-circuit: admittance 
matrices ) p 


1.2 A general realizability theorem, applicable to and characterizing 


completely all finite passive networks, whether having impedance ma 


trices or not. is also proved 


1.3 An effort is made to lay a foundation adequate for the realizability 


theory of both active and passive multi-terminal devices. To this end, 
a large part of the paper is devoted to the serutiny of fundamental 


properties ol networks. 


IH. INTRODUCTION AND FOREWORD 


2.0 Network theory provides direct means for associating with an 
electrical network a mathematical description which characterizes the 
behavior of that network. Typically, this results in shifting engineering 
attention from a detailed, possibly quite intricate, electrical structure 
to a mathematical entity which suecinetly describes the relevant be 
havior of that structure. An essential feature of this shift in foeus is 
emphasized by the word ‘relevant’: only those terminals of the net- 
work which are directly relevant to thé problem at hand are considered 
in the mathematical description. Design work can then be done in 
terms of constructs relating explicitly to these accessible terminals, the 
effect of the Internal structure being felt only by implication. 

The physical origins of these mathematical constructs, and the im 
plications of the internal structure upon them, cannot however be en- 
tirely forgotten, for thev have mathematical CONnSeGuUences which are 
not alwavs immediately evident. Until he knows these limitations 
Imposed upon him by the physical nature or the necessary structural 
form of the networks he is designing—a design engineer cannot make 
free use of the mathematical tools that network theory has provided. 

We give the name “realizability theory” to that part of network 
theory which aims at the isolation and understanding of those broad 
limitations upon network performance, 1.e., upon the mathematical 
constructs which describe that) performance—-which are imposed by 
limitations on the network structure. One may also include in’ the 
province of realizability. theory some of the converse questions: the 
study of those structural features common to all networks whose per 
formance is limited in some specified way 


Realizability theory would have little content were it not that “per 
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formance” here must be construed to mean performance as viewed from 


the accessible terminals only. Were all branch currents and node poten- 
tials in a network available to observation, a mathematical statement 
of performance would be equivalent to stating the full system of dif- 
ferential equations governing these quantities, i.e., equivalent to giving 
the detailed network diagram. 


2.1 With a few important exceptions, the converse kind of problem in 
realizability theory does not lead to a strict implication from functional 
limitations to structural features, because the field of equivalent struc 
tures for a specified performance is very broad. Typically, it is only by 
imposing some general a priori limitations on structure that further 
conclusions can be firmly drawn from a functional limitation. In study 
ing this kind of problem one is rapidly led from those basic issues which 
are Clearly part of realizability theory toward general, difficult, and 
usually unsolved problems of network synthesis. One cannot, and should 
not, draw a sharp boundary here, but Nature so far has provided a 
fairly definite one for us, in that most of these problems have proved 
too difficult of solution. 


2.2 The direct realizability problems, the passage from structural prop- 
erties to functional properties, have been somewhat more tractable. 
Here, again, there is no clear dividing line between general realizability 
theory and the sort of design theory in which, for example, one specifies 
a particular filter structure depending on a limited number of param- 
eters and examines the performance of the structure as a function of 
these parameters. There is an extensive literature at or near this latter 
level of generality, most of it relating to filters or filter-like structures 
(e.g., interstage couplers in amplifiers). 

At a more basic level, the limitations on a network’s structure which 
are commonly met in practice are of the following kinus: 

a. Limitations on the kind of elements appearing, e.g., to passive 
networks, networks without coupled coils, networks whose elements 
have specified parasitics, etc; 

b. Limitations on the general form of the network diagram, e.g., to 
ladder or lattice structures, without limitation to a specified number of 
elements or parameters. 

Here the problems are varied and difficult. We survey briefly the 
present status of some of them. 

93 Networks with two accessible terminals, two-poles, are basic in 
network technology. Fortunately, also, two-poles are unique among 
networks in that there is always a simple way to describe their perform- 
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ance, Except for the trivial limiting case of an open circuit, every two- 
pole has a well-defined impedance, Z(p), a function of the complex 


frequeney parameter p, which describes its performance in a way which 


is by now well understood, Dually, except for the limiting case of a 
short circuit, every two-pole has a well-defined admittance function 
Yip). Even the limiting cases are tractable: every open circuit has the 
admittance function Y(p 0 and every short circuit the impedance 
function Z p (). 

ln other words, by exercising his option to speak in terms either of 
Impedance or of admittance, one can always describe the performance 
ola two pole by using a single function of frequency. 

The descriptive simplicity and practical importance of two-poles led 
early to a fairly complete realizability theory for them. In 1924 R. M. 
Foster’ gave a function-theoretic characterization of the impedance 
funetions of finite passive two-poles containing only reactances. The 
corresponding problem for two-poles which are not at all limited as to 
structure, beyond being finite and passive, was solved by O. Brune’ in 
1931. The effeets of various structural limitations have since been 


studied by several writers (cf. Darlington,’ Bott and Duffin 


24 Technology, and the promptings of conscience, have meanwhile 
urged the study of devices with more than two accessible terminals. 
Here, however, Nature has been less kind, in that no uniquely simple 
method is available for describing the performance of such devices as 
viewed from their terminals. 

Indeed, basic network theory has been remiss here, in not even mak- 
ing available a mode of description which is generally applicable 
whether simple or not 

W. Cauer’ showed that, when one admits ideal transformers among 
his network components, it is sufficient to study networks which are 
natural and direct generalizations of two-poles, namely, 2n-poles,* for 
arbitrary values of n. The corresponding natural generalization of the 
impedance function Z p) ob a two pole is the impedance matric of a 
2n-pole: just as one multiplies a scalar current by a scalar impedance to 


get a scalar voltage, one multiplies a veetor current by an impedance 


matrix to get a vector voltage. 

2.41 Not all descriptive difficulties are resolved, however, by consider- 
ing 2n-poles and their impedance or admittance matrices. For the 
moment, a simple example will suffice to show this: the 2X 2-pole which 


consists simply of one pair of short-cireuited terminals and one pair of 


ind also later here 





FORMAL REALIZABILITY THEORY 22] 
open-circuited terminals is a finite passive 2n-pole (n = 2) which has 


neither an impedance matrix nor an admittance matrix. 


2.42 When one eliminates this kind of descriptive difficulty by fixing 
his attention only upon 2n-poles for which an impedance matrix (or, 
dually, an admittance matrix) is available, the general realizability 
problem for finite passive devices is solved. A partial solution, for the 
case n = 2, was published by C. M. Gewertz in 1933. The solution 
(Theorem 1.1) of the problem for a general value of n has been an- 
nounced recently by three authors, independently: Y. Oono,’” in 1946,* 
the present author, in 1948,¢ and M. Bayard,’ in 1949. The problem 
for reactive 2n-poles is much simpler and was solved by Cauer,” in 
1931. 


°) £ 


.) Intermediate between the fairly specific problems of filter theory on 
the one hand and the general realizability theory of multi-terminal 


devices on the other, lies the study of four-poles as transducers. There 


is a considerable literature on the realization of transfer funetions or 
transfer impedances under various structural limitations. The basic 

. . 14 . ‘ ‘ 
and extensive work of Bode™ on active systems belongs also in this 


category since it is avowedly limited to single-loop structures. 


2.6 Beyond the important result that, by sufficiently elaborate circum- 
ventions, one may avoid the use of transformers in the synthesis of any 
two-pole, (Bott and Duffin’’) little in general is known about networks 
which do not have transformers. 


2.7 The present paper is an essay in the realizability theory of devices 
with many accessible terminals. Ideal transformers are admitted as 
network elements; indeed, their use is essential. This fact is indicated 
by the adjective “formal” appearing in the title. 

The availability of ideal transformers makes it possible to exploit 
the simplification noted by Cauer and to consider only networks which 
are 2n-poles in his sense. The aim of the paper, therefore, is to set a 
foundation for realizability theory which is completely general within 
this framework. 


2.71 The first problem is that of deseription. We observed above an 
example—entirely trivial—of a passive four-pole which had neither an 
impedance nor an admittance matrix. Unfortunately, opportunities 


* Date of Japanese publication This reference, and manuse ipt of Oono , 
were sent by Professor Oono in «a personal communication to R. L. Dietzold, who 
showed them to me in December, 1948, while an early draft of the present paper 
was in preparation. The recent (1950) American republication of Oono!® unfor 
tunately omits reference to the original 

+ Cf. footnote to 1.1 
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for this kind of degeneracy become manifold in multi-terminal devices, 
ind some degree of degeneracy is the rule rather than the exception. 
Consider an entirely practical example: that of an amplifier chassis from 
Which the tubes have been removed.* Here the degeneracy is essential 
and intrinsic; it would be highly artificial to regard it as the mere accel 
dent of a limiting cause. True, given any particular degenerate network, 
there is usually evident a method for representing or deseribing its he- 


havior. What is needed, however, is a mode of representation which is 
applicable generally to any 2n pole Without a priori knowledge of its 
structure or peculiar degeneracies 


s> "74> 


2.72 The mode of representation adopted in this paper, embodied in 
the notions of general 2n-pole (Section 4) and linear correspondence 
Section 6), is an obvious one, and so completely general that it solves 
ho problems other than the elemental one for which it was introduced. 
It pro ides a definite mathematical construct: whose properties one can 
discuss with mathematical precision. This is all that we ask of it 

it alizability theory begins and ends with the study of these proper- 
ties. It would be more accurate to sav that the notion of general 2n-pole 
deseribes a particular, but still very large, class of mathematical en- 
tities; realizability theory consists in the study of certain subclasses of 
the whole class of these entities, the particular subclasses being distin- 
guished by special, and to us interesting, properties. 


2.73 Despite its avowed aim at generality, the paper is oriented toward 
the realizability theory of finite passive networks. It ultimately provides 
i proof of 1.1 and indeed a complete characterization of finite passive 
2n-poles, however degenerate. This characterization Is accomplished in 
a sequence of postulates, each one delineating a property of general 


2n-poles, ie., a subelass consisting of all 2n-poles having this property. 


The | - ol 2n poles having all of these properties is then identified 


vith the class of 2n poles obtained from finite Passive networks 


2.74 Ht we have succeeded here in our hope to set an adequate founda- 
tion for the realizability theory of devices with many terminals, it will 
be because of the nature and organization of the postulates themselves 
Phey describe what at present seem to be individually significant prop- 
erties of 2n poles, ot progressively greater specificity, which in the 
aggregate characterize finite passive devices. By eliminating them in 
various combinations one obtains larger classes of objects. Further re- 


tly this « x unple, ind the prac tical need of an ade quate theory for 
iuthor first to study the realizability theory of passive multi 
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search alone will tell whether or not one obtains in this way the kinds 
ot device which sre significant. kor example, one would like general 
realizability. theorems for structures containing vacuum tubes with 
frequency-independent transconductances, vacuum tubes with non 


vanishing transit times, unilateral devices with specified parasitics, ete. 


2.75 Actually, the postulates as we have given them are certainly not 
adequate for such an ambitious program. Exigencies of the presentation 
have dictated a number of condensations and compromises. It is hoped 
that the basic ideas are still evident even if not isolated indiy idually in 
separate and entirely independent postulates. In any event, it is the 
author’s firm belief that the presentation as given is at least illustrative 
of the kind of approach, and the level of mathematical detail, which 
will be needed if one is ever to provide a truly adequate realizability 
theory: a theory which will cover, for example, the broad range of active 


linear systems which present-day technology allows us to consider. 


2.8 Apart from the network theoretic concepts, Which must be evalu- 
ated by their effectiveness in solving problems—an assessment which is 
by no means vet complete—-this paper is strongly marked by an idio 
svneracy of its author: a consistent and insistent use of geometric ideas 
and terminology. This is based on the personal experience that linear 
algebra achieves logical unity and a freedom from encumbering notation 
when viewed in this way. A general reference covering most of the linear 


algebra (geometry) required here is P. R. Halmos’ elegant monograph’. 


2.9 Fora proot solely of 1.1, which has already been three times proved 


in the literature, ' ' this paper provides an apparatus which is too 
cumbersome. There is even a sense in which 1.1 alone provides a charac 
terization of all finite passive devices, for it seems to be generally ac 
cepted that, by the use of ideal transformers, any finite passive net work 
can be represented as a network which has an impedance matrix to 
which is adjoined suitable ideal transformers. Therefore we cannot claim 
that, in using this cumbrous apparatus to characterize all finite passive 
2n-poles Gncluding the degenerate ones), we have offered anything not 
already provided by a simpler proof of 1.1. 

Three things may be said in rebuttal. First, we have already empha 
sized that the apparatus here exhibited was designed for more problems 
than that to which it is here applied. It is presented in the belief that 
it will prove of further use. 

Second, even in the study of passive networks, it has seemed to the 


author helpful to look at the manifold things which are not passive net 
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works. One gets then a clearer view ot the unique position occupied by 
passive devices among all linear systems 

Third, there is a kind of semantic issue here: the assertion that any 
finite passive network (sic) can be put in such a form that 1.1 applies 
seems to this author to give a kind of circular characterization of such 
devices. A characterization which did not itself involve the concept of 

network seems more satisfying. Logically, there is no circle here, but 
this is a fact requiring proof. A careful reading of this paper will show 
that it provides a proof, This particular subtlety does not of itself Justify 
the lengths to which we have gone. It is, however, no longer a subtlety 
if one wishes to consider devices which do not have a representation in 
terms of something non-degenerate to which ideal transformers have 


been added 


2.91 The present Part | of the paper is so organized that at the end of 
Section S& the reader is in possession of all of its principal results and its 
basic ideas. The remaining Sections, 9 through 20, may then be regarded 
as an Appendix contaming the details of proofs. Indeed, Part IL will be 
largely devoted to further details ot proot, though there will be there 
one important idea not mentioned, save casually, in Part [the idea 
of degree for a matrix, 

In Sections 4 through 11, technical paragraphs have been distinguished 
from explanatory or heuristic ones by starring the paragraph numeral. 

Part [1 of the paper contains the bulk of the proof of 1.1. This proof 
is modelled after that of Brune’ for the realizability of two-poles. One 
familiar with the Brune process will probably find Part I] readable 


without extensive reference to Part | 


Let the reader be warned that the Brune process is not a practical 


one for realizing networks because of its critical dependence upon a 
difficult minimization and balancing operation. The same criticism 
ipplies to the generalized Brune process of Part IT. 
Phe Brune process is of theoretical importance because it does realize 
‘twork with the minimum number of reactive elements. These facts 
lL be brought to light in Part IT. 
Phe proots of Oono and Bayard’ are different from ours. That of 


done wiih follows the srune model, 


INTRODUCTION TO PART 1 


3.0 We keep before us first the problem of finding a mathematical de- 
scription applicable to and characterizing the behavior of all finite pas- 
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sive networks. Second, we seek to make mathematically precise those 


ideas which appear to form the basis of general realizability theory. 
Sections 4 through 7 introduce the immediate mathematical machinery 
for this. Section 8 states the fundamental realizability theorem and 
outlines its proof. At this point the reader has had an introduction to 
the results of the paper. The remainder of the paper is then devoted to 
the technical details of proof. Beginning with Section 12, the device of 


starring the technical passages will be dropped. 


3.1 Cauer” distinguished precisely the class of networks called 2n-poles 
from the class of all multi-terminal networks. He also showed that, by 
the use of ideal transformers, any multi-terminal network is equivalent 
to a network which is a 2n-pole (for some n) in his sense. We shall in 
Section 4 define a class of objects to be called general 2n-poles. This 
¢lass includes all electrical networks which are 2n-poles in Cauer’s sense. 
Its definition abstracts the significant properties isolated by Cauer. 

For the study, alone, of finite passive networks, this definition is 
unnecessary, since one can in fact so put the arguments as to deal only 
With 2n-poles which are finite passive networks, and therefore to deal 
only with concepts already defined in Cauer’. The somewhat physical 
notion of a general 2n-pole is a convenient backdrop against which to 
display the important physical properties of finite passive networks, 
and, indeed, of networks in general. Having it available, we use it 
throughout the realizability arguments. 


IV. DEFINITION OF GENERAL 2n-POLE 


1.0* Network theory establishes a correspondence between oriented 
linear graphs and systems of differential equations. With each node of 
the graph is associated a potential £, = £,(4) and with each oriented 
branch acurrent J, = /;(). These potentials and currents are constrained, 
first by Krichoff’s laws, and second by differential equations which de- 
pend upon the nature of the branches but not upon the topology of the 
graph. 
1.01* A finite passive network is one whose graph has the following 
properties: 

(i) There are finitely many nodes, 1, 2, 


(ii) There are finitely many branches, 1, 2, 


* Technical paragraph as explained in Section 2.9] 
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iii) Let the b-th branch begin at node n, and end at n,. Let vp 
/ I. . Then for each b, one ot 


Rolo, R, 
dv» 


holds, where the matrix Ly», is real, symmetric, and semi-defi- 


ile 


Cauer has shown’ how an ideal transformer can be defined as the 
limiting case of a finite passive network. It is indeed no more nor less 
ideal than an open circuit PR, © or Cs 0) or a short circuit (2), 0 


or (, s 


1.02 We seldom deal with networks in the detail which is implicit in 
iii) above. We are usually interested in the external characteristics, so 
to speak, of such networks as viewed from a relatively small number of 
terminals (nodes). These multi-terminal devices, however, we continue 
to incorporate into larger network diagrams. It is usually clear how 
Kirchoff’s laws are to be applied in these cases, and what the differential 
equations of the resulting system are. We are obliged, however, to make 
these matters precise before we can deal intelligently with the most 


yeneral physical properties of networks 


1.1 We have seen the two kinds of constaint that a multi-terminal de 
vice imposes on the branch currents and node voltages in a network in 
which it is incorporated: the topological ones contained in’ WKirchoff’s 
laws and the dynamical ones described by differential equations, Cor- 


respondingly, there nre Two aspects to the concept of general 2n-pole. 


1.11* In its relation to Kirehoff’s laws, a general 2n-pole is indicated as 


an object with pairs of terminals (7, 7,), 1 <r < n. Each terminal 


ean be made a node in an arbitrary finite diagram constructed out of 
network elements and other general 2m-poles, with arbitrary values of 
” huis diagram IS not wan oriented linear graph, SO we have ho basis 
for the use of WKirchotf’s laws. From it, however, we construet an or! 
ented linet graph, called the rdeal qraph ot the diagram, by the follow- 


n) 
l 


rule 


Phe nodes of the ideal graph are those of the original diagram. Every 


explained in Section 2.91, 
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oriented branch of the original diagram is repeated in the ideal graph, 
similarly situated and oriented. Between those nodes which, in the 
original diagram, were the (7, , T.) of a 2n-pole N, is drawn a branch 
8, , called the r-th ideal branch of N, oriented from 7’, to T. . This is 
done for each such terminal pair 

Kirchoff’s laws now apply to this ideal graph. 


$.12* Consider a particular 2n-pole N. Let £, be the potential of 7, , 
kK’ that of 7. . Define 
v(t) = E, — E,. 

Then v,(t) is the voltage across the ideal branch 8, so oriented that 
v(t) > O when T, is positive relative to 1. . Let k(t) represent the cur- 
rent entering 7’, . Then k,(t) I,(t), the current in 8, , so k,(t) is also 
the current leaving . A . This is the force of the notion of ideal branch 
and the fact which distinguishes a network which is a 2n-pole from an 
arbitrary network with 2n terminals. 

4.13 For example, the network at (a) of Fig. 1 is not a 2 & 2 pole because 
its currents are not constrained to meet the ideal branch requirement. 
The addition of ideal transformers in either of the ways shown in (b) 
or (c) of the figure converts it to a 2 & 2 pole. Of course in a circuit in 
which the currents are constrained externally—as they would be, for 
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(Cc) 
Fig. 1—Conversions of a four pole to a 2 * 2 pole 


. Technieal paragr iph “us explained ith Section 2 9] 
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example, when the 2 & 2 pole is driven by separate generators in the two 


external meshes—these transformers can be eliminated. The definition 
ol 2n pole requires however that in every context the ideal branch con- 
cept is valid 

$2* The se ond aspect of the concept ol general 2n-pole is that it Im- 
poses some kind of constraint-——other than that implied by 4.11 and 
Kirchotl’s laws upon the voltages “UCTOSS and currents mm its ideal 


branches, Define the symbols 


as the n-tuples, respectively, of voltages across (T, , 7) and currents 


into 7, 1 <r <n. These are added and multiplied by scalars by the 
usual rules of vector algebra. If » and & represent stmultaneous values 
of voltage and current in the 2n-pole N-—1.e., values satisfving all the 
constraints then we say that N admits the pair |v, 4}. 

We say that N admits v if there is a / such that N admits the pair 
vp, kj}. This & is said to correspond to v. Dually, N admits & if there is a 
v (correspon ling to &) such that N admits [v, 4] 

The constraints imposed by a general 2n-pole N on voltages and cur- 
rents are completely described by the totality of pairs [v, 4] which N 


admits. We shall define a general 2n pole, therefore, as 
i) a collection of m oriented ideal branches, as in 4.11, and 


Hh) a list of pairs fe, A] of voltages and currents admitted in these 
branches 
Hereafter we shall usually drop the adjective “general.” 
1.21 The definition we have just given is, in a way, a postulational 
form of an n-dimensional Thevenin’s theorem. It postulates that a 
2n-pole is a thingt which, as far as the outside world is concerned, can 
be represented hy a collection ol two poles 8. .1 <9 Son, among which 
there exists a complicated agreement as to what currents and voltages 
will be admitted 


1.22 The passive networks (b) and (¢) of Fig. 1 define 2 X 2 poles, be- 


cause they satisfy 2.01 and clearly permit a complete specification of 
the admissible pairs [e, &]. Any equivalent network would also specity 


iragraph as explained in Section 2.9] 


this level of generality, any thing 
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the same 2 X 2 pole, because—by its very equivalence—it would admit 
the same pairs. The closest association we can make between a 2n-pole 
and a network, then, is to identify the 2n-pole with an equivalence 
Class of networks. 


1.23 The completely symmetric role played by voltages and currents 


in this definition of general 2n-pole will make it possible to take early 


advantage of the well-known duality principle of network theory. We 
shall do so freely. 


1.3* We shall call a 2n-pole physically realizable if its admissible pairs 
lv, A] are the solutions of a system of differential equations obtained 
from a finite passive network, admitting the limiting elements: ideal 
transformers, open circuits, and short circuits. 


V. PHYSICAL PROPERTIES OF NETWORKS 


5.0 There are clearly a great many properties of finite passive networks 
which are not vet possessed by the general 2n-poles now introduced. It 


is instructive to examine these properties physically. 

5.1 In the first place, the dynamical constraints (a), (b), and (¢) of 4.01 
are expressed by linear, time invariant, differential equations. Accord- 
ingly, the 2n-poles of network theory are: 

5.11 Linear, in that the class of admissible pairs [v, £] is a linear space; 


5.12 Time invariant, admitting with each = |[ve(t), A(t) also all 
lo(t + 7), A(t + 7)| for aribtrary r. 

5.2 In the second place, a physical network N cannot predict the future, 
i.e., It Cannot respond before it is excited. This can be formalized in 
terms of the pairs [{v, |] admitted by N, but to do so would require some 


digression. The reasons will be seen under 5.7 below. 


5.8 We have already mentioned the constraints imposed on voltages 
and currents in a network by the topology of the network, through the 
medium of WKirchoff’s laws. These constraints have three important 
properties: 

5.31 They are workless, since they are imposed by resistanceless 
connections, leakless nodes, and, in the formal theory, by ideal 
transformers. 

5.32 Though it seems scarcely necessary to say it, they are the only 
workless constraints. All other constraints are dynamical and have 


powers or energies associated with them. 


* Technical paragraph as explained in Section 2.91 
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5.33 They are frequeney independent, that is, holonomic in the sense 

of dynamics 

5.4 The workless and the dynamical constraints In a physical network 
all defined by relations with real coefficients. The space of admissible 

pairs is then a real linear space 


».5 The positivities specified in 4.01 are characteristic of passive sys- 
tems. They correspond to the fact that the power dissipation and the 


stored energies are all positive. 


5.6 By definition, finite passive networks contain finitely many lumped 
elements Correspondingly, their resonances and anti-resonances are 
finite in number 


5.7 We are accumstomed to dealing with networks which have, in addi- 
tion to the properties listed above, a kind of non-degeneracy, in that 
the list of admissible pairs [v, ’] satisfies: 


5.71 At least one of v or & can be specified arbitrarily—any real function 


is admitted; 


5.72 When the free number of [v, /] is specified, the other is uniquely 
determined 
For these non-degenerate networks, the property 5.2 above is easily 


formalized: if, sav, & is determined by v, then 


v(t) for 
irmpolie s 


f 


’ 


where [e', &'] are admissible pairs, 2 1,2. The general statement of 5.2 
In) olves this condition and some discussion ot the v’s for which N ad- 


mits {[r, O}, and the dual notions 


5.8 The reason for speaking in terms of pairs |v, 4], instead of in terms 
f ‘cause’ and “effeet,” or “impulse” and “response,” is hinted at by 


5.7 above. For the tacit implications of the cause and effect language 
completely obseure the fact that 5.71 and 5.72 are properties which are 


not automatically possessed by electrical networks. In fact, the simple 


four-pole of 2.41-— a pair of unconnected terminals 7) , 71 , and a pair 


of shorted terminals 7 , 7-—has neither property, vet it is a perfectly 


good linear time invariant four pole. Its admissible pairs are 
0), (0, ke)], 


where v and ke are arbitrary real functions of the time 
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VI. LINEAR CORRESPONDENCES 


6.0 In developing the formal properties of 2n-poles which are equivalent 
to the physical ones just listed, it would be instructive to adjoin re- 
quirements piecemeal, much in the order given in Section 5. Space does 
not permit us full enjovment of this luxury, but the reader will find a 
rough parallel between Section 5 and the developments of this Section 
and Section 7. 


6.1 It is well known that linear time invariant systems are best studied 
by the tools of Fourier or Laplace analysis. We make this fact the basis 
of our first step in characterizing physically realizable 2n-poles simply 
by phrasing our whole discussion in the frequency language. The con- 
tent of the following paragraph will be obvious enough, but it does de- 
fine terms to be used later. 


6.11* Let v and k, without underscores, represent n-tuples of complex 
numbers: 


(1) 


») 


(2 
These are to be manipulated by the rules of vector algebra. Let p be a 
complex number. We shall say that a 2n-pole N admits the pair [v, 4] 
at frequency p, if in the sense of 4.2 N admits the pair [v, 4} (with under- 

scores) where v has components 
v(t) = Re(v,e"'), oe cae (3) 

and & has components 

k(t) Re(k,e”'), ne ae (4) 
Also analogously to 4.2, we say that N admits uv at frequency p if 
there is a k such that N admits [v, 4] at frequency p, and that this / 


corresponds to v (at frequency p). Similarly, N admits k at frequency p 


if there is a (corresponding) v such that N admits |v, k| (at p). 


6.12* Let V denote the aggregate of all n-tuples (1), and K the aggre- 


gate of all n-tuples (2). These are then complex linear spaces. 

6.2* As our first step toward characterizing realizable 2n-poles, let us 
consider a linear corre spond nee L between V and K described by the 
postulates: 

Pl. There is a set [, of complex numbers and for each peD, a list 


L(p) of pairs [v, kh], veV, keK. 


* Technical paragraph as explained in Section 2.91. 
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eL(p) and [e, k|eL(p), then 
yi + go. ay + aol |eL(p 

for any complex numbers a, , a 
6.21* Given such a linear correspondence LL, we can always describe a 
Qn pole N, by: 

N, admits |v, A] at frequency pif and only it lv, h:\eL(p ; 

Phat is, we can alWays interpret the pairs [v, 4] generated by (3) and 
1) from the |v, kjeL(p), tor each peP, , as the voltages across and 
currents in a set of nm ideal branches. We call N; the 2n-pole associated 
with 1 


4) yo 


22* We call , the frequency domain of L (or of N;,). 


6.23 From here on, the words ‘*2n-pole” can with some strain be re- 


garded as suggestive but unnecessary. We in fact deal with linear corre- 
spondences having properties as vet unspecified and shall show how 


physical networks can be constructed which admit the pairs [v, AjeL(p). 


\ctually we use freely the concept of general 2n-pole and thereby avoid 


some elaborate circumlocutions 


6.24* We identify two correspondences L, and Le as being the same it 
1) their frequency domains differ only by a finite set, and (it) for each 
p where both are defined the lists Ly(p and L(p are the same 


6.3 The simplest linear correspondences are those generated by ma 
trices. For example, let Z(p) be ann X nm matrix with, say, elements 
Z,.(p) Which are rational functions of py l <r, s <n. Let I, consist of 
all the values of pat which Z(p) is defined. For pe?, , define L(p) as 
the class of all pairs 


obtained by letting / range ovel K, where for each oe vis defined by the 


matrix equation 
Z p i. i) 


Phis kind of matrix equation will be used throughout to symbolize the 
nm component equations 


The list of pairs Lip) defined by (5) clearly satisfies Pi and F2. it 
can therefore be used to define a 2n-pole N, . It is easy to see that N;, 


is explained in Seetion 2.91 
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in fact is non-degenerate in a sense similar to that of 5.7, for the current 
amplitudes / can be specified arbitrarily, and the resulting voltage 
amplitudes v are then fixed by / and p, by (6). 

Z(p) is called the impedance matrix of the 2n-pole N, . It is also 
sometimes called the open-circuit) impedance matrix, because each 


Z,.(p) 1s, by (7), the voltage amplitude across (7', , 7) when the current 


amplitudes at all terminals save (7, , 7.) are zero—i.e., when all pairs 
save the s-th are on open circuit, 
6.31 Dually, the pairs 

[v, V(p)e} 
defined by an admittance matrix ¥(p) as v ranges over V define a linear 
time invariant 2n-pole which is non-degenerate. 


VII. WORK AND ENERGY 


7.0* A linear correspondence satisfving Pl and P2 is something which 
abstracts the properties of linearity and time invariance. Most of the 
remaining properties of physical networks involve the mention of work 
or energy. These concepts enter our picture by way of the scalar product 
v, k) between a voltage n-tuple (1) and a current n-tuple (2), of 6.11. 
This scalar product is defined by 


(v, k) i 


l 


7.01 If p = tw, one easily calculates from (3) and (4) of 6.11 that 


2 Re(v, k) = lim ; =f ipo (tk, ofa 


roe 


That is, when p = tw, the real part of 2(v, /) measures the average total 
power dissipated by the system of currents /,(f) against the driving 
voltages v,(t). 

When p is not a pure imaginary, the interpretation of the sealar 
product (v, k) is not so clearly physical as this. The reader will ulti- 
mately observe that our significant statements about such products can 
all be reduced to statements applicable when p iw, 1.e., When the 


power interpretation is valid. 


7.1* An important concept in what follows is that of the annihilator of 
linear manifold (Halmos’, par. 16). Let V; © V be a linear manifold. 


* Technical paragraph as explained in Section 2.91 
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innihilator (V is the set of all k such that 
veV implies me. () 


isa linear manifold in K 
Dually, given K, © K, (Ky) is the linear manifold of all veV such 


thicat 
leeK, implies . (). 
| hye 


york of the idea of orthogonality 


hilator cor cept is the analog hh oul general reometric lrame 


It clearly suggests a connection with 
Vorkless Conustraimets. 


»* 


F pe - Phi complex conjugate olf an vt tuple / or | is defined in) the 


ObDVIOUS WaVvi dl 


Phis conjugation operation clearly has the properties 


ae -F ae 7 b# 


where aq and 6 are scalars and € and 9 are (consistently) elements of V 
or K. Furthermore, at once trom (1) of 7.0, 


21* A linear manifold will be called real if it contains, with any 


n-tuple also the conjugate of that n-tuple 


2* A real manifold is spanned by real n-tuples. This will be proved 
\ppendix, Section 2O 


Phe annihilator of a real manifold is real. For let Ky be real and 
be real nw tuples which span K, . Then i ve(K,) every 


v, f 0), 


and conversely, But then also 


is explained in Seetion 2.91 
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7.3* Given a linear correspondence L, we make several definitions: 
V.(p) is the set of all veV such that there is a & with [v, kleL(p 
K,(p) is the set of all keK such that there is a v with |v, kleL(p). 
Vio(p) is the set of veV,(p) such that 

[v, OleL(p). 
K jo(p) is the set of keK,(p) such that 
[O, kleL(p). 


7.31* The postulate P2 implies that for each peP,, Vip), Ki(p), 
Vio(p) and Ky,o(p) are all linear manifolds. 

7.32 V.(p), for example, is the set of veV such that N; admits v at 
frequency p. 

7.4* We now postulate 


P3. There exist fixed linear manifolds V, ¢ V, K, € K such that 
(A) For every pel’, , Vi(p) = Vz (Kyo(p) )° 
(1) For every pel, , Ki(p K, (Violp 


7.41* We may henceforth write Vio , Kio , for Vio(p), Kio(p), knowing 
that, under P3 


Vio = (Kz), 
Kio (V,). 


7.42 Linear correspondences satisfying P3 abstract the properties men 
tioned in 5.3. The equalities V,(p) V.,K.(p) = K, guarantee the 
frequency-independence of the workless constraints. The equalities 
Vilp) = (Kin(p)), Kilp) (Vio(p)) in a sense guarantee that the 
only constraints imposed upon admissible currents and voltages (as 
opposed to constraints relating currents and voltages) are those which 


arise from open or short circuits, 1.e., are workless 


7.43. An illustrative consequence of P3, for example, is that if L satisfies 
P3 and if N, is such that all of the current amplitudes can be specified 
arbitrarily, then indeed the voltages are determined by the currents 


This will appear as a consequence of 8.1. It is a very general theorem 


about networks of a kind that this author, at least, has not heretofore 


encountered. 


* Technical paragraph as explained in Section 2.91. 
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Continuing toward realizabilitv, we introduce 


Pi. Tf pel, , then pel. If |v, kleL(p), then [b, kleL(p). 


This postulate embodies most of the reality properties of networks. 


It has as an immediate consequence the 
7.9O1* Lemma: If L satisfies Pl, P2, P38, and P4, then all of 


Proof: By PA, veVi(p V, implies veV,(p) V.. Hence V;, is 
real. Then K, (V,) is real, and dually. 


7.6* The three remaining postulates on L refer to scalar products, 
They are concerned with the energy questions related to passivity, 
rather than with the workless constraint questions. 

PS. Uf [u, yleL(p) and [v, kleL(p), and if 


\) w and v are real, or if 


(| J and h are real, 


(u, k) (v, 7). 


7.61 This is the property which provides the reciprocity law. In its 
presence, the relations in P3 may be weakened to 

Vilp Vi D (Kiv(p 

K/(p ms a (Vio(p) 

fact will appear as a consequence of the lemma of Section 12. 
Lemma: A consequence of P2 and P8(A) is that if 
kelp), 

for any ueV,, 


kor bv PZ we have that 
kes] (O, ky kajeL(p), 
hence | loeK Then however, by P3 \), ueV, implies ue(K, 


Pechnieal paragraph as explained in Section 2.91, 
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that 
O = (u, ky — ke) = (u, ky) 


Q.E.D. A dual result follows from P3(1). 
7.71* The result of 7 


fixed by v alone when we know that [v, Aje(p). This means that, given 


.7 above means that the scalar product (v, /) is 


veV, , there is a unique function F (p) defined for pel, by 


(v, k) 


where [v, kjeL(p). Dually, 
J,(p) 
is defined for each fixed keK, . 
7.72* (P6.) The complement of IT, is finite and 
(1) For each veV, , F’,(p) is rationai 
(A) For each keKy, , J,(p) is rational. 
7.73* (P7.) (A) Re(p) > 0 implies Re(F,(p)) > 0 


(1) Re(p) > 0 implies Re(J/;(p)) > 0. 


VIEL. THE FUNDAMENTAL REALIZABILITY THEOREM 


8.0* We can now state our fundamental realizability theorem: If a 
linear correspondence L satisfies Pl, --- , P7, the associated 2n-pole 
N, is physically realizable. Conversely, given a physically realizable 
2n-pole N, the associated linear correspondence satisfies Pl, --- , P7. 

8.0L Actually, the postulates Pl, --- , P7 are not unique nor even en- 
tirely independent. Many changes may be rung on them. We indicated 
one above. At the expense of apparent asymmetry, the (A) or (1) por- 
tions, in various combinations, can be deleted or weakened. We shall 
not pursue this subject further at this point, but must come back to it 
in Section 12. 

8.02 We close this Section by outlining the proof of 8.0. The details are 
then contained in the remainder of the paper. 

8.03 The proof that Pl through P7 are necessary for physical realiza- 
bility will be a direct one: it will be shown that, considered indiy idually, 
each network branch and each ideal transformer satisfies the postulates. 


* Technical paragraph as explained in Section 2.9] 
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By an application of Kron’s method (described by Synge’’), it will then 
be shown that the Imposition of Kirchotf’s laws preserves the postulates. 
This work is most efficiently performed after the full machinery of the 


sufficiency proots Is available, and will be done in Section 10. 


8.04 The sufficiency of Pl through P7 can be deduced—and we will do 
<0~-from the lemmas to be quoted below. \part from Section 19 on 
necessity, the remainder of the paper Is devoted to the proots of these 
lemmas 

S.1* Lemma: WL is a linear correspondence satisfying Pl, P2, P38, and 


P4. then there exists a fixed real nonsingular matrix W such that 
S.11 The list Ly (p) of all pairst 
W'v, W’ 


where |p, kleL(p), describes a linear correspondence Ly satisfying 
Pi, P2, P3, and P4. 


S12 The 2n pole N N,; " associated with Ly consists of 
Some number r of open-circuited terminal pairs (7; , 
ee 8 
ii) Some number s of short-circuited terminal pairs (7', 
PF, 
iil) A set of m r s terminal pairs (741, T'41), °°° 


’ 


Tesay Aven): 
8.13 Either m = O, or the terminal pairs in (iii) are those of a 2m-pole N, 
which has a nonsingular impedance matrix Z,(p). 


This lemma, and the following, will be proved in 13.2. 


8.2* Lemma: If L satisfies P5, P6, and P7, then Z;(p) is a positive 


realf matrix, that is, Z,(p) satisfies (1), --- , (iv) of 1.1. 


3* Lemma: If a 2m-pole N, has a positive real impedance matrix, then 
l 


oe 
N, S physically realizable 
This is the sufficienev half of the matrix realizability theorem 1.1. 


Part 1] will be devoted to lis prool 


S.4* Lemma: If Ny is physically realizable, then N can be constructed 


from it by the use of ideal transformers. 
This is Cauer’s Pransformatiou Theorem about which we shall Say 


more in Section & 


raph as ¢ xpl ined in Seetion 2.9] 
respectively the recipros al and the transpose of 


nolog by now traditional 
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8.5* The sufficiency half of 8.0 is now clear. By 8.2 and 8.3, N, is physi- 
cally realizable. Clearly then Ny is, simply by the adjunction of the 
necessary open and short circuits. Finally N is by Cauer’s theorem, 8.4. 
8.6* We can see now how to prove the necessity of positive reality for 
the realizability of a positive real matrix Z(p). This is the necessity half 
of the matrix theorem 1.1. Let Z(p) be the matrix of a realizable N. 
Then N has an associated linear correspondence L satisfying Pl, --- , P7, 
by the necessity half of 8.0. The pairs of L are the pairs 


[Z(p)k, kl] 


generated as / ranges over all n-tuples. By definition, then, the pairs of 
L w are 
[WZ (p)k, Wk]. 


As k ranges over all n-tuples, the nonsingularity of W implies that 
Wk does also. Let U W". Then the pairs above are the same as 


[(UZ(p)U'k, k] 


as kh ranges over all n-tuples. Hence Lw has the impedance matrix 


rr , yj . 
UZ(p)U", where U = Wis real and nonsingular. Because Ly has an 


impedance matrix, r = O in 8.12. 

Now by 8.1 and 8.2, Z,(p) is positive real and the matrix ('Z(p)U’ 
of Ly is Just Z,(p) bordered by s rows and columns of zeros. It is then 
easy to see that (Z(p)U” is positive real, and finally also that Z(p) is. 
These last two facts will be proved formally in Section 16 


IX. CAUER’S TRANSFORMATION THEOREM 
9.0 Cauer’s transformation theorem’ is the cornerstone of formal reali 
zability theory. In one form, the theorem reads: 
9.1* Let Z(p) be the impedance matrix of a physically realizable 2n-pole 
N. Let U be a real, constant, nonsingular matrix. Then 

UZ(p)l”" (1) 
is again the impedance matrix of a physically realizable 2n-pole, Ny . 
N, can be constructed from N by the use of ideal transformers 
9.2* A superficial generalization of this theorem can be obtained at once 
from Cauer’s proof. It asserts that if N is physically realizable and is 
described by the linear correspondence L, then there is a physically 
realizable 2n-pole Ny, obtainable from N by the use of ideal trans- 


* Technical paragraph as explained in Seetion 2.91 
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lormers, which Is described by the linear correspondence L " whose pairs 


at each p are the pairs 
[Wo'e, WK], 


where |v, kleL(p 

We refer to Cauer’ for the proof. [t is straightforward. 
4.21 We shall use the second form (9.2) of Cauer’s theorem in our 
realization process. Notice that it is in a sense a “‘physical”’ theorem, 
about the way one physical network is related to another. It is used in 
this wav: we shall always solve a realizability problem by finding some 
network N which is easily realized, and then a W such that Nw , which 
is now realizable, proy ides a solution to the given problem. 
4 .22* We shall call the 2n-pole Nw a Cauer equivalent of N. 
0.3 Although Cauer’s theorem will be applied, in a sense, only a posteriori, 
its effect is fundamental. For it implies that formal physical realizability 
is a property of matrices which is invariant under the operation (1) 
or a property ol correspondences which is invariant under (2). There is 
an extensive classical literature on the properties of matrices Invariant 
under operations like that of (1), and the effeet of Cauer’s theorem is to 
make these results all available to formal realizability theory. 


O31* [tis worth observing here that we are already well set up to use 
Cauer’s theorem: 

Lemma: Ut Lis a linear correspondence satisfying Pl, --- , P7, then 
the correspondence Ly of 9.2 also satisfies Pl, «++ , P7. 

Proof: Let M Ly. Pl and P2 for M are obvious, with Py ae 


By definition of JV, 
Viulp Wo'V.(p wv, 
K w(p W’K,(p W’K, 
V wolp Wo! Vio(p WV 
K wo(p W’'Kio(p) = WK 


where WS for a manifold S consists of all n-tuples W '» where 


Hence in P83, 


V ulp) Vu W ‘V, 
 e p Bx W’K, 


for fixed manitolds Vy, Ky as defined. 


cal paragraph as explained in Seetion 2.91 
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Now if veV,o , then 


for every keK, = (Vio). Then, however, by direct calculation from 


Section 7.0. 
(Woy. W*k) (). 


where W* is the adjoint, Le. transposed conjugate matrix of HW. But 
because Wis real, W* = WY’. Hence if veVio , then 


(Ww ty, I) 0 
for every keW’K, = Ky. Hence 
K y = (W 'V, CV wolp 


By this and its dual, P3 is completed for J. 

The reamining postulates for / follow from those for L by the simple 
equality 

(vk) = (We, W’k) 

already established, combined with TP yy Fr; 
9.32 For fixed Z(p), the matrices (1), as (° ranges over a group, form an 
equivalence class. Classical matrix theory treats of such equivalence 
Classes. This author’s predilection is to regard this theory from a geo- 
metrical point of view. In part this prejudice may be justified by the 
ease with which that slightly more general object, a linear corre- 
spondence, can be treated by geometrical methods. In any event we shall 
begin our program of proofs with a brief introduction to the geometrical 


approach. 


X. GEOMETRICAL PRELIMINARIES 


10.0* We now wish to consider V and K as complex n-dimensional 
linear spacest respectively of voltage vectors v and current vectors k, 
The distinction here is in point of view. A vector v is regarded as an 
absolute geometrical object; an n-tuple [v| la, , «** , a,| is regarded 


| 
as a set of coordinates for the vector v, relative to some coordinate basis. 


Given a fixed coordinate basis, there is a one-to-one correspondence 

between vectors v and the n-tuples [ve] which represent them in that 

basis, a correspondence which preserves the operations of vector algebra. 
* Technical paragraph as explained in Section 2.91 


t For nN reference concerning the ideas in this section, see Halmos?® Ch ipote rs 


I and II 
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10.01 The effeet of attaching a geometric identity to vectors, rather than 
tow tuple sis to make it possible to choose coordinate bases freely and 
as convenient, Without elaborate constructions or even interpretations. 
We « ! discuss properties olon tuples (and other objects, eg. 
matrices) Which are Invariant under the kind of operations exemplified 
hy 1) and (2) of Section 9 as properties of a single qeomelric object, 
rather than as properties shared by an extensive class of concrete ob- 
jects which are converted into each other by the group of operations, 
10.1 This change in point ob view need not change formally anything we 
have said to date; it simply erects a conceptual superstructure, or pro 
vides a conceptual foundation, depending on the reader’s personal 
attitucl 

We shall support this statement by going through the important ideas 
of Sections 4,6, and 7 and examining their geometrical meanings or 
counterparts. It is convenient to consider first and at some length the 
notions of scalar product and complex conjugate. The geometric struc 
ture will then be complete enough to permit a rapid survey of the 


remaining ideas 


10.11* The geometrical counterpart ot the sealar product introduced in 
7.0 is a numerically valued function o o(v, hk) of two vector variables. 
Its first argument » ranges over V and its second argument / ranges over 
K. ‘The function o(7, /) is linear in ve and conjugate linear in k: 

Lol ag u. | t ba(v, k) 


(1) 
if dal(v, k) + bale, &). 


We denote this funetion a(v, k) by the simple bracket notation (v, /). 


10.12 With this sealar product, the geometry of V and K is that of a 
space K and the space K* V of conjugate linear functionals over K. 
Phis is analogous to the real geometry of space and conjugate space 
discussed at length in Halmos®. In fact, in the introduction to Chapter 
HT of Halmos’, the modifications introduced by the conjugate linearity 


ol over K are treated in detail 


10.13* Because of its importance, we quote here a paraphrase of the 


results covered in Halmos’, par. 12 
f fie) is anv numerically valued homogeneous linear funetion of 


there Is a unique vector 4 eK such that 


} 


a , 


is EXP uned in Seetiron 2.91, 
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(ii) If g(k) is any numerically values homogeneous conjugate-linear 
function of eK (i.e., if g(k) is linear in k) then there is a unique v,eV 
such that 

g(k) = (, k) 
for all keK. 


10.2* The annihilator (V;) of a manifold V; C V is, as in 7.1, the set of 
all keK such that 


veV, implies (v, k) = 0. 


10.21* It is shown in Halmos® that to each basis v,--- , 2" in V there 
exists a unique dual basis kh’, «++, k" in K such that 


(v’. k*) 
where 6,, 1s the Kronecker symbol: 6, Oifr ¥ s,4,, 
10.22 If 


ere 


[k] = [bi --+ , bal 


are the n-tuples representing v and /& relative to a pair of dual bases, 
then it is easily computed from (1) and (2) that 


(vy. k) = > a,b,. 1) 


r=] 


Therefore the concrete scalar product of 7.0 is indeed the geometric 
scalar product here considered, when we restrict our pairs of bases in 
V and K always to be dual in the sense of (2). 
10.23* We shall use the words ‘coordinate frame” or simply ‘“‘frame”’ 
to denote a pair of dual bases in V and K. Any basis in V (or K) specifies 
a frame by the uniqueness result quoted above. 


10.24 We shall henceforth deal always with coordinate frames, in fact, 
ultimately, real coordinate frames, rather than arbitrary pairs of bases. 
This means in classical language that we are considering as “geometrical 
properties” all properties which are preserved under the group. of 
linear transformations which leave the bilinear form (4) invariant 
The properties related to physical realizability will turn out to be 
invariant only under the subgroup of real linear transformations pre- 
serving (4). 


* Technical paragraph as explained in Section 2.91 
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10.3* Conjugation is an operation which to each veV associates a 


vector ® uniquely determined by v with the properties 


au + th au -t bi, 
h are any complex numbers and 4d, b their conjugates. 


ven any such conjugation operation in’ V, and given any 


eK, detine a function gle) by 
Gilt (v, (6) 
for veV. Then gy(v) is linear inv, by (5) above and (1) of 10.11. There- 
fore, by 10.13, there is a unique vector keK such that 
Qala (v, k). (7) 
10.32* Directly from (1) of 10.11 and (6) above, if 7 = ak + bf, then 
aokv) + bajo). 
From (7), therefore 
a(v, k) + hiv, 0) 
eV. Comparing this with (1) of 10.11, we see that 


ak + bf. 


Phe second item of (5) above then holds for vectors keK. 
Phat / ; follows easily: We have from (6) and (7), written for the 


vector hk, that 
v, k) v,k 
We also have, by writing (6) and (7) for vectors @ and k that 
Paking complex conjugates of these two numbers, and using v 
from (5), we have 
(10) 
Phen (9) and (10), which hold for all veV, identify & and k by 10.13. 


1lO.34* We have now showed in (5), (8) and (10) that this complex 


conjugate satisfies the formal properties of the conjugate for n-tuples 


se | ) 


introduced in 7.2 


Pech paragraph as explained in Section 2.9] 
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10.35. The abstract scalar product of 10.11 turned out in the end to be 
no more than the concrete one of 7.0 when we restrict our attention to 
n-tuples derived from vectors by the use of coordinate frames. In a 
similar way, it is not hard to show that there always exists a coordinate 
frame in which the abstract conjugation now introduced has the form 
of 7.2. This will be done in the Appendix (20.2). 

10.36* Our need for writing out the components of vectors has now 
almost vanished. Henceforth we shall use subscripts to denote particular 
vectors, e.g. vu, rather than components. 

10.4* A vector will be called real if it is equal to its own conjugate. 
A manifold will be called real if it contains with each vector also the 
conjugate of that vector. V and K are then real. A basis will be called 
real if it is made up of real vectors, and a frame will be called real if its 
bases are real. Any frame in terms of which our conjugation operation 
takes the form of 7.2 is real by definition because its basis vectors in 
that frame have components which are 0 or 1. The vector O is real, 
similarly. 


10.41* The basis dual to a real basis is real, for if 


then by (10) of 10.3 and the hypothesis that v, = 0, , we have 


= 6, 


so the k, satisfy the same equations as the hk. The uniqueness of the 


basis dual to 7; , --- , v, then proves that k, = ky, 1 <s <n. 


10.42* Any vector v can be written 


10.5* It is shown in Halmos", par. 34, that if veV, keK are represented 


|, , [A], in some other frame, 


by [v], [Kk] in some coordinate frame, and by |/ 
then there is a nonsingular matrix [W], which (a) depends only upon the 


- Technical paragraph us explained in Section 2.91, 
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two frames, and (b) relates these n-tuples as follows: 


(11) 
[K] (W}*{A]. 


jl 


It is easy to show that if [W] has real elements, so that [W]* Ls¢ 
then the two frames involved above are either both real, or else neither is 
real. Also, conversely, if both frames are real, then necessarily the [|W] of 
11) has real elements and [Wy]" [Wy’. 
10.6* Some further important geometrical notions must be mentioned 
before we proceed. 

If V, and V2 are disjoint linear manifolds in V—1.e. linear manifolds 
having in common only the single vector 0-—we write 


V; t V, 


for the linear manifold consisting of all vectors v ¥) + vw, Where 
eV; , 1 1, 2. The cirele around the plus sign is used to denote the 
disjointness ol V; and V, , 


It is shown in Halmos’, par. 19, that if 


V V; t V, 


K K, t K, 9 


where K, V.), Ko (V,) and the dimension of K, is equal to that of 
V;,2 |, 2. We eall (13) the decomposition dual to (12). We some- 
times write K, Vt to denote the K, dual to V, in the decomposition 
13). It is shown in Halmos’, loc. cit., that there exists a basis v,,°°° , 

in V and its dual k,,--- ,k, in K such that, if r is the dimension 


of Vi, 
is a basis for V, 
is a basis for V2 
_k, is a basis for K, 
,k, is a basis for Ke. 


Furthermore, if 7, v, is any basis in V satisfying the first halt 


of (14), its dual basis satisfies the second half, and dually. 


We shall show in the Appendix that if any one of V,, Ve, Ki, or 


Pechmesl paragr iph is ¢ xpl ined in Section 2.91 
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K, is real, then they all are, and that in this case the bases (14) can be 
chosen to be real. 
Similar considerations apply to decompositions into more summands: 
if 
V . V; us V2 


then 
K . K, $ K, u 
where 


* 


vi =K,=N Vv} = 


pytt 


XI. GEOMETRICAL CORRESPONDENCES 


11.0 With the geometry of V and K now in hand, we consider the 
geometric aspects of our network theoretic concepts, 

The definition in Section 4 of general 2n-pole describes a concrete 
thing and stands unaltered in our geometric view. The definitions in 
6.11 of the terminology typified by ‘““N admits |v, k] at frequeney p”’ 
are unchanged except that we should now explicitly indicate that we are 
discussing concrete n-tuples of complex numbers by placing brackets 
around the vector symbols, thus: [v], [A]. In other words, a 2n-pole is 


described by a concrete relation between n tuples. 


11.1* All of the postulates Pl, --- , P7 are stated in a language which 
now has been given an absolute geometric meaning. In this meaning, 
P| and 2 describe Pa | geometrical linear corre spond nee between vectors 
veV and eK. This is the geometric counterpart of the concrete 


notion of a linear correspondence between n-tuples. 


11.11 An impedance matrix, as in 6.3, deseribes a particularly tightly 
knit linear correspondence, namely a linear function from K to V. 
The geometrical counterpart is an impedance operator which for each 
p is by definition a linear homogeneous function which assigns to each 
vector keK a unique v Z(p)keV. That is: an operator is a functional 
relationship between vectors and as such has a geometric identity. 
11.12 It is easy to provet that, given an impedance operator Z(p), 
and given any coordinate bases in V and K respectively, there is a 
matrix [Z(p)|, with elements Z,,(p), 1 < r,s <n, such that relative to 


these bases the coordinates kh, of a vector k and the coordinates », of 


v Z(p)k are related by (7) of 6.3. We call [Z(p)| the matrix of Z Pp) 


* Technical paragraph as explained in Seetion 2.91 
t Cf. Halmos?’, par. 26 
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relative to the given pair of bases. A strong analog of this observation is 


contained in the following lemma 


11.13* Lemma: (i) Let L be a geometrical linear correspondence. Fix 
any real coordinate frame and let [Z| be the linear correspondence 


whose paired 7 tuples are 
((e, (ell, 


where 


. lel. p 


ii) Alternatively, let [Z| be a (concrete) linear correspondence be 
tween n-tuples. Interpret the n-tuples related by [Z| as representing 
vectors im some real coordinate frame. Let L be the geometrical cor- 
respondence Whose pairs, expressed as n-tuples in this frame, are those 
of the concrete correspondence [L}. 

In either ease, (1) OFr (il), the vyeometric correspondence L satisfies the 
geometric postulates PL, ---, P7 if and only if the concrete corre- 
spondence {L| satisfies the concrete forms of these postulates. 

The proof of this lemma consists essentially in reading the postulates 
carefully. We shall not reproduce it 
11.2 Our position is now this: We have on the one hand geometrical 
objeets, veetors v, kk, operators Z(p), ¥(p), and geometrical correspond- 
ences L. On the other hand, we have concrete n-tuples [v], [A], matrices 


Zip)\, |\¥Cp)|, and linear correspondences [L]. Given any pair of bases 


in Vand K, in particular, given any coordinate frame, each geometric 


object: generates a corresponding concrete object which represents it 
relative to those bases or that frame. Conversely, given a concrete ob- 
jeet [E], we can choose a frame in Vand K and find that geometric object 
Whose coordinates in the chosen frame are given by (é]. 
11.21* The concrete object, linear correspondence, defines a linear time 
invariant 2n-pole by 6.2L. To complete the picture, we might say that a 
geometrical correspondence L. defines a Cauer class of 2n-poles. 
11.22* This terminology is motivated by the following observation: 
if {7} and [/); are linear correspondences representing L in two distinet 
real frames, then there exists a real nonsingular matrix [W] relating the 
v], (Ae fe[L)(p 
and the 
= [A|iJe[L],(p 


h as explained in Section 2.91, 
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by the formulas of 10.5. This means that [] and [1], are related like the 
[L) and [Ly] of 9.2. The 2n-pole associated with [L], therefore is a Cauer 
equivalent of that associated with [ZL]. 


11.23 The observation of 11.22, combined with (ii) of 11.13, gives an 
alternative proof of 9.31. This proof is deceptively free of calculation, 
but of course the calculations are concealed in the extensive geometrical 
developments of Section 10, many of which are there offered on faith. 


XII. THE FUNDAMENTAL LEMMA 

12.0 This section is devoted to the statement, and the proof in part, of a 
lemma which, on the face of it, looks like an exercise in manipulating the 
postulates. In fact, the content of the lemma, and most of the details of 
its proof, are essential in what follows. To postpone them would force us 
into needless duplication of effort. 

Lemma: Let L be a geometrical linear correspondence satisfying P1, 
P2, P4, P51), P61), P71) and the following weak form of P3(1): 


P3’/(1): If pel; 7 then K,(p) = K, = (V, (p)). 


Then there is a frequency domain [, C I,, differing from I, by a 


finite set, such that L satisfies all of the postulates for pel. 

The statement of the dual result is evident and will be omitted 

The proof that L satisfies P3 will be given in this section. Verification 
of the remaining postulates will follow in paragraph 16.6. 

We assume that the properties of positive real (PR) functions are 
known. They are summarized for later use in Section 15. We make 
occasional advance references thereto. 


To the proof: 
12.01 First, K, is a real manifold and for pel, 
K, C (Vin(p)). (1) 
This, with P3’(1), gives P3(1) for L. 

Proof: Ky, is real, as in 7.51. Consider now a pel’, and a veVio(p); 
then |v, OjeL(p). Consider any real jeK, ; then there is a weV,(p) 
such that [u, jJeL(p). Now O and / are real. Hence by P5(1) 

(vy, 7) = (ua, O) (). 
Therefore any real jeK, has a vanishing scalar product with every 
veV,(p). Since K, is real, it is spanned by real j and (1) follows. 


12.1 By the dual of 7.7, if we know that 


(v, kleL(p), 
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then the value of (v, k) is determined by k. This makes it possible to 
state P61) and P7(1) for L (we take P6(1) to include the hypothesis 
that I’, has a finite complement 


12.11 It deK, , then J,(p) is PR 


Proof: if kis real then 
Ji(p v, | v, A 
where, of course, [v, kjel(p). Then however [d, kleL(p), by P4. Hence 
by 12.1, (2) gives us 
Jip) = Jilp 


From this and P6 1), P71) we conclude that JS, (p) is PR for uny real 
eK 
Now, given anv keK, , we have keK, by 12.01. Then 


/ k, + al 


are real and in K, , since Ky, is a linear manifold (see 


yk Jel (p), 


2. Then 


Ip 


for any pel Since each summand in (3) is a PR funetion, it follows 
that J,(p) is PR for any keK 


12.12 Let K, be the set of all vectors keK, such that 


P p 0 for every pel’, . 


Notice that we do not assert that Ky is a linear manifold. 
If fkeK, then eK, and (3) above applies. Then 


, ky) 4 I 0) 


this and the PR property of each summand, we conclude 





FORMAL REALIZABILITY THEORY 
12.13 We wish now to show that K,; C Kyo(p). Consider a real jeK, 
a real keK,. Let 
(u(p), jleL(p), 
[v(p), kleL(p). 
Then, given any real A, by P2 
[u(p) + Av(p), 7 + AkleL(p). 
Then, because /eK, , 
(u + Av, 7 + Ak) (u, J) + Alv, 7) + Au, | 
Since j and k are real, by P5(1) this can be written 
(wu + Av, 7 + Ak) (u, 7) + 2A(v, 7). 


Choose any p; such that Re(p;) > 0. Then P7(1) implies that the 
side of (5) has a non-negative real part at p = p,. The right side, 


suitable choice of A, can have any chosen real part unless 


Re(v(pi), J = (. (6G 


Hence P7(1) implies (6). Now (o(p), 7) is a rational function, by P6(1) 


applied to the other members of (5). By (6), this rational function has a 
vanishing real part throughout the right half p-plane. Hence it is an 
imaginary constant: 

(v(p), J a. (7) 
Then 

(v(p), 7) = (ve(p), 7) = —1a. 

But [v(p), kleL(p), so [v(p), kleL(p) by P4. Since also [v(p), kleL(p), 
by 12.1, we have from (8) that 

(v(p), J) va. 
Comparing this with (7) written for p, we have a = 0 and 


(v(p), J) 0 for pel ,. D) 


Now v(p) was determined by (4) wherein k is real. For any keK, , 
k = ky + tk, , where k; and ky are real and in K, (12.11). A correspond 
ing v(p) satisfying (4) can be written 


v(p) vi(p) + wol(p), LO 
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k,jeL(p), 7 1, 2. Then (9) holds for each of 


P2. where ‘ p 9 
p), v2(p and therefore also for the 7 p) ol (10). 
We have showed now that for any peP, and any keK, , the v(p) of 
1) has a vanishing scalar product with every real jeK,. Since K, 


is spanned by real j, 


prelK, LO - (11 


v(p ~OleL Pp). 
Comparing this with (4), and applying P2, 
v(p v(p a (Q)| QO, hjeL(p). 
Since fis now anv vector in Ky, we have 
K, CK, \p) ae K, 

for everv pel, . 
12.15 We can now also show that Vi(p) © (K,). We return to 12.13 
and read (9) thereof as originally derived for real j and i. Applying 
P5(1), we have from (9) that 

u(p Tk 0 for pel, . (13) 


sy the argument immediately following (9), (13) also holds for any 
eK, , provided / is real. As in 12.11 any jeK, can be written 

, Where j; and js are real, and the corresponding 

u(p) ui(p) + tol p 
where [vu . JrleL(p). Therefore, finally, (13) holds for any u(p) 
satistving (4 Lew, any up eVi(p and any keK, . Therefore 
Vi(p) & (K, (14) 

for any pe 
12.2 We now fix our attention on a specific real poe, 
12.21 


v, h\eL(p 


v, k lel, p 


In particular K p is real 
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12.22 We can now show that K, is a real linear manifold. Consider a 
real keKyo( py). Then (0, kleL(po) and by 12.1 

J, (po) = 0. 


Then by 12.11 (and 15.12), Ji(p) = 0, so keKy . Since Kyo(po) is spanned 
by real & (12.21), we have 


Kk. (p ) c K, Pe 
Comparing this with (12) gives us 
K 10( po) = K, . (15) 


Since Kyo(po) is a real linear manifold by definition and 12.21, we see 


that Ky is. 
12.3 Let us now write, by (12) and (15), 
Kk, K, t K, (16) 


where K. is an arbitrary fixed manifold disjoint from K, and with it 
spanning K,. All three manifolds are real (12.21, (15), 10.6). 
Choose a K; disjoint from K, such that 


K =K; 0K, @K,. 
Let the decomposition of V dual to (17) be (10.6) 

V=V;@0V:0V,. 
Then V; = (K,. @ K, (K,) = Vio by 12.01. Hence 

V=Vwn@OV:@ V;. 
By (14) and the definitions, 
Vio C Vi(p) C Vw @ Ve. 
12.31 Consider a real po. Then by P3’(1), (15) and (16) we have 
K (po) © Ki(po) © Ko ® Kuo(po). (20) 


This is now an expression dual to (19). We shall prove next that, given 


any keK,(p))MK2(= K:), there is a unique rmeVi(po)NV. such that 


fo , AjeL(po). (21) 


Dually, given any veVilp INV: , there is a unique Ir eK, (p, INK, such 
that 


[v, hk JeL (po). 
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- a standard one in algebra and depends only upon P2, 


1 heK p NK there Is some veV, p such that 


P } el. p 


f 


17  OleL Po) 
so, applying P2 to this and (22), 


| 


, kleL(p - (23 


Hence eV, p NV 
reV (poe and 


satisfies (21). Suppose now 


[v3 , kleL(p 
Then using this with (23) and P2 


[1 ve. OleL Pp 


eV Now V,(p AV: is a linear manifold and contains 


veV, NVi(p INV: (). 
Therefore i 
The dual argument completes the proof, 
12.32 The argument actually exhibited in 12.31 uses only P2 and (19), 


hence the e of (21) is unique whether or not po is real. Indeed, this is 

true even when /eK, 

12.33 The result of 12.31 establishes a bi-unique linear mapping between 

K. and Vi(p AV. . Hence these two manifolds are of the same dimen- 
ince K and V 


K* are of the same dimension by construction, 


Vi(p)NV: = V: 


Vi(p Vio @ V2 


snow introduce a real frame in Vand K which provides real 


K, , K; , Ks and in V, , Ve, Vio of (17) and (18). Let / 


l. 
Mi, °°* hm 
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be the basis vectors spanning K,. By 12.32, there are unique vectors 
W(p), ss, U»(p) in Ve such that 
lv (p), /, lel, p?. 
Lete,,--: ov, be the (real) basis vectors in V2 dual to the / 
| 6 | 
Since the u, p) are all in Vz we have for each peP, 
Us\p) S a (pit 
where the coefficients a,,(p) are calculated by (24) to be 
a (p (u.(p), | 
12.41 Because the /, are real, P5(1) implies that 
(s,-(p (u,(p), ks) (us(p), hy) (,.(p). 2/) 


By the reasoning just following (8) and by the uniqueness of the 
u,(p)eVe , since V2 is real, we have u,(p u,(p). Then 


(Qys(p) (u,(p), i) 


12.42 We have by P2 that 
[u-(p) + Au p), ky + AksleL(p), 
for any X. The identity 


(up + u,, k, 
(29 


, ky 


holds in fact for any vectors u, , u,, hy, hy. Using (27), (28) and P6(1), 


it exhibits a;.(p) as a rational function. 


12.5 Consider the m & m matrix [Z;(p)| whose elements are the a,,(p). 
the s-th column of this matrix consists of the components of u,(p). 
The rank of the matrix is by definition the dimension of the space 
spanned by these columns. 


12.51 Now the rank of [Z;(p)| can be expressed in terms of the vanish- 
ing or not of its various minor determinants. There are finitely Many 
such minors and each is a rational function. Each is either identically 
zero or else vanishes at only finitely many points. Hence the rank of 
[Z,(p)], except at these finitely many points, and at the p in the comple- 
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ment of I isa constant. We call this constant the nominal rank of 
Z\(p 
12.52 Let I, consist of all pel, where |Z,(p)| has its nominal rank, 
Phen ry has a finite complement By the realitv result of 12.41, it 
pel : then pel 

It is clear that at any pel, the rank of Zi p does not exceed its 
nominal rank 
12.53 By construction, the vectors (p), ,U,(p) all lie in Vi(pyAV. 
By the reasoning of 12.33, at any real pel, they span V.. Hence the 
nominal rank of {Z,(p)| is m. Therefore, for any pel; ,{Z,(p)| has rank 
mand the w(p), , iu» (p), wing in Ve, still span V. . Therefore for all 
pel’; 


V.(p)NV, = V 


V, p V, 


a fixed manifold, for all pel’, 
12.54 It is clear by its construction (ef. Halmos ’, par. 26) that {Z;(p)| 
deseribes the mapping of 12.32 from K. to V Vip) VW: by 

Z\(p)|lk] 


Here the m-tuples (e] and [A] are the components of 7, and / relative to 


} 


the bases now available in Ve. and K 
12.55 We repeat 

K, © Ki(p) © K, K, © K,. (12 
Fix a pel, anda keKyo(p)NKe. Then [0, kleL(p). Since VeVe , it follows 


from 12.54 that [Z,(p)] annihilates &. Suppose m # O. Since the rank of 


I 


Z\(p)\ is m, it follows that / 0. Hence for pel, 


K.o(p NK (), 


K, \p K, K, 


a fixed manifold. This, with the result of 12.53, proves that L satisfies 
, when m # 0 


O then V O,K Q and (31) follows from (12) and (16). 


p)| is of dimension m and rank m for any pel, . Therefore 
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the correspondence of 12.32 and 12.54 between Ve and Key is bi-unique 
for any pel’, . This extends 12.31 to any pel’, 


12.57 If m = 0, .e., if V. = Kp 0, then V,, K,»)' and the fact 
that ZL satisfies all the postulates is trivial because all scalar products 
(v, kh) for veV, Vi» and keK, K,» are zero. If m # 0, we have vet 
to show that ZL satisfies P5(A), P6(A), PTCA). 

12.6 Since now L satisfies P3, 7.7 


(with 12.1) that if pel, and 


as given Is applicable and we find 


lv, kleL(p), 
then (v, /) is fixed by either v or /&. Furthermore, 
(v, k) (v +, h Ie, 
for any veVio , kyeKyo . 
12.61 If pe, and {v, kleL(p), then veV,, kK, . By (30), 


(16), therefore, there exist voeVio. AkyeKyy such that 
) = (hk — ko)eK, . Then by P2 


fu, jleL(p). 


By 12.6, then, any value assumed by a scalar product (ve, /) with 


lv, kJeL(p) is also assumed by a product (vu, 7), where (32) holds and 


ueVe , jeKe . 


XIII. SUFFICIENCY OF THE POSTULATES 


13.0 We suppose that 1 satisfies the postulates of 12.0. Then the results 
of Section 12 are applicable. The ones of first importance are contained 
in the facts from (15), (30) and (31), that 

V, = Vio us V, 

K, K. t Ew. 


where the choice of Ky was governed only by the requirement that the 
second of these formulae hold. 


ee ae * 
13.01 Considering Ky. and V. as separate spaces, Ve. = Ks by 10.6. 
Let VM be the geometrical linear correspondence between them with 
frequency domain I; and pairs described by 12.31 and 12.56 (or 12.54). 


That is, as vectors in Vo and K, 
lv, AeA (p 
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if and only if, as vectors in V and K, 
lv, kleL(p). 
i the real frame of 12.4 let us renumber the basis vectors so that 

span Vio, 
span V2, 
span V,. 

Then 
span K, , 
span Kyo . 

We sav that such a frame reduces L. 


13.1 Let us now interpret the s-th components of {v] and [A] in this frame 
respectively as the voltage across and the current in an ideal branch 
3 of a2n-poleN, 1 Ss <n 

By construction, the vectors veV, in this frame have components 
a a, QO, since by, *** , Uram Span V,. At the same time, 
the components bey aii, o¢* yO, of [hk] may be chosen arbitrarily without 


altering the fact that , [Ale L](p) because of 12.06. Therefore, the 


ideal branches 3,441, , 6, can each be realized physically by a short 
circuit 

Ina dual way, since / ide ft Ks span K, , any eK, has components 

, », all zero in our chosen frame. Furthermore, the components 
a, of [ve] can be chosen at will. Hence the ideal branches 
3, can each be realized physically by ah open circult. 

Let Ny, now be the 2m-pole whose ideal branches are 8,41, *°° 5 Bram - 
Let the pairs [fe], [/]] admitted by N, at each pel, be the [fv], [Al], 
where [r, AJeM(p) (13.01). The representation just found for N> shows 
that Nis physically realizable if and only if N; is. 


[3.11 The matrix {Z;(p)] of 12.54 is the impedance matrix of the 2m- 


pole N 

13.12 We now show that [Z(p)] is a positive real matrix. The displayed 
formulae of 12.41 show (ii) and (iii) of 1.1, and 12.42 shows (1). Now 
suppose that {r, AJeM(p). Then, as vectors in V and K, [v, kleL(p) by 
definition of AM p). Then, however, if & is fixed 


} 


Ji P) (v, k 
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isa PR function (12.11). Regarding v and k in V2 and K, let 
[bras , “*-e = k}. 


Then by (1) of 7.0 


(v, k) = > > Ast( P) Deg eDere 
s=1 


t=] 
and this has a non-negative real part if Re(p) > 0. This is (iv) of 1.1: 


13.2 We can now prove the lemmas 8.1 and 8.2. Given a linear cor 
respondence [L] which satisfies Pl, --- , P7 by 11.13 we can interpret 
{(L.] as the concrete correspondence representing a geometrical cor- 
respondence L in some chosen real frame, and L satisfies Pl, --- , P7. 
Then by the results in 13.01-13.12 there exists a real frame in which the 
representative [L]; of L has the properties claimed in 8.1 and 8.2 for Ly. 
But we saw in 11.22 that [L] and [LZ], are related by a real matrix W like 
the L and Ly of Section 8. Q.E.D. 

13.21 With the proofs of 8.1 and 8.2 we have reduced the sufficiency 
claimed for Pl, --- , P7 in 8.0 to the sufficiency of positive reality of 
{Z(p)} claimed in 1.1, by the argument outlined in 8.5. 


XIV. OPERATOR-VALUED FUNCTIONS OF DPD 


The next three sections are directed principally toward the proof of the 
matrix theorem of 1.1. They do however, contribute to 12.10 and to 
the necessity proof. 

14.0 We continue to use the geometric language. The reader who re- 
gards this as unduly pedantic is free to place a concrete interpretation 
upon every argument, for all of the arguments are either frankly based 
on matrix representations or upon the three identities: 

14.01 (Zj, k) = (Z*k, 7) for all 7, keK. 

1402 Zk = (Zk) for all keK. 

1403 Z’ = (Z)* = (Z*) 

Geometrically, the first and second define Z* and Z, and the third 


defines Z’ in two ways. The equivalence of these two ways is a theorem 
based on (10) of 10.33. 


14.04 These identities are obvious for matrices using 7.0 and 7.2 


14.05 The symbol Z will always denote an impedance (operator, matrix, 


scalar), and Y will always denote an admittance. An impedance oper- 
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ates from K to V, an admittance dually. The operators in Halmos’ 
are physically dimensionless, in that they operate, e.g., from V to V. 
This difference is scarcely noticeable 

We shall regularly omit the duals to concepts or proofs given in terms 


ol impedances In doing so, we adopt the rule that the dual to an 


CXPresslon 


Zk, | 


7 
14.1 An operator is called symmetric if Z Z’. Such operators have 
three useful special properties: 
4.11 Uf Zissvmmetrie and 7 and / are real, then 
Zj, | Zj, | Z)*k, j) = (Z'k, j) = (Zk, j) 


by (10) of 10.33, 14.02, 14.01, 14.03, and hypothesis. 


L4.12 Let t thy , Where hy and dy are real (10.42). If Z is symmetric 


(Zhe, ky) 

(Z (tks), hy). 
Cf. the similar identity in 12.11 
14.13 The symmetric operator Z is completely defined by the quadratic 
form 

(Zk, fb (1) 

isa funetion of real keK. For 14.11 permits the formula (29) of 12.42 
In any real frame, where wu Zk. Vhe matrix elements of {[Z(p)] in 
that frame are then defined by that formula in terms of values of (1) 


for real 


The form (1) specifies any Z (symmetric or not) if fis allowed to 
range over all of K (Halmos’, par. 53). 
14.2 Let Z Pp) TOW be an impedance operator depending Oh p. We say 


that Po ~ *% is i pole ot order m ot Z p il 


(UJ lim. (p po) (Z(p)k, ke) 
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exists for every keK and is not identically zero. By 15.13, this limit 
((k) defines an operator Ro , the residue* of Z(p) at po , by 
(Rok, k) (( ke) for keK. 

The changes in (2) required to define a pole at p © are obvious 
14.21 A pole py of order m of Z(p) is a pole of some matrix element of 
[Z(p)], of order m, in any frame, and no element of {[Z(p)] has a pole at 
po of order exceeding m. For the elements of {|Z(p)} are defined by the 
values of (Z(p)k, hk), by 14.11 and Halmos’ loc. eit. 


XV. POSITIVE REAL FUNCTIONS 
15.0 Let fi ~P) be “a scalar function ol the complex variable p. Following 
Brune’ we define f(p) to be positive real if 

(1) f(p) is a rational function of p, 

(i1) f(p) tp), 

(iii) Re(p) > O implies Re(f(p)) > 0. 

The property (1) of being rational is of course on a quite different 
level of ideas from the other properties, but it saves words later to in 


clude it specifically in the meaning of positive real. 
We abbreviate the words positive real to PR. 


15.01 The open region of the complex plane consisting of all finite p 
such that Re(p) > O—the right half plane—we denote by I, 


15.1 Brune, loe. cit., established a number of properties of PR functions 


f(p) which will be useful to us here: 
15.11 f(p) has no poles in TP, 


15.12 If Re(f{(p)) = 0 for some peP, , then f(p) = O for all p. 


- ~~ l 
15.13 If it exists, . is PR. 
I(p) 


15.14 If f(p) has a pole at p = po, it has one at p = p 
15.15 Hf f(p) has a pole at p = tw, that pole is simple an | 
2p 
(p) = 4 gat + hi), 
where r > 0, and fi(p) is PR. 


* Properly, Fy is a residue only when m 1. There is no convenient name 
available for general » 
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15.16 Tf fp) has a pole at p ©, that pole is simple and 


f p . 


} 


where > YU, and f,(p 
15.17 We shall use all of these in the next section, save 15.13. Our aim 
is to prove properties analogous to 15.11, , 15.16 for PR matrices 
iid operators 

The reader familiar with the Brune process” for realization of a 2-pole 
vill remember the Importance of the properties 15.11, --- , 15.16 for the 
success of that process. Correspondingly, we must establish the analogs 


of these properties to Implement the general Brune process for 2n-poles. 
pro} | | 


XVI. POSITIVE REAL OPERATORS 


16.0 An operator Z p from K to V will be called positive real (PR) if in 


ome real coordinate frame the matrix {[Z(p)] is a PR matrix in the sense 
of 1.1--that is 

Z(p | has rational elements Z p 

ZirdkpP Z,s(p 

Z <p Z.(p 

lor anv real /eK and any pel’, 

Re(Z(pyk, | pe. 
We intend in this seetion to establish for PR operators the properties 

listed below. By subtracting 0.9 from the designation of each property 
one obtains the designation of the analogous property of a PR sealar 


function, stated earlier 
1.0L Z p has no poles in r'. 
16.02 Tf Re(Z(pyk, | 0 for some pel, , then Z(p)k 0 for all p. 


16.03 If it exists, Z"(p) = ¥(p) is PR. 


16.04 Tf Z(p) has a pole at p = po, it has one at p = po. 
16.05 If Z(p) has a pole at p tw) , that pole is simple* and 
“P 7 tn 
Z p - R T Z\(p), 
p T WwW 
Where Fis real, symmetric, and semi-definite, not zero, and Z,(p) is PR. 


one 
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16.06 If Z(p) has a pole at p 2, that pole is simple and 


Z(p) pR + Z,(p) 


where h I’ R,R > Oand Z(p) is PR. 


16.07 There is property of rational scalar functions f(p), whether PR 
or not, that is essential in the Brune theory: the existence of a finite 
integer, the degree of f. Each step in the Brune reduction of f(p) leaves 
an unreduced portion which is of lower degree than the funetion upon 
which the step was performed. The finiteness of the original degree of f 
then guarantees the termination of the process in finitely many steps. 

There exists also for rational matrices (and operators) a concept of 
degree. This degree plays the same role in the general Brune process for 
2n -poles as the degree of a scalar function does in the process for 2 poles. 
To define this degree and develop its properties requires an excursion 
into classical algebra. Since we shall not need these ideas until Part II 
we defer further discussion of them to that part 


16.1 If Z(p) is PR it follows at once that the matrix [Z(p)] is PR in any 
real frame. 


Proof: Two such matrices are related by 
A(pyh = (UZ @ILY 


where ( is real, by 11.22 and the argument in 8.6. The PR properties 
of [Z(p)| are oby iously preserved by this operation. 


16.11 If Z(p) is PR, then 
Z(p) Z'(p) = Z*(p) Z(p). 
Proof: Use 16.0 and 14.03 in a real frame. 
16.12 If Z(p) is PR, then for any given /eK the function 
Jip) (Z(p)k, k) 


is a PR sealar function. It follows that the limitation in 
to real & is a simplification, not a restriction. 

Proof: J,(p) is independent of coordinate representation. 
real frame, (i) of 16.0 implies (i) of 15.0. 

By 14.01 and 16.11 


Ji(p) = (Z*(p)k, k) = (Z(p)k, k) = Jel). 


This is (ii) of 15.0. For any k, 14.12 and (iv) of 16.0 imply 
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ersely to 16.12, 1f 7 p) is sVimetric and Jy (p Is PR tor 
then Z(p) is PR, and J,(p) is PR tor all k. 
p) is rational so 4) of 16.0 holds in any frame by 14.13. 
of 16.0 holds 

bv (10) of 10.33 and 14.02 


Jdp) = dip Z(p)k, | 


Hence Z(p Z(p) by 14.18. This is Gi) of 16.0, and (it) there holds 
Py hypothesis 
16.2 Proof of 1.01: By 15.11 and 16.12, /;(p) has no poles in TP, 
his is 16.01 by the definition 14.3 of pole 
1.21 Corollary: Any PR Z(p) can be considered as defined throughout 
I : for anv &k, 7 pP) Is defined throughout i hy 6.2. lor each P, 
asa funetion of hk, J,(p) defines Z(p) (14.13). 
16.3 Proof of 16.08: In any frame |Z (p)| Z(p)|’ [}'(p)} consists 
of rational elements, by direct caleulation of the inverse matrix. Ina 
real frame | ¥Cp Z (p)\issymmetric and real for real p by the same 
argument (both facts are also deducible geometrically). Hence we have 
the duals of (2), Gi) and Gil) of 16.0 for ¥(p). Clearly Y(p) is defined 
throughout P 
Now suppose that for some veV and some poe, we have 

Reiv, Vipe) < 0. 

Phen there is a veK such that : Z(podk. Therefore 
Re Zipok, | Re(Z(podk, hk) « Q). 


s impossible, we have the dual of (iv) of 16.0 for ¥(p) and 


16.04: This is immediate from 15.14, 14.3, and 16.12 


of 16.05 and 16.06: Suppose rf p has a pole at p lw 


Phen (Z(p does and that pole is simple by 15.15 and 16.12. Phen 


coh write 


p 
s regular at p lw). Now Z)(p) has a pole at p 


ya similar argument gives 
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where Z;(p) has no pole at ta or —tw . The symmetry of Z and linear 
independence of the terms above then imply the svmmetry of Ry, Rk, 
and Z\(p). 

For any keK, now, 

(Z(p)k, kh) = ! (Rok, k) 


p lw 
Applying 16.12 and 15.15, 
(Rok. h (P,/ x i, a 0 


for all k. Hence RP R, R (say) and FR is semi-definite. Also, 
(Z\(p)k, hk) appears as the residue f,\(p) in 15.15 and is therefore PR. 
Then Z\(p) is PR Vv 16.13. With Ry and AR, identified, (1) above is the 
expansion given in 16.05. We have now proved all of 16.05 save the 
reality of FR. But 

») 

2p 

»R 

pte 
is PR, by 16.13, hence is real for real p. Therefore Ff is real. 
The proof of 16.06 is similar. 


16.6 To prove 16.02 we appear to digress somewhat, by first com- 
pleting the proof of the fundamental lemma of 12.0. It was established 
in Section 13 that the matrix [Z)(p)] deseribing W(p) in the chosen 


basis is PR. The case in which it is nonsingular (i.e., m # O, ef. 12.56, 
12.57) remains to be examined. 


16.61 Tf [Z,(p)] is nonsingular then its inverse is PR (16.3). Then for 
any veVe , 
(v. A = (v, ) pr (2) 


is PR (16.12 dual). By 12.61, for any weV,, the values of the function 
F’,(p) are the values of (2) for some veV.. Hence F,(p) is PR. This is 
P6CA) and P7(A) tor L. 

16.62 To settle P5 for L in 12.0, consider pel, and 


vr, kleL(p), — u, jel p), 


where uw and v are real. Then, say, 


where meV po) , (ye Ve. But then 
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and, because Vy,» and V2 are real, i ' , , and these vectors are 


real, Using similar reasoning for 7, 


(ma, , (pyr), (3) 


by 12.61. The equality (1, / (v, 7) now follows from (3) and the duals 
of 16.11, 14.11. Hence we have P5CA) for L and 12.0 is proved. 
16.7 We now prove an important 

Lemma: Let Z(p) be a PR operator from K to V. Let Ty, be the set 


ol p where Z p) is defined and has a rank equal to its nominal rank. Let 


L be the correspondence with domain TP, and pairs 
[Z(p)k, k], keK, . 


Then 1 satisfies Pl, --- , P7. 
Proof: L satisties Pl and P2 (6.3). T, satisfies P4 by the argument of 
12.52. Then L satisfies P4, for by 16.11 


Zip)! Z(p k. 


L satisfies P51) by 14.11 and 16.11. Ty, satisfies P6 by 12.52. Then L 
satisties POUL) and P7(1) by 16.12. The fundamental lemma, 12.0, 
now proves that L satisfies all the postulates. 
16.71 We call a correspondence satisfying all the postulates PR. 
16.72 Proof of 16.02: Suppose Re(Z(pok, hk) Q for some poel, . 
Jecause this function of pis PR (16.12) we have 

Ji(p) (Z(p)k, hk) (). 
Hence eK, Kio (12.12, 12.55). Hence [0, kleL(p) for every pePy. 
Phat is 

Z(pyhk Q for peD,. 
16.73 Corollary: If Z(po)l QO for some poel'y , then Z(p)h 0. 
lor the hypothesis here implies that of 16.72. This is the analog of 
15.12; the result of 16.02 is stronger. 
16.8 An important consequence of 16.7 is the 

Lemma*: Hf Z(p) is PR and of rank m, then there exists a real coordi- 


nate frame in which the matrix [Z(p)| is an m X m nonsingular PR 
matrix {|Z p | bordered by ZCTOS 


* Proved by Cauer 
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Proof: Consider the PR correspondence L defined by Z(p). Then 
Vio = 0, because Z(p)0 = 0 for every pel’, . Consider the real frame 
of 13.02. [Z(p)] in this frame takes any of ky, _41, °** , ky into 0 because 
these span Kyo . Within Ke , [Z(p)] must describe the same operation as 
the [Z,(p)] of 12.54. Because [Z(p)| is symmetric the lemma follows. 


XVII. THE JUXTAPOSITION OF CORRESPONDENCES 


17.0 This section and the next will consider ways of constructing new 
correspondences from old. This will provide the basis of the necessity 
proof of Section 19. 

17.01 It is obvious that if two physical networks are set side by side 
and their accessible terminals regarded as the terminals of a single 
larger network, that enlarged network is again a physical network. 
This is the gist of the present section. 


17.1 Suppose that 


V Vi @ Vz, K K, @ K, 


where K; = Vv: and all spaces are real (10.6). Let Ey project on V 
along Vz (Halmos’, par. 33) and F, 1 — FE, project on V, along V; . 


rr * . ) 
Then £; projects on K; along K,, 7 # ¢@ (Halmos’, loc. cit.). It is 
»* 


easily verified that F, E;, Ej 
dimensionless operators. 


* . ° . 
kh; , from the analog of 14.02 for 


Considering V,; and K; as separate spaces, let L,; be a geometrical 
linear correspondence between them with frequency domain I, , 7 1, 2. 
Consider the correspondence L between V and K defined by 
(i) The frequency domain [, = nA. 
(ii) [v, kJeL(p) if and only if [Ew, Ez kleL,(p), = j, 2. 
In (11), of course, we regard Ew and Erk as elements of V; , K;. 
17.11 L so defined is called the juxtaposition of LZ; and Le . 
17.2 Lemma: L is PR if and only if each of 1, and Le is PR. 
17.21 Proof of “if: It is clear that L satisfies Pl and P2. Further 
notation is now simplified if we put Ll; = M, L. = N. Consider the 
manifolds 


Vu © Vw, V wo r Vwo, K wv t Ky, K wo ¢ Kee. 


where Vy © V; is the manifold of voltages admitted by L, M con- 
sidered as a correspondence between V,; and K, , and V wo the manifold 
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of voltages reV, such that |v, OleL,(p) for all pel). Dual definitions 
for Ky, Kyo, and symmetrical ones for Vy, --- , Kyo need not be 
repeated 

It is clear from these definitions that the four manifolds above are, 
in the order listed, the manifolds 


V_, Vo, K., K, 
for L. Now, for example, 
K, K w @ Ky K wo)" (Kyo)” 
by 10.6. This last manifold, in V, is (Vy @® Ve) A (Vy @ Vi), byP3 


for Mand NV, and by 10.6. But by direct calculation 
t V NV, t V, vw t Vy Vi . 


The dual of this result then completes P38 for L. 


: * 
P4 for LZ is immediate because the /, and /; are real. 


Phe duality of the decompositions of Vand K implies the identity 


v, | tw, kyk) + (Rw, Eek) 
that is hk DEY De QO, and dually. This is Halmos*, par. 33). All 


of Pd, P6, and P7 for L follow at once from this identity. 


17.22 The “only if? of 17.2 is a special case of the result of Section LS, 
Its proot will be deferred to 18.4 


17.23 It is obvious that the notion of juxtaposition and the lemma of 


17.2 extend to Juxtapositions of more than two correspondences, 


17.8 Keven without the “only if” part of 17.2, we have enough for the 
following characterization of PR correspondences : 

Theorem: \ correspondence Lis PR if and only if it is the Juxtaposi- 
tion ol 

i) a correspondence defined by a nonsingular PR matrix between 
a ViandaK, = Vy, 

i) a correspondence consisting of short circuits: that is of pairs 
0, | for all keKe and all p, 

ll) a correspondence consisting of open circuits: that is, of pairs 
v, O} for all veVy and all p 

Proof Ht Lis PR, the decomposition indicated is that of 13.1, 13.11, 

12. It Lis the juxtaposition indicated, then it is PR by 16.6 and the 
‘Win 17.1, provided the short and open circuits are PR correspondences. 
Phe verification of the postulates for these latter is easy and will be 


omitted 
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17.31 The labor of considering PR correspondences instead of matrices 
has vielded the disappointingly simple result of 17.3. We have already 
been warned of this, however, by our knowledge of the properties of 
physical networks (2.9), 


XVIII. THE OPERATION OF RESTRICTION 


18.0 In addition to juxtaposition, which is an operation on correspond- 
ences clearly motivated by physical considerations, there is an operation, 
here called restriction, which has important use in the next section. 
There the physical meaning of the operation will become clear. 
IS.1 Let Vand K V* be a pair of dual spaces. Let U and J U* be 
another pair. Suppose that ( is a given fixed linear operation from J to 
K: given any jeJ, there is a unique k())eK, written 
I()) G7, 

such that if /;, Jes? , 2, then 

ayky + ash Clan 4 
for anv complex scalars a; , de. 


8.11 Let (v, k), denote the scalar product between V and K, and 
(u, j)y that between U and J. Given C, and any veV, let us find that 
unique vector u(v)eU for which 


(uly), 7)e = C9) (1) 


for every je. That such a vector u(v) exists and is unique follows from 
10.13 when we notice that the right-hand side of (1) defines a function 
conjugate linear inj. Now for fixed 7, the right-hand side of (1) is linear 
inv, hence so also is the left side. That is, there is a linear operation C* 


from V to U such that 
wl 


The following chart illustrates the situation: 


, = 
cf! Cc 
U J 


18.12 We suppose now that C takes real j into real /, i.e., that ¢ 
real. Then by (1 


(C*v, 2 = (C*v, 7) 
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By comparison with (1), we have 
Cs C*i 
Hence C* also takes real vectors into real vectors and is real. 


18.2 Now let L be a PR correspondence between V and K. We define 
one, say M, between U and J, as follows: For each pel, , let M(p) 


consist of all pairs 


C*p and 
lv, CyleL(p). 


This definition ean be illustrated by enlarging the chart of 18.11: 


The u’s corresponding to jeJ can be constructed by going around 
through C, L and C*. This then defines a direct mapping from J to U. 


18.21 We call the Vv defined by 18.2 a restriction of L, since its pairs 


are images under C* and C’ (which is not defined over all of K) of a 
restricted set of pairs drawn from L. 
18.22 Clearly there is a dual operation defined by an operator D from 
U to V. We might distinguish the operation of 18.2 by calling it a 
current restriction, its dual by calling it a voltage restriction. 
18.23 The restriction M of ZL is defined by lists M(p) which exist for 
any pel’,. The frequency domain of M has not yet been specified, 
however 
18.3 Theorem: If L is PR, then there is a frequency domain Ty for M 
such that 17 is PR 

Proof: Pl and P2 for M are evident at once, for any pel,. The 
remainder of the proof is divided among 18.31, -+- , 18.37 below. 
IS.31 For P8, let Jue be all yeJ such that CyeK,. Then, given jeJu, 
for each pel’, there is a v such that 

lv, CjleL(p), 


whe hee 


[C'*p, gJeM (p). 
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Therefore J y(p), the space of currents admitted by M at frequency p, 
coincides with the fixed J y at each peI’, . 

Clearly J w is a real linear manifold. 
18.32 Consider now U yo(p): if [u, OJeM(p), then there is a v such that 
u = C*v and 


[v, CO} [v, OleL(p). 


t 


Hence veVio(p) Vio for each pel’, . Therefore, for each pel; , 


U wo(p) CC Va . (2) 


Now suppose, conversely, that pel, and veVy Vio(p). Then 
[v, OleL(p). Now 0 = CO, so [v, COleL(p). Hence [C*r, O]eM(p), so 
C*veU yo(p). This proves the inequality opposite to that of (2), so for pel’, 


U wo(p) C*V, Uw, 3) 
a fixed space. 


18.33 Now consider (U ue) If Jel U wo) , then 


for every “eU yo. That is, by (3), 
(C*v, J)2 = (v, Cy) = 0 


for every veVio. Therefore Cje(Vio)” K,, and jeJ yw by 18.31. That 
is, we have proved 


and, combining 18.31 with this and (3), 
J u(p) = Ju 2 (Um(p)) (U wo). (4) 


This is the weak form P3’(1) of 12.0 for A/Z. It is as far as we can go 
with P3 at the moment. 


18.34 Consider P4. If for pel’, we have 
[w, le (p) 


then [v, CyJeL(p) and u C*v. But then [é, CyJeL(p) and a 
18.12. Then however 


[w, jleM (p) 


by definition of 7. This is P4. 
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18.35 Consider P5(1 


eL(p). Since C 


by Pd5(l) for L. This with (5) for r # s proves P5(1) for MV. 
IS.36 Fix a yeJ y and for each pel, a u(p) such that 


u(p , le M p 


vp), CyleL(p), 


lor some p Phen as in (5) above 


(p), J nip), ©) 


Pol) and P7(1) for 1 then imply that P6(L) and P7(1) hold for JZ, 


using I’, for Py in P6. 


18.37 We now have 


Vsatisfving the hypotheses of 12.0. Therefore 
there is: 


i Ty such that J satisfies all the postulates. This is 18.3. 


18.4 Proof of “only it” wm 17.2: 


.2: Suppose that L between V and K is 
the juxtaposition of L,; between V,; and K,, Le between V2 and K.. 
U V, and J K, . Let C be the identity map from K, to K: 
K, , then () is just 7 considered as a veetor in K. Then ( is 
is easily computed that C* is ky . 
Consider the restriction Woof L based on this C. [ts pairs for 
pel y Cl, are all the pairs [w, 7] such that 7 K*jeK, and u 


where 


v, JleL(p). 
But then 


i, JI Kv, E*)| 
and this isin Ly(p) by (6) and the definition of juxtaposition. Therefore 
the list V/ p) is contained in Ly, p 


Suppose that 


uw. lel. Pp). We have |0, Olel. p) so by P2 and the defi- 
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nition of juxtaposition 


MH, J\eL. p 


But then / BB}, t Ku, and by definition of 


lu, jie Mp 


Therefore for every pel y, Mp Li(p). Therefore there 
queney domain (Ty) for 1, sueh that 1; is PR. 


XIX THE NECESSITY PROOF 


19.0 Fortunately for this section, those parts of network theory which 
we require have recently been very succinctly stated by J. L. Svnge™ 
We shall paraphrase them here, referring the reader to the source” for 
details of definition. 


19.01 First, we observe that in Cauer’s definition’, which we = shall 
repeat dn detail below, an ideal transformer with m windings is a 2m-pole 
whose terminal pairs are the termini of the respective windings 

A system of m coupled coils is a 2m-pole with similarly defined terminal 
pairs. 
19.02 Given a 2n-pole N which is a finite passive network, let us adjoin 
ideal transformers as in Figure I(b). We then draw the ideal graph of 
this network. Adjoin to the graph ideal generator branches y, , 
¥.. 7, between 7, and 7. , 1 <r<n. Let 8, be the ideal branch repre 
senting the transformer winding between 7, and 7), 1 <r <n. Enu 


merate the remaining branches of the graph 8,41, °°° , 8». 

19.03 The branch y, is in a mesh with 8, and no other branches. Let us 
call this the r-th external mesh. Anv basic set of meshes must include 
each of these. 

19.04 Let (4, ---, ¢, be the currents in the generator branches, 
hy, s+, hy the currents in the branches 3; , --+ , 8) and 


» Ry, ot, Wi 


w, be the voltazes across the generator branches, 


, » the currents in the 3) , - , By and 
lw = lw Le. ? Sere | 


ls» 


19.05 Let us choose a basic set of meshes, let j;, --° , J. be the respec 
tive mesh currents, and 
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Me} 
’ } 


be thre tuple of mesh voltages. We SUppose that Be * os 


reve! respectively to the nm external meshes. (Cf. 19.08.) 


19.06 The results of Synge can now be stated as follows: 


There exists a real constant matrix [Cy] of s columns and b + n rows 


in fact, elements which are +1, 1, or O) such that for any [7] 


[C4] | (1) 


t 
asetof branch currents satisfving Kirchoff’s node law, and for any [w] 


is a set of mesh voltages satisfving NWirchoff’s mesh law. Furthermore, 
viven any [f| which satisfies the node law, there is a [j] such that (1) 


holds 


t 


19.07 If we interpret the [¢], [jy], ete., as representations in real bases 


1 
] 


* 
I*. 


then [Cy] is real and [C,]’ IC’ 


19.08 The matrix |()] has the form 


where [Cy] is an nm & n diagonal matrix (having diagonal elements +1, 
in tact 

Proof: By construction, ji, -** , J, are mesh currents in the external 
meshes, These are then equal, suve for sign, to the currents f,, °°° , fy 


in the generator branches. 
19.00 By 19.08, (1), and the detinitions in 19.04, 
[C}lyl, 


; 


and by 19.07, {¢ 


10.1 Let us suppose that we have enumerated the branches 8,4: , 
in 19.02 in sucha way that 8,,.,-°°,8-areall the two poles in the 
graph ) m , ba are all the branches containing coils which are 
magnetically coupled, and By.) , -°* , 8 the remaining ideal branches of 
ideal transformers 
Let |Zu(p)| be the (d n) XK (d n) impedance matrix relating the 
voltages across the branzhes 8 , By to the currents in them when 





dies 
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we consider the individual two-poles and the system of coupled coils 
as separate unconnected networks. Then [Za(p)| is composed of a 
(e — n) X (ec — n) diagonal matrix in the upper left field and a (d — ¢) X 
(d — c) matrix in the lower right, with zeros elsewhere. 

19.11 The diagonal part of [Z(p)| has elements drawn from the follow- 
ing list: 

(i) f(p) p 

(11) f(p) = dp 

(ill) f(p) = Ap 
where Pp, 6, \ are non negative constants, possibly different for each 
branch. 


19.12 It is shown in texts on electromagnetic theory that the matrix 
representing a system of coupled coils is of the form 


PIG), 
where [G] is a real, constant, symmetric, and semi-definite matrix. 
The lower right field of [Z¢(p)| is then such a matrix. 


19.13 It is obvious from this description that [Za{p)| is PR. It therefore 


describes a PR correspondence between (d — n)-tuples of current and 
voltage. 


19.2 We must at last consider ideal transformers in detail. Let V; and 


K, be m-dimensional spaces represented as aggregates of m-tuples. 


Let p:, po, *** , Pm be m real numbers. Let Vy consist of all m-tuples 
la} la, , °-+ , dmleVy such that 


ay a» 
pi pe 

We interpret these relations as follows: 
(a) If any p, 0, then a, = O 


(b) If any two p, , ps are not Zero, then 


(c) If only one p, # 0, then a, is arbitrary. 
Let Ky consist of all m-tuples [b] = [hb , «++ , bmleKy such that 


ta 


{ l 


pili + pobo + -°* + pmbm 0. 


V, and K, are linear manifolds. 
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Let {/ be the concrete linear correspondence defined by the list 
iL p which COUSIStS lor each complex p ol all palrs [la], [b]] where 
aleV>. (b|eK, 

The correspondence deseribed by [7] is what Cauer” defines as an 
ideal transformer. He shows, loc. cit., how it can be defined as the 
limiting case of a physical transtormet! 

‘| hye re is also il dial kind ol dey ice. deseribed by FY Y correspondence 


admitting all [bJ/eK, for which 


and all {aleV, for which 
\ (ly 


This also is an ideal transformer obtainable as a limiting case of a 


phy sical one 


1921 The correspondence Lyris PR 
Proof: We observe that V, €¢ . for let laleV, ; lhleKy , and let 


(be the common value of the a p Then 
a,b (Xp h {(Xp,b,) (). 


Phe postulates are now all easily proved. We omit the details. 


10.3 Let Vand K be 6-dimensional spaces. We interpret the 4-tuples 
and {hk} of 19.04 as representing vectors veV, keK in a real frame. 
Let L be the correspondence between V and K formed by Juxtaposing 
1) the correspondence described by |Za(p)] relating components with 

indices in the range n + | tod, 

the several correspondences described by ideal transformers, 

relat y components with indices in the ranges lL ton and d + | to h. 

Lis PR because it is the juxtaposition of PR correspondences. 


31 Let U and J be s n-dimensional spaces. We interpret the [1] 


of 19.04 as representing weU, jeJ in a real frame. 


Let © be the operation from J to K whose matrix in our chosen 
frames is [C|. Then C* operates from V to U with the matrix 
cr ( By these definitions, ( is real. Let .W be the correspond- 
ence between U and J obtained by restricting 1 with C. Then there ts a 
frequency domain Poy such that J is PR (18.3). 


194 By 19.09, [AZ] in our chosen frame is the correspondence estab- 


lished between mesh currents and mesh voltages by. the network of the 
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2n-pole N. When this network operates as a 2n-pole, the only mesh 
voltages which are not zero are those relating to the external meshes, 
since there are no internal sources of voltage. We must now account 
for this. 


19.41 Let Ve, Ky be n-dimensional spaces. Choose a real frame and let 
D be the operation which takes 


la,, °°: , a,jeV 
into 
lay, «++ ,a,,0, °°: , OfeU (4 
in the frame of 19.31. Then D is real and D* in the chosen frames takes 
bh, °-+, ble) 5) 
into 
[br, --+ , OnleKe. (6) 


19.42 We interpret the n-tuples (3) and (6) as voltages and currents 
in the external meshes of N. Their relations to (4) and (5) are con- 
sistent with this interpretation. 

Let us restrict 1/ by D, to get a correspondence M, between V. and 
K,. In our chosen frame, the passage to [.1/;| corresponds, by (3) and 


(4) of 19.41, to considering mesh voltages in N which vanish for every 


internal mesh, and, correspondingly letting the mesh currents adjust 
themselves to this situation. We of course observe only the external 
mesh currents (6). 

19.43 Mo owas PR. So, therefore is MM, (18.3 dual). Since [.1/,] is the 
correspondence established by the physically realizable 2n-pole N, 


the necessity of Pl, --- , P7 for formal realizability is established 


XX. APPENDIX TO PART I 
20.0 We must prove 7.22 and those assertions of 10.6 which are not 
covered in Halmos’. These concern reality. 
20.1 Let V, be a real manifold and 
V=V,90V:, K K, @ K, 


where K, = (V2), ete. The basis (14) of 10.6 exists by Halmos’, par. 19. 
We show that it can be chosen to be real. We have linearly independent 


vectors 
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where the first r span V,, the last n 
’ u, + 1, , 
,w, are real (10.42). Since V,; is real and a linear manifold, 
+ », eV; , ] 
r. Among the 2n real vectors 
,U,,? pg ©* 5 hy gpa yo Se (1) 


Zr are in V,, and they span V, because the v,, 1 << s <r, can 

- constructed from them. The whole list (1) spans V, because from it 
all the vw, 1 << s <n, can be constructed. Since the v,eV2 do not use 
in their construction any of the first 2r vectors (1), it follows that the 


last 2(n r) vectors in that list must contain a set spanning V2. The 


reality of the veetors (1) then establishes the existence of a real basis, 
Say, 


; / / (sy 
Se ee (2) 


which provides a basis in V,; and V2. 
20.11 We now have 7.22. The unique dual basis 


ee 


to (2) is real by 10.41. Hence all of V;, Ve, Ki , Ke are real. The proof 


of 10.6 is then complete. 
20.2 If in a real basis (2) (dropping primes) 


Qyv, + 


then by 5 of 10.3 





[ 
Hy . 


Phe geometrical conjugation of 10.3 is therefore simply the concrete 


one of 7.2 in any real basis. This proves the remark of 10.35. 
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An \ pplication of Boolean Algebra to 
Switching Cireuit Design 


BY ROBERT E. STAEHLER 
Manuscript received January LO, 1952 


This paper discusses the application of switching (Boolean) algebra to the 
development of an all-relay dial pulse counting and translating circuit em- 
ploying the minimum number of relays. An attempt ts made to outline what 
appears lo be the most promising method of obtaining beneficial results from 
the use of the algebra in the design of practical switching circuits. 


INTRODUCTION 

The demands made upon telephone switching svstems in regard to im- 
provement Sh handling Capacily, speed, flexibility and economy are Con- 
tinually increasing. In order to meet design objectives enabling the 
fulfillment of these demands, switching circuits have of necessity become 
more and more complex and intricate. As certain types of relay switeh- 


ing circuits increase in complexity, the problem of control and output 
contact network design becomes more and more laborious and time con- 
suming. This is especially true in those circuits in which an attempt has 
heen made to achieve the ultimate in efficiency and economy in that the 
number of relays used therein approaches the absolute minimum neces- 
sary to provide the required number of distinet output combinations. 
In this type of near-minimum combinational or sequential relay cireuit 
there are numerous parallel control and output contact paths which 
thread through the same relays repeatedly, thereby causing the indi- 
vidual relay contact loads to become relatively large. Thus the designer’s 
problem becomes that of first developing a workable control and output 
contact network and then manipulating and minimizing contacts within 
that network so that the maximum number of contacts used on any one 
relay is within that permissible on any commercially available relay 
having th Hhecessary speed characteristics 

Even in those combinational and sequential relay circuits which are 
not near-minimum and therefore probably have fairly light individual 
relay contaet loads, there are, of course, advantages to be gained by 
using the least number of contacts possible. Although the initial cost per 
idditional contaet (assuming that a few added contacts per relay will 
hot impair the relay speed or space characteristics to an extent that the 


circuit requirements are not met) is almost negligible, there are other 


20) 
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economic savings possible, Since each contact must be connected to the 
remainder of the contact network, minimizing contacts and consequently 
soldered connections means a saving in wiring time and labor. Further- 
more, if the designer will manipulate the contacts so that the relays can 
be chosen from a comparatively few standardized codes, which are in 
large demand, it is possible to avoid the expensive stockpiling of numer- 
ous special designs having only a limited demand. In addition, using 
the least number of contacts minimizes the focal points of most relay cir- 
cuit failures which are the contacts themselves (.e., dirty or Worh Ccon- 
tacts). 

It might also be noted at this point that electronic combinational or 
sequential circuits usually require electronic gating networks to perform 
functions which are completely analogous to those of relay contact net- 
works. Hence, the same problem of minimization exists. However, in 


electronic circuits, gate minimization is even more advantageous since 


the cost per additional electronic gate is much higher ‘han the cost per 


additional relay contact. 

It is rather obvious that the multiplicity of paths in most combina- 
tional and sequential circuits can cause their design to become an ex- 
tremely difficult and time consuming problem if the contact paths are 
developed with the aforementioned considerations in mind. 

The cireuit designer’s usual approach to the solution of such contact 
minimization and manipulation problems is that of inspection. The 
method of inspection presupposes a background of considerable experi- 
ence in that the designer must recognize certain contact network arrange- 
ments that may allow further rearrangements and thereby he must 
mentally develop his own rules. In order to cheek on any of his manipu- 
lations he must repeatedly redraw the network during this inspection 
design process. It is evident that this is often a long and tedious method 
and, depending on the skill of the designer, may or may not result In an 
optimum or even adequate solution. 

Suitable contact network arrangements often appear only after con- 
sideration of several alternative schemes and the rearrangements of the 
network interconnections of these schemes. Realization of this makes it 
quite e\ ident that any means ol obtaining and comparing these various 
schemes quickly and with a mathematical accuracy which does not 
require continuous checking of network paths permits a more rapid and 
complete exploration of the particular problem. Switching algebra, first 
codified by C. FE. Shannon', is the systematic application of G. Boole’s? 
~ 1G, E. Shannon, A Symb jsis of Relay and Switching Circuits, Trans 
AIEEE, 57, 1938 


2G. Boole, The at fied a s of Logic (Cambridge 1847) and An In 
tigation of the La Oud London 1854 
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nd is just such a means. [t ts 

vate the complex combinational and 

to determine sutistactoryv contact arrangements o1 

vith no minimum of time and effort It should 

with anv tool, satisfactory results dle- 

nuitv al dl logical rensoning of the user. 

ident from the following development, switch 

ts presel t stute is not to be considered entirely 

the most beneficial results, should be applied, when 
coninmetion with i spectiol techniques so that the latter 
nv limitations in the algebra techniques which have not been 


omple ted Vate miatized is vet due to the newness ot this field. 


Lh blem of solving the contact requirements of a minimum relay 


iw and translating eireuit recently developed as a com- 
vy register ol the No. 5 Crossbar Svstem will be 
practical use now being made of 

ind of indicating exaetly where the upplication of the 


desig proble mM 


COUNTER REQUIREMENTS 


funetion of the originating register is to receive pulse 
senting digits from a telephone dial or similar calling device 
i record of the digits inn torm suitable for use Hy in external 
il pulse counting and transhiting cireuit, an integral part 

is oriented with respect to other parts ol the 

im of Fig. 1. The L relay is the pulse detecting 

When the subseriber’s switehhook contact is closed due to the 
ot the pho e. the originating register is connected to the line and 
~ Oper ited Phereatter i follows the brenks snd mikes ol 

dint and feeds these re } ited dial pulses Into the counter. 

ire coul ted thre Viare ti inslated tou new code Ih switch- 
dvantageous to translate from the basic dial ten pulse 

voout of five” se lf-chee king code In this latter code 

ircuit, Will result in either one on three relavs 

storage circuit rather than two and thus an 

detected Lhe output of the translator is fed Vila a 

the register or storage circuit. The slow release RA 

Which holds between the indi 

only at the end of the pulse train. 

the translating cireuit and thereby transfers 


the storage cireuit. The RAT relay 
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in operating terminates the output from the translator and simultane- 
ously releases the relays in the counter to prepare it for the next digit. 

Specific requirements imposed by the originating register circuit neces- 
sitate the counting of one to eleven pulses; the use of a driving source 
consisting of a single break-make (or transfer) contact with ground on 
the armature spring; and outputs as follows: 

1. Count of | through 10: ground on two of the 0, 1, 2, 4, 7 output 
leads in the combination corresponding to the count. 

2. Count of 10: ground on the ZO lead. 

3. Count of 11: ground on the 0 lead only (this is a trouble-detecting 
feature). 

In addition, the design of the steering and register-storage circuit 
requires that no output leads be connected together until the second 
pulse is received. Furthermore, each relay is limited to a combination of 
simple make and break contacts not exceeding a total of twelve. This uti- 
lizes the maximum number of springs obtainable on presently available 
relays and also avoids the larger armature gaps imposed by transfers 
which would result in a reduction in the relay speed of operation. Speed 
requirements also do not permit the use of shunt release in the circuit 
operation, 





oY RELAY COILS AND |, 
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FOR PRODUCING 
REQUIRED OPERATING 
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Fig. 1 The schematic of a portion of a dial pulse register circuit for counting 
decimal code pulses and translating them to ‘‘two out of five’’ signals In the 
symbolism used in the illustrations a cross indicates a “make’’ contact and a 
vertical bar indicates a “break’’ contact 
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rTHRORY OF A MINIMUM RELAY COUNTER 


The counting cireuit unde consideration does not contemplate the 
use of any circuit elements other than relays that react to the beginning 
or end of a pulse. Therefore it must establish a distinet combination of 
relays operated or released during and between successive pulses. The 
minimum number of ordinary ‘“‘two-position”’ relays, R, required to 
count P pulses can be obtained from the expressions (1) 2P < 2” if 
the counter is to lock up during, or recycle after, the last pulse or 
a s2 | if the counter is to lock up after the last pulse. 

The usual counting circuit used for determining the number of pulses 
ina dial train is required to count ten pulses, however there are certain 


advantages in regard to trouble indications if the counter counts eleven 


pulses. In either case the minimum number of relays necessary, accord- 


ing to the preceding formulae, is five. It should be noted that the ease 
with which this minimum number can be attained depends upon whether 


the Input is derived from a single, double or transfer contact source. 


DETERMINATION OF OPERATING SEQUENCI 


Having determined that the minimum number of relays necessary is 
five, the first step in design is to develop an operating sequence pattern 
from the resulting 2° or 32 possible relay combinations. These combina- 
tions may be utilized in any order deemed desirable to obtain the 23 dis- 
tinet combinations needed to differentiate between eleven pulses (22 for 
the eleven makes and breaks plus an all relays normal combination). In 
this phase of the design switching algebra is not involved. The optimum 
sequence to meet a particular set of requirements can only be determined 
1\ repeated trials guided by ah intimate knowledge of objectives. 

Initial studies, made by Joseph Michal, of various possible sequence 
patterns for a five relay circuit, including those having a three relay 

ring” followed by two auxiliary relays and those having a two relay 
pulse divider followed by three auxiliary relays, resulted in the conelu- 
tion that the latter approach was the most fruitful. The sequence pat 
tern adopted is shown in detail in Table I. The pattern is extended 
through 12 pulses, and it can be seen that the nature of the sequence is 
such that this employs all 32 combinations of the 5 relays. Several of 
these are transient and oceur during part of a pulse or inter-pulse inter- 
val. Examination of the tail end of the sequence indicates that it will be 
simpler to design on the basis of a full 12 pulses than attempt to block 
at the end of the IL pulses specified by the requirements. If trouble con- 
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TABLE I 


SEQUENCE OF OPERATION 


g Relays 


seizure 


Ist pulse 


2nd pulse 
3rd pulse 


ith pulse 


5th pulse 


6th pulse 


7th pulse 


Sth pulse 
9th pulse 


10th pulse 


lith pulse 
12th pulse 


Total of 25 32 combinations used 
Note: 0 1s used to indicate that the relay listed at the head of the column is 
operated, and 1 is used to indicate that the relay is released, 


ditions introduce pulses beyond 12, the cireuit will without difficulty 
recycle through combinations corresponding to pulses 11 and 12. 

Table I also indicates the leads which must be grounded in order to 
provide the translations to the “two out of five” and “single lead” codes. 


The characteristics of this circuit may be summarized as follows: It 


contains only five relays which is the absolute minimum necessary. It 
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uses all 32 of its available combinations. Its control and translating job 
IS complex enough to indicate the need for a considerable number of con- 
tacts and hence the need for extensive contact manipulation to minimize 
and distribute these contacts 

It, is apparent that a great deal of time would be necessary to accom 
plish this manipulation by inspection methods, thereby indicating the 
need for an additional tool such as switching algebra to assist the de- 


signers in this task. 


\LGEBRAIC METHODS APPLIED TO CONTROL CIRCUTT 


The sequence of operations of Table [is used as the starting point 
in the application of the algebra. The exact calculations necessary to 
develop the control and translating cireuit by this means are shown in 
detail later. However, the individual steps in the solution might well be 
outlined here. First, the design of the control and translating networks 
will be regarded as separate problems. In theory these can be integrated 
together, but the resultant network is likely to be so complex that under- 
standing and maintenance of the circuit would suffer. Each of the two 
networks can be individually considered as a multi-terminal network of 
the single input type. That is, the control network is an associated set of 


contacts which connects a single ground input to the windings of five 


relays, and the translating network is an associated set of contacts which 
connects a single ground input to the six output leads. Since switching 
algebra is directly applicable to two-terminal networks rather than multi- 
terminal networks, the approach to this particular problem is of neces- 
sity somewhat indirect. 

The most satisfactory method of attack is to develop first a two- 
terminal network for each of the output paths of the multi-terminal 
network under consideration. The two-terminal networks can be ex- 
pressed algebraically and manipulated into their simplest form by means 
of the switching algebra theorems to be given later. The individual net- 
works can then be inspected carefully, either in algebraic or circuit form, 
with the objective of combining them in the most advantageous fashion. 
It will be found, in general, that the simplest network configurations do 
not readily combine and that further manipulation is necessary to obtain 
an economical circuit. It is at this point that the algebra achieves its 
vreatest utility, since its application permits the simple and rapid chang- 
ing of a given two-terminal network into a large variety of different 
forms with mathematical assurance that cireuit equivalence is main- 
tained Inspection of the networks in the several forms provides clues 
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as to the preferable combining forms and often indicates additional 
manipulations that might be desirable. 

This network development is a combination of mathematics and inte 
gration by inspection, [It is characterized by repeated trials of alternative 
forms and at no stage is there any definite assurance that the optimum 
circuit has been attained. However, the ease of manipulation provided 
by the algebra greatly enhances the probability of designing a better 
circuit than would be possible by inspection alone. In combining the 
two-terminal networks, care must be taken not to introduce “sneak” 
paths which improperly connect outputs together. The algebra usually 
offers means of introducing one or two additional contacts which permit 
combining networks and vet eliminate the adverse effects of the sneak 
paths. 

The above procedure will now be carried out in detail with the switeh- 
ing algebra theorems that are used in all the following algebraic manipu- 
lations noted at the margin by the number which corresponds to the 
number of the theorem in the complete listing in Table I. This table is 

TaB_Le I] 
SWITCHING (BOOLEAN) ALGEBRA* 


Detinitions Postulates 


Addition 1) X = OorX 1, where XY is acon 

(-+ = AND Series tact or a network 
Multiplication 2a) O-O 0 

OR = Parallel 2b) 1+ 1 
(3a) 1-1 

Circuit States 3b) 0 

0 Closed Cireuit hy l Q= 

1 = Open Circuit th) 0 


Theorems 


Sia 
Sb 
ea 
Ob 
10a) 
LOb 


y o 
+ ¥'+Z 


X 


teprinted from The Desian of Swett hing Circuits by Keister, Ritchie and 
Washburn with the permission of D. Van Nostrand Co., Ine 
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taken from The Design of Switching Circuits by Keister, Ritchie and 
Washburn*. The development of the algebraic expressions from = the 
sequence of operations table will be in exact parallel to the methods 
suggested in the aforementioned text 

The symbolism adopted in the following development is basically that 
of using the notation A for all the make contacts on the A relay, and A’ 
for all the break contacts on the A relay. Contacts or groups of contacts 


In series are related by the S\ mbol of addition E and contacts or groups 


of contacts in parallel are related by the symbol for multiplication ay 


Which may or may not be explicitly written, as in ordinary algebra. 
Therefore (A + B’) symbolizes a series contact path that is closed 
when the A relay is operated and the B relay is released, while (AB’) 
symbolizes the parallel contact path that is closed when either A is 
operated or B is released. Switching algebra includes only two numerical 
values, O and 1, with the quantity O assigned to represent a closed path 
and | to represent an open path. For the tabular notation of Table I, 0 
is used to indicate that the relay listed at the head of the column is oper- 
ated and 1 is used to indicate that the relay is released. 

\s stated earlier, the present application of switching algebra utilizes 
the sequence of operation chart of Table I. The operate and release 
combinations for controlling the A, B,C, D and E relays can be selected 
from this table by observing where each relay to be controlled changes 
state. Fon example, the operate combination for relay Dis relay combti- 
nation & and the release combination for relay D) is relay combination 24. 
It is not necessary to include the contacts of a relay in its own operate 
and release combinations. Note that the A and B relays which serve as 
iu pulse divider can be controlled solely by the L relay and contacts on 
1 and B without reference to C, D, bE. However the C, D and E relays 
are internally controlled by all five counting relays. The development 


of all these control paths uses the following abbreviations: 
(X) operating combinations for the XY relay 
(\ releasing combinations for the NX relay 
holding combinations for the NX relay 
X make contact on the X relay 


furthermore as expressed hy theorem (Ga and 6b) the negative of a 
contact network is defined as a network which Is a closed path under 
all conditions for which Y is open, and is open under those conditions 


*)D. Van Nostrand. 1951 TI e Bell Te le phone Laboratories Series 





BOOLEAN ALGEBRA AND CIRCUIT DESIGN 289 


for which X is closed. Hence h(X) may be obtained from r(X) by noting 
that ACN) is the negative of r(X). Therefore the entire control path of 


any relay can be expressed generally as 


f(N) = g(NX)IN + ACX)] 


Thus for the A relay 


Also for the B relay 

g(B) 

r(B) 

h(B) 

and 

CL. 
(L 
(L 


The schematic forms of the A and B control circuits as represented by 
the above algebraic expressions are shown in Fig. 2a and 2b. Since the 
general requirements of the basic problem specify only one transfer on 
the L relay, only simple makes and breaks on the A and B relays and no 
shunt release paths (to avoid reduction in speed of operation), the combi- 
nation of the above specifie circuits is not possible without recourse to 
double windings. Another factor which affeets the practical form of the 
circuit is the finite transit time of the L relay armature spring. Switching 
algebra presupposes Instantaneous action of relay contacts and in certain 
cases, When the use of a break-make transfer is required, additional con 
tacts are necessary to cover the open contact interval. The final circuit 


form, conceived by F. K. Low, is shown on Fig. 2¢ and uses an added A 
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contact. Algebraic equivalence of this circuit with the original is shown 
below 


i+ BL 


(L". 
L 
L’ 4 
L", 
(1, 
L’ 


For the ¢ relay which operate 4 


combination cycle 
1)’ 
y+ D+ 
A+ B’'+D 
NBDE + ABDI 


LY\) 
4)} 


[(C + A'BD'E + A'BDE") 
A'B(D' + E'\(D 4 


(3b) 
(3b) 


The schematic circuit which the above represents is shown in Fig. 3a. 


Circuits of this type which use certain contacts more than once can some- 














Fig. 2. Pulse divider of counting circuit 
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times be drawn in bridge form with a consequent saving of contacts. 
One method is to manipulate the expression into a form which is known 
to be the series-parallel equivalent of a bridge. However, following 
usual algebraic procedures it is often difficult to recognize where this is 
possible. In the present case a method developed by G. R. Frost (not 
yet published) was used effectively. This resulted in the bridge circuit of 
Fig. 3b which has the series-parallel equivalent: 
fC) = (C+ (D+ E’)\(D + BA’ + B+ (D' 4+ EF’ 

+ B] 


By use of theorem (3b) this is seen to be equivalent to the previous 


expression for f(C). 
For the D relay which only operates and releases once in the entire 
cycle 
g(D) = (A + 
r(D) = (A 
h(D) = (A 
























































Internal control of counting circuit 
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By noting that A + Cis the negative of A’B’C’, this can be 
reduced to 

f(D '+ h’)(D + A'B'C’) (3b, 12a) 

However, this transformation introduces a hazard caused by the 

transit time of A relay contacts in passing from relay combination 9 to 

10. Therefore the original expression will be used for relay D. The control 


path Is shown in hig. 3C 


Kor the & relay which operates and locks only once in the cycle 
gik A'+ B+ CC + D) 
h(h 
(FE) = (A' + B+ (C+ DE + E) 
A'+ B'+C'+ DE (da) 


This control path Is shown inh hig. 3d. 
\part from the problem of developing the required contact network, 


the practical problem of what operating power must be given to the 


relays in order to meet speed requirements must be dealt with. Since 
the use of low resistance windings in series with protective external 
resistors is called for to obtain the speed required, it appears that the 
use of two windings per relay might prove advantageous. By operating 
on the low resistance winding while locking on the high resistance wind- 
ing, the current drain may be reduced (thereby saving a fuse) and 
furthermore some code reduction may be made possible as shown later. 
If double windings are used, two of the external series resistors may be 
eliminated by combining the control network so as to make certain that 
only one of the low resistance windings on the C, D, or FE relays is 
energized at any one time. This would permit the use of one common 
external resistance with the aforementioned relays instead of three. 
Keeping these practical considerations in mind, further savings may 
be made by combining the control cireuits as shown in Fig. 3e. Although 
there is in this cireuit a possibility of contact stagger on the A relay 
contacts causing the C and J) low resistance windings to be energized at 
the same time, this will not be harmful since, when the stagger occurs, 
both relavs are firmly locked operated by: their high resistance holding 


windings 


IRANSLATING CIRCUTT 


Phe translating circuit is particularly adaptable to switching algebra 
manipulation. Table T shows the combinations which prevail at the end 
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of each pulse and the necessary “code” leads that must be grounded at 
these times. Reference to the block diagram of Fig. | shows that the 


output of the translator is not activated until the slow release RA relay 


releases after the last pulse of a digit has been received. Therefore the A 
relay can be eliminated from these combinations since at the end of 
every pulse the A and B relays are either both operated or both released 
and hence only one is needed to indicate the condition of both. Table IT] 
lists the numerous combinations which must close a ground path through 
to each of the five code leads and the ZO lead. At the conclusion of the 
algebraic manipulation, the A and B contacts may be redistributed 
evenly since they perform interchangeable functions in translation. 
The objective in the design of the translating circuit is to obtain the 
most economical contact network subject to a spring distribution that 
TasLe III 


TRANSLATION 


Counting Relays 
Output Lead 
Grounded 


ZV 


Inconsequential combinations that m: be simplification 
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fits in with the control contact network. Again, this is a multi-terminal 
network problem and the procedure is to design two-terminal networks 
that combine most readily. Since it is impractical to illustrate all the 
repeated trials that led to the final design, each network will be designed 


separately with the understanding that some of the steps are imposed 


by the form of all networks viewed collectively. 

The procedure adopted for developing the “0” lead network is as 
follows. First set up the miniature table repeating the portion of Table 
[iL that corresponds to the “0” lead. These parallel combinations should 
then be manipulated algebraically to obtain the greatest simplification 
possible. It is rather easy to apply some of the algebraic rules by observ- 
ing the condition of the relay in the several combinations in the table. 
\ simple “shorthand” rule to follow is: if in the table of combinations 
describing a particular two terminal network, all possible combinations 
of certain relays appear in conjunction with a single combination of other 
relays, the network contacts on the former relays may be neglected. In 
other words when 2" different combinations of any number of variables 
m, are identical in all but » columns, contacts on the corresponding n 
relays are not required. This procedure is carried out below. 


D 





(BB + C’ + D’ 


==’ + C + B) 





Bsc + £) 
-(C’ + D’ + E) 


Thus we have the following algebraic expression for the “OQ” lead, 


Which can be simplified as shown. 


(3b) 

(3b) 

This is shown on Fig. da. A somewhat different manipulation of the 
equation permits placing the network in the bridge form of Fig. 4b. 


he algebra equation, v1IVen below, ean easily be shown to be the 


equivale nt of the original 


C" + BB + D‘)| 





BOOLEAN ALGEBRA AND CIRCUIT DESIGN 295 


In certain cases the use of theorem (14b), normally employed to reduce 
the contacts of a particular relay to a single make and break, can pro- 
duce simplifications difficult to accomplish otherwise. This is shown 
below, with the theorem applied with respect to relay F since & tended 
otherwise to be heavily loaded. 


D’ + 0)] 
(C* 4 ))’ + 1)] 
(l4b 


T BD’) hk’ T ig: F: r} (Ya, 10a, 3b, 9b, db) 


By modifying the first factor of the final expression in accordance with 
theorem Ila, this equation can be put in bridge form as shown on Fig. 4c. 


[e+ C’ + BOB’ + DYER’ + (B’ + C)\(B + C’ + Dd) 


The above equation uses the same contacts as the previous expression, 
and although the right hand member is in a slightly different form, the 
expression Is equivalent to the one obtained earlier. 

When it is known that output conditions are inconsequential for some 
relay combinations, these inconsequential relay combinations may be 
combined with valid combinations to eliminate contacts in the network. 
Inconsequential means that the output during these particular combina- 
tions does not affect the proper functioning of the circuit. Four such 
combinations are listed in Table IIT. Only those inconsequential combi- 
nations which will combine readily with the actual combinations, thereby 
resulting in a reduction in the number of contacts, are to be used. Al 
though the use of all the all-relays-released condition may be helpful in 
certain cases, it will not be used in the cireuit under consideration since 
its use makes the requirement that no tie shall exist between output leads 
until the second pulse is received hard to meet 

With this in mind the ‘0’ lead network is agadin examined. Note the 
use of another “shorthand” rule which states that if a part of the 2’ 
possible combinations is used in closing a path, the negative of the unused 
part of the 2" possible combinations is equivalent to the original com- 
binations. Thus if in the case at hand three of the possible four combina- 


tions of the B and C relays occur in series with the same combination of 


the D and E relays, the expression used is that for the series path of the 


D and FE relays plus the negative of the missing combination of the B 
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and © relays. In the following tabulation the combinations below the 
horizontal line are inconsequential, 


+ D’ + BE) 


+ KH’ + B’D’) 


D+ 0)) 


D’ + BI(B+ D+ 1) 


l 
(14b on F) 


D' + B)} 
(9a, 9b, 10a, 10b) 


D')(B+ D)[k’ + CD’ + CB + B’'D" (Sb, 9a, 4b, 5a) 


1D) y+ J))\[k’ 4 t (12b) 
D t Ph’ + ( ' ‘(B+ C0) (13) 


Fig. dd shows the schematic of the above expression. It is possible to 
put this in a bridge form without other changes because of the manner 
in which the front and back contacts of D are related to the other con- 
tacts. Comparison of all the circuits of Fig. 4 indicates that they all use 
the same number of contacts although final decision should be post- 
poned until all the output circuits are obtained and the ease of combina- 
tion of the different circuits can be compared, 

he procedure for determining the remaining code leads is carried out 


on the tollowing pages. 
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0 
0 0 0 O (C + D + E) 





resulting in [i + C + DI|[k” + B] which is shown on Fig. 5a. It will later 
be found advantageous, in combining, to include the B’ term in the first 


factor, giving the expression: 


(+ B’' 4 " + (k’ +B shown on Fig. db 


I ee 


“S40 
Pia 


1. 

a 
0~ tania 
| dl 
1 
ma = 


————— 


hence 


C+ BCD’ + E)\C’ + B’ 
For later ease in combining, this is changed to: 
C+ BOB’ + D’' + k&)IIC’ 


shown on lig. ode. 
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+ B/C’) 


hence 


+ fb] 


whir hy Is shown Ol hig a 


7" li ad 
B 





hown On big De 


I’) which is shown in Fig. 5f 
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The final contact savings are achieved by combining the various out 
put paths. The combined translation circuit that appears to be as re- 
duced as possible is shown In Fig. 6. Note that certain forms of the 
individual output paths combine more readily than others. For example 
Fig. de and Sb combine more readily with the remaining paths than 
Fig. 4b and 5a. Note also that sometimes it is not the most reduced form 
of the individual output paths that permits efficient combining. This ts 
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Fig. 5—The “1, 2, 4,7, and ZO” leads of the translating circuit 
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exemplified by the use of Fig. Sb rather than Fig. 5a. Although various 
forms of all of the output leads were tested for efhicient combination only 
the form used is shown for outputs other than the “O" and “1 leads. 

It is essential to scrutinize the final network for possible sneak paths. 
Sometimes to avoid these sneak paths It Is necessary to add one contact 
on one relay to allow savings on others. Here again the inspection tech- 


niques go hand in hand with switching algebra and the need for both is 


obvious Phe algebra obtaims the Various torms which are capable ot 


_| 





















































is} | RAI 


| i. ( ornbine d a inslating ecireuit 








different degrees of combination very quickly and efficiently. The inspec- 
tion method is then necessary for the actual combination of these forms. 

Phe additional RA relay contact is necessary to assist in avoiding in- 
terconnections between the output leads until after the second pulse is 
received. The final assignments of either A or B relay contacts are 


chosen to equalize the loud on these relays 


rik COMPLETED CIRCUTI 


Phe final form of the counting and translating circuit is shown on 
hie. 7. The relavs are all double wound to gain the benefits of current 
drain reduction. One additional advantage of using double windings Is 
the relay code reduction made possible since now only two codes are 


necessary, One code serves the A relay and one other code serves the 
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B,C, D and EF relays. In comparison to this total of five relays and two 
codes the circuit in present use in the latest crossbar system requires ten 
relays and seveh Ce des, 


AN ALTERNATE DESIGN OF THE PULSE DIVIDER 


To illustrate the application of algebra where the apparatus contem 
plated puts less premium on contact minimization but more on stand- 
ardization and winding minimization, certain modifications of the 
proposed circuit are considered. 

In the event that new apparatus developments make possible the 
construction of relays that meet the necessary speed requirements even 
though winding impedance is increased, it appears possible af the pulse 
divider is redesigned) to use only one code having a single winding for 
all five relays. The use of added contacts might be allowable if the new 
type of relay carries more springs than the present relay. 

The redesign of the pulse divider to use single windings can be accom- 
plished by manipulation of the basic algebraic expressions derived earlier 
for the pulse divider. 

Thus for the A relay 


BA + BL + B)I lla) 


j 


By attempting to manipulate the B relay control circuit into the same 
form, one obtains 


The schematics represented by the above algebraic expressions are 


shown in Fig. 8a and 8b. The circuit of Fig. Se is obtained by combining 


the first two circuits so that only a single transfer is needed on the L 
relay. Note however that it is necessary to make the lower two B trans 
fers have continuity action to insure proper functioning. Fig. 8d shows the 
pulse divider drawn in conventional form, 


CONCLUSION 


As far as is known, the dial pulse counting and translating cireuit 
described herein requires fewer relays than any other circuit with similar 
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functions at present employed in Bell System standard switching equip- 
ment. The previous dial pulse counter used in the latest crossbar system 
required a total of ten relay s. Thus the present design represents a 
considerable saving in cost and space To a certain extent this result can 
be ascribed to the use of switching algebra during the circuit develop- 
ment 

Relay circuits designed on the basis of utilizing a large proportion of 
the possible combinations permitted by the component relays usually 


require heavy spring pile-ups. Since general purpose relays are limited 


inthe number of springs which they can carry, this type of circuit usually 


J + 4 a | 


F 
























































Fig. 8 


entails considerable design effort to make most effective use of the avail- 
able springs. Application of switching algebra to this aspect of the design 
problem can often provide crucial assistance. 

It is recognized that switching algebra, in its present state of develop- 
ment, does not permit complete mathematical statement and manipula- 
tion of multi-terminal networks as represented by the counting and 
translating cireuits. It does provide, however, facilities in manipulating 
two-terminal networks into a Variety of forms from which can be selected 
those that combine most readily. This can result not only in a saving of 
time, but also in improved circuits which might not be realized by other 
design techniques. Unfortunately the algebra in its present state does 
not indicate when the optimum circuit has been attained. To some extent 


this is caused by apparatus or circuit considerations to which, since it 1s 
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concerned solely with contact networks, the algebra does not apply. 
Thus, there is still considerable room left for the ingenuity and judgment 
of the switching circuit designer. 

As a result of the experience in designing the dial pulse counter and 
translator, certain observations on the use of the algebra are believed to 
be valid. Although switching algebra may be used in the design of the 
simplest circuits, the most noticeable benefits are obtained by the appli 
cation of the algebra to the design of those circuits in which the control 
and output paths are complex and interrelated. The particular minimum 
relay counting and translating circuit under discussion is an excellent 
example of this type of circuit. A secondary advantage of the algebra is 
its compact notation and its value as an efficient circuit “bookkeeping”’ 
method. 
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Interaction of Polymers and Mechanical 
Waves 
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New techniques of Mason, MceSkimin, Hopkins and co-workers for gener- 


alion of shear waves over the fre quency range 2 K 10° to 2.4 XK 10° cps 


have been used to study mechanical properties of chain polymers. Polymer 
solids, melts and dilute solutions, representing the main states in which 
plastics and rubbers are fabricated or used, were explored to find the char- 
acteristic relaxation times, rigidities and viscosities of various chemical 
sfruclures Polyisobut jl ne, he ved rubbe sf polydime thyl siloxane vin yl chlo- 
ride-acelate copolymers and plasticized nitrocellulose were compared with 
polyethylene and polyamides as examples of the range of solid properties 
encountered 

1s melts, several polyisobutylenes, polybutadiene, polypropylene, poly- 
propylene sebacate and poly-a-methyl styrene were investigated as models 
for varying degrees of chain substitution. Chain rigidity in, for instance, 
polyisobutylene, seemed to reflect visco-elastic over-all configurational changes 
up through the kilocyele range, but nearest neighbor interactions took over im 
the megacycle reqion, leading to moduli of 10° dynes cm” even for syrupy 
fluads 

In dilute solution, pol psobutyle ne, polystyre ne, natural rubbe rand poly- 
hutadiene microgel exhibited characteristic dynamic viscosities and rigidt- 
fics depending linearly on concentration. Pre sumably, this reflects mechan- 
ical properties of isolated chains. Some possible models were suggested for 


the fre quency de pe ndence of such prope rtves. 


INTRODUCTION 


The “equilibrium” mechanics of polymers, the giant molecules of 
plastics and rubbers, have been quite elegantly developed in the range 
of high strains (‘kinetic theory” of elasticity Mever, et al.). However, 
the molecular displacements as these strains, and, indeed, much smaller 
ones, occur, are little understood.” Nevertheless, it is essential to know 
about detailed motions in connecting chemical structure with physical 
properties. Only in this way can there be obtained from the chemical 
industry compositions which will serve properly in telephone apparatus. 


306 





INTERACTION OF POLYMERS AND MECHANICAL WAVES 307 


Other studies have treated one way of getting at these mechanisms by 
relating stress relaxation, creep, viscosity, ete. to a distribution of 
molecular relaxation times (and energy barriers), as originated by 
Kuhn.’ * Another approach is to strain polymers with periodic waves 
over a very wide spectrum of wavelengths, eventually going to tre 
quencies Comparable with those of the thermal vibrations of significant 
groups or segments in the macromolecules. The resulting dispersion or 
resonance phenomena can then be examined. Hence a mechanical radi 
ation field can interact with the masses of elementary structural units, 
as the usual electromagnetic field interacts with atomic and group 
charges. In general, direct interpretations of this kind must be done 
with shear waves, and, at least, not only with longitudinal or ultrasonic 
Waves. 

This kind of study is now proceeding using Waves generated and fol 
lowed by piezoelectric ervstals connected in as actual electromechanical 
circuit elements (A. M. Nicolson, 1919). Recent schemes of Mason and 
co-workers cover the frequency range from 10 & LO’ to 60 10° CDs, 
as reported in the paper by Mason and MeSkimin in the last issue, 
while a tuning fork method used by I. L. Hopkins has been applied to 


soft” polymers (rubbers) over the range 10° to 10° cps the general 


range of J. D. Ferry’s work at Wisconsin on concentrated polymer 
solutions). 

The relation of these studies to the scientific and technical exploita- 
tion of plastics and rubbers is in knowing what a particular chemical 
composition does to strength, stiffness, ease of molding, impact tough 
ness, ete. That is, are there qualities of the interaction of saturated 
aliphatic groups that make polyethylene or polyisobutylene have some 
glass-like as well as liquid-like, or rubbery, nature even at room tem 
perature? If so, conditions causing brittle failures must be watched for. 
How is the storage of molecular strains in injection molded plastics re 
duced by increasing molding temperature (when the kinetic theory stiff- 
ness per chain actually increases)? These and many similar problems 
may be generalized under the headings below a each case the chemical 
structure of the macromolecule appears to be reflected in relaxation 
times which combine in different ways to give flow or rigidity, toughness 


or brittleness. 


Extrusion and Molding 


Non-Newtonian flow leading to “frozen-in” stresses, subsequent dis 
tortion and irregular shapes of plastics” and rubbers,” implies energy 
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storage in the sheared molecules. The dynamic shear studies will confirm 
this. Also dispersion of carbon black and other pigments is restrained 
by elastic qualities of “liquid” polymers (.e., instead of “ mixing’, com- 
pounds just microscopically deform and later re-form.) Likewise, the 
efficieney of compounding’ and extrusion’ depend on how quickly the 


molecules relax atter straming 


I mpac i Stre ngth, By ittleness and Ts nacity 


Poughness, mechanical shock resistance, ultimate elongation and 
strength reflect the facility with which the polymer molecules can be dis- 
placed without breaking the piece. Thus, they accommodate to the stress 
hy motions presumably similar to those described above. (The situation 
is complicated when crystallites are also displaced.’) In any case, time 


10 
sensitivity in the range 10° see upward exists. 


The discussion by 
Morey is a valuable survey of these ideas, and explicitly notes the 
significance ot multiple relaxation processes on damping of shock waves. 
evidence of the relation of simple changes in chemical structure to the 
principle relaxation times effective in these physical properties of plastics 
and rubbers is thus another part of the dynamics studies. The “brittle 
point’, or volume-temperature transition of amorphous polymers, 

apparently reflects directly the correspondence of the time of experiment 


14, 15 . 
with dominant relaxation time of the polymer. A few measurements 


on plasticized polymethyl methacrylate (from which, however, no actual 
rigidities were calculated) indeed indicate abrupt stiffening as a function 
of frequency at a given temperature.’ However, the changes measured 
were too small and indefinite to indicate any particular molecular relax- 
ation. Other work’ with plasticized polymers is nevertheless concordant 
with the current findings that molecular relaxations and not long range 
order determine embrittlement. The converse of this is, of course, that 
us some “transition” is approached, hysteresis, heat build up, flex crack- 


Ing and tatigue are greatest. 


reep, Stress Relaxation and Recovery 


Keven these “long time” qualities of plastics, such as found in cold 
flow, apparently result from integrated displacements of rapidly oscil- 
lating segments of the chains. A most interesting analysis of stress 
relaxation in rubbers employs Kuhn’s suggestion of a particular distribu- 
tion of relaxation times.” The present point is that, again, these relaxa- 
tion times refleet processes which should appear directly in reaction of 


the polymer with high frequency shear waves. 
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From these aspects above the current results of dynamics studies will 
be reviewed. 


POLYMER SOLIDS: OVER-ALL MECHANICS 


Solid polymers will denote rubbers and plasties in the state in which 
they are technically used. This is usually their most complex form, with 
inter- and intra-molecular factors undistinguished. Thus, separation and 
identification of the main relaxation processes are difficult or impossible. 
However, it is interesting to consider typical values of modulus and 
Viscosity as related to chemical structure, in the range of frequencies 
corresponding to extrusion rates, and stresses in actual use. 

These values of dynamic modulus and viscosity are distinct from the 
usual quantities in the literature. The usual expressions are for longi- 
tudinal (sound) waves, and give dynamie Young’s modulus ~ 

E* ki ih, 
ky measures the out of phase part of the force-displacement relation, 
and ky = w: (“effective viscosity coefficient”). Now, the general elastic 
constants are \ + 2u, with AX = Lamé’s constant and u shear modu- 
lus. Here, 


A+ Qu = K + $y, 
with A = bulk modulus. Alternately, 
3K _ dd + 2p 
A+ u ae +p 


However, in general the present results lead to the simpler shear modu- 


ky 


lus w. Further the energy losses studied are expressible directly as the 


usual shear viscosity 


Previous comprehensive studies of the dynamics of rubbers over sig- 
nificant frequency ranges have yielded loss factors either written as 
BE. Fy (see above), or as a function of the shear VISCOSITY based on 
Stoke’s assumption that the compressional (dilatational) viscosity fs 
zero. But as Nolle’’ and Ivey, Mrowea and Guth” clearly recognize, 


recent work has strongly manifested the presence of compressional vis 


cosity in simple liquids” as well as polymeric ones Hence, the pres 


ent understanding relating molecular structure to viscosity, plasticity 
and visco-elasticity is unsuitable for interpreting mechanical wave mo 


tion more complex than in shear, unless shear constants are also known. 
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Phis sums up to mean that the chemical interpretation of basic poly- 


requires shear Wave measurements. Nevertheless, fas- 


mer mechanics 
cinating evidence of the existence of fine-structure relaxations in polymer 
come from longitudinal wave investigations. ?* 7 ™ ” Also, 

pioneering shear wave studies of Ferry and collaborators” ’ on con- 
centrated solutions of polymers have suggested intrinsic relaxations of 
the chain molecules in a highly plasticized “semi-solid” state. 

Phe more simplified findings cited below will be seen to unify ap- 
proaches in this fleld. Comment must first be made, however, on formu- 


lation of experimental results in dynamics of polymers. 


Krpression of Dynamic Properties 


Alternate and equivalent EeXPresslons have been thoroughly surveyed ; 
all represent Combinations of either Maxwell (series) springs and pistons 
elasticity and viscosity) or Voigt (parallel) springs and pistons. Obvi- 


ously, there is no physical separation ol elastic and VISCOUS elements in 


a polymer molecule, so the irrelevance of the detail of the model need 


not be emphasized. However, the models lead to convenient formulation 
of relaration times which dielectric studies, in particular, have shown 
have clear connections with chemical structure. In this chapter, some- 
times one and sometimes the other model, or combination, will be used, 
vith the symbols shown on the next page. 


Other symbols are sometimes used,” but should be easily identified 


in terms of the above 


Ie ith . and Soft Plastic S 


In Pable I, the shear modulh ol rigidity, bi, and of VISCOSITY, yw’, are 
shown as calculated for the Kelvin-Voigt model, for polymers having 
the indicated units of structure. The frequencies are from a few hundred 

few thousand eveles, hence, in the range of much technical use, 

got tires ~3800 cps) and rates of shear during processing.” ’~ Data 

are from a general study by I. L. Hopkins” of the Bell Laboratories, 
based on a tuning fork transducer introduced by Rorden and Grieco.* 
Phe strains emploved were always small, in the range 0.3 to 1.5 per 
cent; wand uw’ were essentially independent of strain, except for some 
loaded rubber stocks. The iv values clearly trace the magnitudes to be 
expected in going from the most typical rubber (hevea) to the semi- 
rigid plastics (vinyl chloride-acetate copolymer and plasticized cellulose 
\s anticipated from steady stress observations the “plastics” 


10’ dynes em”. Increase of « with frequency is also greater as 
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the “plastics” range 1s approached ; a relaxation region is implied Figs. 
1 to 4 show the dispersion of rigidity with frequency in more detail 
Especially striking in Figs. | and 2 is the small temperature dependence 
(at least between 27° and 66°C) of uw. Because of experimental uncer 
tainty, w cannot be said to be actually higher at the higher temperatures 
in accord with straight kinetic theory, but at least it is strongly tending 
that Way, as also noted for lower frequencies studies on natural rubber.’ 
Nothing like this appears for the plastics; in plasticized nitrocellulose 
the 100 cvcle rigidity decreases 10-fold from 27° to 66°C. This is, then, 
the second general dy namic quality which reflects the low van der Waals’ 
(dipole, dispersion and induction) forces in hevea rubber and polydi 
methyl siloxane, as well as their intrachain flexibility. Interchain forces 
in polyisobutylene (Butyl rubber) are low too, but barriers to flexibility 
because of sterically hindered-CH,; groups come in. Table I and Fig. 3 
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TABLE | 


shear Modulus, wu, dynes/cm? Shear Viscosity 
27°¢ , 


Poises, w’, 27°C 


» 100 S000 
5000 cycles , 
. ycles ycles 


Hevea rubber * | 6.5 X 108 350 10 


vl chloride ~92% ) -acetate ‘ ‘ 4 25 000 2.000 


plasticized by ~36% di 


Iphth ilate 


lose nitrate x ‘ . SO O00 4,500 





CH.ONO 


imphor plasti 
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Fig. 2—Viscosity and shear modulus of polydimethyl siloxane (cross-linked 
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———= VISCOSITY 


—— SHEAR MODUL 


TY IN OYNES PER CM* 


400 6 1000 2000 4000 6000 10000 
FREQUENCY IN CYCLES PER SECOND 
ind shear modulus of polyisobutvlene (cross-linked Butyl 


emphasize that thinking about the mechanics of a particular chemical 
structure must include the spatial relationships of groups within the 
Chains, as well as between them 
Phe dynamic viscosities in Table Tare also in accord with the se 
quence of structures. Their frequency dispersion again connotes varying 
relaxation processes. Natural rabber’s low inner friction, for both com- 
pressional” and shear waves”) is famous in its low hysteresis heating. 
Phis unique property is geopolitically crucial, because adequate truck 
and bus tires cannot vel be made ot any other rubber.) Indeed, it is 
striking that at 100 eveles, a prece of gum rubber has a local viscosity 
of only 350 poises. The silicone rubber gum also has high elastie efficiency, 


nd its temperature coefficient of viscosity is very low (see Fig. 2), like 


the thermal coeflicients for familiar silicone liqquids. It is exciting to 


speculate in Figs. Land 2, whether more precise measurements which 
Hopkins is now undertaking will confirm the apparently negative tem- 


perature coefficients of viscosity at some frequencies. ‘Kinetic theory 
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viscosity” arising from transfer of momentum among thermally agitated 
chain segments, does not seem to have been considered in the theory ol 
pertect rubbers. As in gases, it would require an increase of viscosity 
with temperature. 

In polyisobutylene, however, the dynamic viscosity leaps upward in 
both magnitude and temperature dependence. It should be emphasized 
that this is, again, for a cross-linked (Butyl) gum—an infinite network 
like the hevea gum, with presumably infinite macroscopic viscosity 
The striking thing is that this internal y Iscosity Is not greatly dependent 
on the network, at the degrees of “cure” used in rubber technology. For 
instance, recent studies over the frequency range 20-600 eveles, on high 
molecular weight, M, 1.2 X 10°, linear polvisobutyvlene,” give, at 
25°C and 100 cycles, uw’ = 4800 poises, although the steady flow viscosity 
of this polymer at this temperature is greater than 3 X LO’ poises.” 
Then, the infinite network (neteady tow — *) Butyl polymer of Fig. 3 
has at 27°C and 100 cvcles yu’ SOOO poses. At LOOO evcles agreement 


appears to be about the same, and is tolerable considering the several 


100,000 - ——— —— — ——!100 1000 
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per cent of compounding ingredients in the Butyl gum, and possibility 
of a small fraction of low molecular weight uncured polymer in it. 

\lso, wide variations in the degree of cure of Butyl gums were studied 
Without large changes in w’.” In this regard, the particular sample of 
big. 3 had an equilibrium swelling ratio volume swollen polymer in 
evelohexane at 25°C volume insoluble part of dry vulcanizate) of 4.84. 
his indicates an MW, value (average molecular weight between cross- 
links) of <20,0002° Actually many of the dynamic properties can prob- 
ably be found in individual chain units or segments even smaller than 
this. This is a significant point in engineering applications where plastics 
may be cured to reduce creep but where it is desired to retain typical 
“chain” properties to increase impact toughness. That is, usually some 
optimum condition for this compromise can be found. The later seetion 
on liquids will suggest that physical properites typically associated with 
chain polymers can indeed reside in even shorter chain sections than the 


V.’s observed in usual gum rubbers 


killed Polymers 


Marked effeets of carbon black and other pigments are of course 


. », 24 
familiar in both steadv and alternating mechanics of rubbers. 


Brief comment on their influence on dynamic shear properties and thus 
relaxation mechanisms involved may be directed toward plastics, also, 
however. Phus, technologically it would be desirable to load thermo 
plastics with considerable volumes of “inert”? fillers, just as is done with 
rubber. But, almost invariably strength and toughness decline, instead 
of improving, as in the rubber case. A reason for this appears in inves- 
tigations by Tlopkins when carbon black (a standard type of reinforcing 
black) was added to Butyl rubber of the sort deseribed in Table I. It is 
that stiffness seems to rise more rapidly than internal viscositv—-e., a 
given strain results in proportionately higher stress than the aecompany- 
ing infernal viscosity provides means for dissipating the stress (as on 
impact). Hence, the brittleness which fillers normally engender in ther- 


moplast 1 ~ may represent this change in MVS. yu’ balanee. Table TL illus- 


Tassie Il 


11.000 
35 000 
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trates some values for Butyl compounds. The swelling ratio (SR) for 


the compound containing 28.6 wt. per cent filler has dropped to 3.2 
implying also considerable reductions in MW, (since theoretically 


(] SR)” as nM , ‘ Thus, the apparent chain segment between 


a 
M, 
cross-links is shorter than in the unfilled stock (the two were cured to 
give closely similar degrees of primary valence cross-linkage) and corre- 
spondingly the steady-pull modulus is higher. Yet, the internal friction, 
while also higher, seems to reflect new relaxations from interaction with 


the filler, and total shock-absorbing power has declined 


M icrocrystalline Polymers 


The preceding studies at comparatively low frequencies indicated (1) 
magnitudes of shear rigidity and internal viscosity characterizing rub- 
bers and soft plastics. By familiar shifts of temperature or frequency, 
they would also apply to polymers known as hard, amorphous plastics 
at room temperature such as polystyrene and polymethyl methacrylate. 
(2) Dispersion of w and yw’ with frequency affirm that the intrinsic or 
fine structure relaxations have times <10° to 10 * see, and so refer to 
chemical units much smaller than the average molecules in the usual 
technical rubbers and plastics. A way to get at what sizes and habits 
these units might have will be by investigation of low molecular weight 
polymer liquids. But, while still in the seetion on solids, it is recalled 
that microcrystalline polymers such as polvethylene, polyesters (Teryl- 
ene), polyamides (nylons), cellulose esters, polyvinylidene chloride, poly- 
acrylonitrile ete., have mechanical properties dominated by their crystal- 


40, 41 rp 
Che amorphous volumes are clearly those 


line-amorphous ratios. 
which donate the flexibilitv, toughness and shock-resistance of these 
plasties and textile fibers.”’ “’ An interesting point is, how “viscous” are 
the disordered chain segments? In an over-all sense, all kinds of dissipa 
tion including crystallite friction, analogous to grain friction in metals, 
scattering of longitudinal waves, and stiffening by low temperatures can 
occur In these polyphase systems. Thus, effects of chain orientation as 


well as Jateral order crystallinity ) have been detected in dynamics 


‘ 26. 41 rev . 
studies. The intrinsically bquid-like or amorphous components of 


this behaviour—and the things which will correlate most simply with 
dipole concentration and other chemical features— are most accessible 
to study at very high frequencies. For, in these polymer solids, unlike 
the essentially continuous and homogeneous amorphous ones first dis 
cussed, the mechanics reflect small regions having widely divergent 
properties. Thus, methods developed by I tS MesSkimin ol Jel] Tele 
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phone Laboratories (deseribed in the last issue), have been used to probe 
for elemental reactions at the upper end of the frequencies presently 
available. Both longitudinal and shear waves were used. In polyethylene, 
a Wavelength for the shear waves was 0.0074 ¢m., at f 8.55 xX 10' 
eycles, and in polvhexamethvlene adipamide (the usual 6-6 textile rv- 
lon), the shear wavelength was 0.0125 em., for f 8.67 & 10° eveles, 
all at 25°¢ 

The important consequence of these experiments so far has been that, 
despite the small strains involved, the VISCOSITY appears to be a  poly- 
mer’ viscosity, rather than an inner friction involving Just a few liquid- 
like atoms per unit. Thus, polyethylene of “equilibrium” crystallinity 
and average molecular weight corresponding to an intrinsic viscosity in 
xvlene of {7 O.S9 (at SO°C), was measured over the range from 0 to 


50°C. The results from both longitudinal and shear wave measurements 


TABLE III 


are given in Table IIL. These viscosities are expressed in this case for a 
Kelvin-Voigt model, of rigidity and viscosity in parallel. The rigidities 


associated with these viscosities are about 3 X 10° dynes, em’, or not far 
from the value under steady pull of about 1 X 10°. 
Now this suggests that the rigid plastic polyethylene retains, even 


under mechanical impulse of microsecond duration, a shock-absorbing 
capacity reflected in a shear viscosity of 5-15 poises, and a compres- 
sional viscosity of 3-8 poises. The former, u’, may roughly correspond 
to the liquid viscosity of a paraffin-like chain of from 50 to 65 c-atoms 
in length. Thus, the dynamics measurements seem to relate to basic 
premises of polymer structure. These are that the amorphous regions 
Whose existence is shown quite independently by X-ray scattering, den- 
sitv, heat Capacity, ete indeed take up and dissipate sudden stresses 
Which the microcrystallites, despite their great strength, would be too 
brittle to sustain 

Phese results give hope that further probing of the dynamics of liquid- 


like elements in rigid plastics will eventually lead to precise adjustment 
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of molecular weight, chemical structure (degree of branching in poly- 
ethylene), crystallinity, etc. These quantities, when fitted to a given 
pattern of w, A, vw’ and X’ at proper frequencies would yield plastics of 
optimum serviceability under the multitude of stresses encountered in 
use, 

A similar liquid-like structure-even where the (crystalline) rigidity is 
much higher” and mobile chain segments smaller—apparently oceurs in 
polyamides. Presumably the hydrogen bonding and dipole interactions 
are very imperfect in the disordered regions,” and there the chain inter 
action is reminiscent of polvethyvlene. For instance, in polvhexamethy1 
ene adipamide, measurements in the 8 to 30 megaevele range do indicate 
that the Lamé elastic constant \ is about 5.6 & 10° dynes em’, but 
only about 3 xX 10° for polvethvlene This reflects over-all stiffness 
dominated by erystallites. Nevertheless, the compressional viscosity, \’ 
is 17-6 poises (going from 8 to 30 me) for the polvamide, hut only 5-2 
poises for polvethyvlene. Of course, since there is dispersion in both cases, 
these relative magnitudes might be quite different at some other fre- 
quency or temperature (all above are at 25°C). Yet it remains that the 
nvlon, despite its hardness, also has a liquid-like component more vis 
cous than that of polvethvlene. Similar relations appear in the shear 
VISCOSITIES, w’, also determined for these two svstems. For the 6-6 poly 
amide, yw’ goes from 19 to 7 poises over the 8 to 30 me interval while 
polvethvlene changes from 15 to 5. These quantities indicate again, as 
with the polvethvlene, that polymer liquids” rather than just a few 
small groups of atoms are the important mechanical elements even at 
frequencies of 10°. Now polystyrene, an amorphous polymer, also has 
rigidities of about 10 dynes em* but the vw’ and X’ values at room tem- 
perature are far below 5 to 20 poises, and glass-like brittleness (although 
not so bad as silica glass) results. 

So far, then, the two characteristic extremes of polymer mechanics 
have been discussed: (1) the soft rubbers, whose dynamics at low kilo 
evcle frequencies imply, at ordinary temperatures, predominantly over- 
lapping combinations of relaxation processes Whose relaxation elements 


involve many segments per molecular chain; and (2) the hard, micro 


ervstalline plastics, whose behav hour Is predominated by relaxation proc 


esses having times of 10 "to 10° see because the longer period (slower 
displacements have been relaxed out at the temperatures of normal use 
(Likewise, interconvertability by temperature between these two ex 
tremes is presumed. Also, a certain correspondence between dielectric 
and dynamic relaxations in these classes is indicated."") Next, it is in 
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teresting to see what are the simplest structures (particularly in terms 
of molecular weight) vielding these effects. In other words, what kind of 
liquid really exhibits “polymer mechanics?” No detailed answer to this 
can be given below, but results on some polymer liquids of low average 


molecular weight will indicate that the mechanisms in rubbers and plas- 


tics are probably more general than previously supposed. 


POLYMER LIQUID MECHANICS 


By techniques described in detail elsewhere’: ™ a series of polvisobutyl- 
ene liquids have been investigated. These polymers were made by ionic 
catalyzed mass polymerization at reduced, but not very low, tempera- 
tures. While no great care to purify the monomer was used, such poly- 
merizations require fair purity to go at all. Seemingly, the resulting liq- 
uids do represent a polymer homologous series, although head-to-tail 
sequence of the monomer units, some single ethyl rather than paired 
methyl side groups, ete., may differ slightly from the higher molecular 
weight forms in Butyl rubber and polvisobutvlene gum. Whatever are 
these details, it appears that the polymers represent a linear hydro- 
carbon chain, with essentially two methyl groups on every other chain 


atom 


polvethvlene, with the nominal chain CH.-CH:s-CH 


ou lesser extent polystyrene, 


CH.-CH 


have chains in which rotation about the bonds is less sterically hindered. 
V hae final sec tion, on isolated polymer chains (1 dilute solution), will 
cousider this aspect further. However, some results will be reported be 


low on a low molecular weight poly-a-methyl styrene, which may be 


consicl re a tru turally i CTOSS between the rubber, polvisobutvlene, and 
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the plastic, polystvrene, 


% 


Other studies in progress on liquid polybutadiene, polvisoprene, poly 

propylene, and polypropylene sebacate from which turther information 

about intra-chain stiffness may be derived, will also be noted 
Properties of the polvisobutyvlenes studied are summarized in ‘Table 


IV, some additional molecular weights in this range appear as extra 


points in some of the high-frequency graphs. The molecular weights MV, 


° ° . or) 4 ila 
are “intrinsic viscosity’? averages and, with reasonable estimations 
M, - 
ol the V ratio, cheek with CPV OSCOPIC number average, M, P values on 


such materials, which are in turn listed in the table as expressed by melt 


_ P . 7 { ore 
VISCOSITY relations of Fox and Flory. These molecular weights repre 


TaBLe I\ 


dO F O.37 | 3 x 10 <x 1 

1660 7 30 x 10° § X Ww 

10” x WO 

i) 2520 = =1070 216 : x 10" 5.5 aA.: x 1 
56 3350 = 1720 737 ‘ x 10° 17 x 10! 
74 H7O 2530) 1840 5 x 109 : LW 
147 S240) 4850 4600 5.3 « 10° ; 2 x 1 


sent reasonable averages rather than absolute values for these hetero 
geneous polymers. The DP values are just the number of isobutylene 
units per average chain. The 7 values are the steady flow viscosities at 


low rates of shear—usually determined by a falling ball. 


Rigidity and Viscosity Magnitudes 


The properties of these liquids ranging from polymer A having only 
forty times the viscosity of water to FE, which begins to approach fluid 


ities of technical poly mer melts (polvamide s, for instance), were explored 
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in the kiloevele range with shear waves generated by torsional crystals 
and in the megaevele region by shear waves with the reflectance method 
and by longitudinal (ultrasonic) waves from a pulse propagation tech- 
nique. The results have been expressed in two ways. First, in earlier 
reports a trend corresponding with experiment was given by two 
Maxwell elements arranged in parallel This result: is too simple com- 
pared to the distributions of relaxation times previously proposed for 
high molecular weight) polymers to reproduce detailed observation. 
Nevertheless, perhaps because of the smaller molecules involved, there 
seems to be clear indication that two prinezpal relaxations predominate 
the mechanical reactions of these liquids over the range of frequencies 


of present interest, LO” to 10° eps. For example, for polymer D, these are: 


Second Relaxation 


5 ™& 108 eveles 
10° dynes/em? 


In accounting for the second main relaxation, a hysteresis component 
had to be introduced whose significance has been suggested.” ) 

Second, specific values ol shear rigidity ai Maxwell) and M (Voigt i. 

var viscosity w’ (Maxwell) and w’ (Voigt) as well as the constants for 
related compressional Wave systems, A + 2y (elastic) and Xo -+ 2p’ 

viscous) have been calculated for particular frequencies. Unlike in the 
first way of expression, these latter quantities are all highly Trequencey 
dependent. However, they deseribe conditions at various frequencies of 
interest, and are thus often worthwhile. 

Both ways of looking at the data lead, as implied by the figures above, 
to the proposal that typical polymer stiffness (shear rigidity of ~10' 
dynes em’) is present at M, ~ 1600, with DP, ~ 30, or an average 
chain length of about 60 carbon atoms. This appears when the straining 
is done in 10° to 10° see. In the 10°" to 10°” see range, rigidity occurs 


4 
loreveh an average chain length of 20 atoms as shown in Table IV. 


STRUCTURAL FACTOR IN LIQUID MECHANICS 


Phe main relaxations in the kiloevele range in polyisobutylene liquids 
seem to lead to quasi-configurational elasticity. This is where the kinetic 
theory tendency for a most probable separation ol chain ends is retarded 
by viscous interaction of seements between and within the chains. Hence, 
the middle dashed curves of Fig. 5, showing shear elasticity for some ot 


the polymers of Table IV, decrease exponentially with increasing tem- 
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perature. While pure kinetic theory elasticity would give a modulus 
increasing linearly with increasing temperature, these systems, like all 
practical rubbers and plastics, actually grow softer with rising tempera- 
ture when deformed dynamically. It is striking, nevertheless, that a 
modulus of ~10' dynes cm” seems characteristic of the visco-elastic 
energy storing of these simple polymer structures. As noted below this 
is 10° less than the erystal-like, close-packed, stiffness found for these 
sume molecular frequencies above their second principal relaxation time. 
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Thus, it seems that the former, 10’, modulus is typical of the struc- 
tural arrangements in polymers said to be above their second order 
transition temperatures while the second, 10°, modulus reflects in- 
teractions below the freezing-in. 

These conclusions obtain regardless of the particular expression of the 
data. But, for comparison, curves are shown on Fig. 5 for a polyiso- 
butvlene A” in which the dynamic modulus » at 20 ke is computed for 
both Maxwell and Voigt elements. The two points denoting the steady 
flow viscosity of polymer A” rank it with respect to the others in the 
series 

\pparently even very fluid polymer melts, chain molecule plasticizers, 
and small segments of long molecules must be expected to show appre- 
ciable rigidity when stressed rapidly. Referring to the introduction it is 
reasonable that rough extrusions, frozen-in molding stresses and the 
like are so easily produced. The lines of Fig. 5 are not, of course, implied 
to be linear over any considerable temperature range. In the region rep- 
resented, the temperature coefficient for viscous flow is about 16 keal 
for the B,C and D liquids (about 12 for A). This agrees roughly with 
the steady flow values found for very high molecular weight polyiso- 
butylene.’ *" The temperature coefficient for the rigidity is less, as 
would be expected, since the whole center of gravity of the chain need 
not be displaced, but only local segments 

This quasi-configurational elasticity is increased by molecular weight 

although kinetic theory elasticity of chain segments in a network is de- 
ereased by increasing segment length). The log uw vs density at 25°C 
plotted in Fig. 6 indieates that the chief influence is the number of 
chains per ec, since the points for all the molecular weights now lie on 


a single line. It should be repeated that the elasticity modulus plotted, 


wu, is again for a roughly frequeney-independent or “absolute” model.” 

he same is true for the three solid lines on Fig. 6, showing w in the 
second, o1 10’ evele, relaxation range. Here effects of detailed liquid 
structure come out; the three average molecular weights no longer lie 
so nearly on a single line. This elasticity is presumably from the erystal- 
like interaction of nearest-neighbor segments. If temperature is adjusted 
so that densities are the same, it is seen that the lower average molecular 
weight liquid has the Adgher elasticity modulus. This difference 1s not 
large, and should not be interpreted as showing an equal segment inter- 
action, for a polymer of lower specific volume (B compared to D). 
Rather, it emphasizes in this relaxation range, approaching the “glass” 
behaviour, that the relaxation rate is vastly more temperature dependent 


than the specifie volume change alone, and structural variations in the 
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Fig. 6—Principal shear elasticities of polyisobut ylenes as related to their densi 
ties 


packing of segments in the liquid (coordination number, ete.) become 
important. 


sip Soft” and ‘‘ Hard”’ Liquid States; Second Order Transitions 


A recent noteworthy study” of volume-temperature and viscosity 
temperature changes in polystyrene (with a note on polyisobutylene 
brings out many points in common with ideas of polymer liquid struc 
ture indicated by the dynamics work. Particularly, the fact that 
according to steady state measurements, the ‘local configurational ar 
rangement of the polymer segments’ below 7, remains fixed accords 
with the postulations from dynamics work. That is, above the second 
main relaxation, it seems to be just the interactions in these fixed ar 
rangements which cause the glassy (or “erystalline’’) dvnamie modulus 
of 10° to 10" dynes em’. Further, the point that 7 is not an isoviscous 


. 4 
state for polymers” agrees with the dynamics result that macroscopic 
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viscosity of the polymer has relatively little to do with the actual values 
of dynamic viscosities. These would be at frequencies where the response 
of the polymer liquid to the mechanical field is determined only by 
motions within the local fixed arrangements mentioned above. 

7 illustrates this, where on one scale the macroscopic viscosity is 


big 
plotted according to the familiar log-log relation with molecular weight. 


‘Two extremes of average molecular weight, VW, and M, are used for the 
liquids, to show that the molecular weight distribution does not alter 
the general conclusions. (M, is an upper limit weight average figure.) 


On the other seale, the dynamic viscosity u’, in this case for a single 


element frequency-dependent Voigt model, shows low values and marked 


curvature. These betoken the relaxation in which molecular weight, 


through its effect on free volume and other structural factors, is signif- 
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icant for displacements superficially quite different from those in macro- 
SCOPLC VISCOSITY. 


The compressional VISCOSILYV, ' is also plotted, for the same model, in 
Fig. 7. It is, within experimental error, zero for polymer A’, as deter 


mined by shear and compressional wave studies at 8 me frequency 
This is a rare case, then, where the attenuation of sound waves through 
a liquid has been quantitatively accounted for by the shear viscosity. 
But, as soon as the average molecular weight rises to LOOO or so, \’ 
comes in clearly, and the new mechanism for dissipating compressional 
or dilatational stresses is developed. As this presumably represents di 
rectly free volume or coordination number changes in liquid strue- 
ture,” its detailed study near 7, and in connection with brittle 
points of rubbers, may eventually be especially fruitful. 

Another depiction of influence of average molecular weight in these 
liquids on dynamic viscosities occurs in Fig. 8. Here, the X- 
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Fig. S—Dynamie viscosities of polvisobutylene liquids as function of 
molecular weight (25°¢ 
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for, again, a crude model attempting to show compressional viscosity 
over the whole jrequency range, while the other viscosities are Voigt 
expressions at 8 (or 14) me. Extremes of molecular weight averages are 
shown 

Comparison of the “soft” or quasi-configurational rigidities, expressed, 
like the w of Fig. 5 as relatively frequency independent yu, , with the 
“hard” or glassy rigidities is given in Fig. 9. The \ and w values are for 
the Voigt model at 8 mc. The graph does not show the bend-over of the 
‘soft’, we, curve with molecular weight, but that happens more grad- 
ually. The “hard” rigidities \ and w quite readily show this inflection. 
As before, the relaxing segments must be <100 chain atoms, according 
to the behaviour of the molecules at room temperature. 

Concerning influence of molecular weight on engineering ‘brittle’ 
points’ of such importance in rubber technology, the present studies 


agree With earlier proposals. Thus, although the 7’, or v-7 second order 


: . ' 
transition point always decreases with decreasing molecular weight, 
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the brittle pot tends to increase as the molecules get smaller. This was 
supposed to be because the ultimate elongation (doubtless due to viseo 
elastic or quasi-configurational elasticity and not kinetic theory elas 
ticity as sometimes said) declined with chain length, so that the speci 
men broke at lower and lower strains,” even though it was not really In 
the glassy stute. Now, the results above demonstrate that the shear 
modulus at a given temperature does fall off for average molecular 
weights of polvisobutvlene below ~5000. Hence, the mechanical em 
brittlement of the low molecular weight samples is not because they are 
stiffer, but because they are weaker 


POLY-@-METHYL STYRENE: A “PLASTIC” LIQUID 
Most of the liquid studies have been on polyisobutvlene polymers, 
made “hard” or ‘‘soft’” by temperature or frequency, but under use 
TABLE \ 


w’, 14 me, Poises uw, 14 me dynes/em 
Poly-a-methyl 
styrene 
Maxwell Maxwell 


111.4 


Ht 2 


conditions, considered rubbery. If one of the methyls In polyisobut ylene 
is replaced by a phenyl, poly-a-methyl styrene, a hard plastic is pro 
duced. Low molecular weight polymers of this composition are, however, 
liquuids at room temperature. Hence, it is interesting to compare their 
reaction to mechanical waves with that of polyisobutylene liquids of 
similar macroscopic viscosity. Table V lists a few properties at 25°C, 
Polymer I has roughly the steady flow 7 of polyisobutylene B. Also, 
the Voigt w’ at 14 me is similar: 15.2 poises compared to the 23.6 poises 
of the polv-a-methyl stvrene. However, the Voigt u is already 50 times 
higher for the phenyl substituted chain. (This shows the shift of the 
second relaxation range, of course, where nearest neighbor interactions 
rule.) Even more striking, the temperature coefficient of a ue , calculated 
as the second principal shear viscosity in a frequency independent 
model,” as before, is about 24 keal, compared to about 12 for Le tor 
polvisobutvlene in its similar model. Thus, although there was no ap 
parent difference between the mechanical properties of these two poly 
mer liquids in their room temperature state, their dynamics diverged 


remarkably. This was when they were studied with shear waves whose 
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lrequencies approached the glass-into-rubber relaxation times. Clearly, 
again, individual interaction of chain-like chemical units and not any 
micellar or other special aggregation of them, predominates polymer 
mechanics 

It still remains, however, to separate interactions of the basic units 
within and between chains. Most likely, the model plastic vs rubber 
liquids just discussed differ in the high frequency region substantially 
only because of inter-chain forces between phenyl vs methyl groups. 
However, especially in the high-frequency region, questions of intra- 
chain structure, such as the steric hindrance of adjacent pairs of methyl 
groups in polyisobutvlene, restricted rotations about bonds, ete., come 
in. Obviously where configurational or quasi-configurational displace- 
ments are important, as in all cases of elongation >20 per cent (this is 
certainly an upper limit), flexibility of single chains needs to be under- 
stood. This is built deeply into chemical structure; plasticizers pre- 
sumably may change over-all configuration as well as modify interaction, 
but they are impotent to vary flexibility. Accordingly, problems of rub- 
ber, usable in the Arctic, and of wire and cable insulation bendable at 
low temperatures always come back to whether the polymer chain bonds 
have free rotation, Some examples of the combinations of effects within 
and between chains can indeed be shown in several other polymer liquids 
which are rubber models. 

This influence of small changes in chemical structure is compactly 
illustrated by comparing a few other hydrocarbon polymer liquids with 
polyisobutylene. Also, rather dilute dipolar groups have been introduced 
in the linear polyester liquid polypropylene sebacate, whose structure 1s 


otherwise like that of hydroearbons.’* In Table V1, liquids of the given 


structure with some (unknown) distribution of molecular weights, were 
studied with shear waves at 77 and 142 ke at a temperature where each 
had the same steady flow viscosity. The figure chosen was 700 poises, 
and the temperature range required to adjust to it in the series was 
10.9° to 85°C, meaning that the liquids had comparable consistencies at 
ordinary temperatures 

Despite these similarities under steady stress, the retardation times, 
r’, vary three-fold, with the highly substituted hydrocarbon chains, 
polvisobutvlene and polypropylene, the highest. Despite the intermo- 
lecular action of the dipoles in polypropylene sebacate, the low polymer 
has a short retardation time, although its “brittle point’? with decreas- 
Ing temperature is far above that. of polybutadiene or even polviso- 
butvlene. Presumably the flexibility around C—O-—C bonds rather 
compensates for mecreased dipole interaction. Where both low polarity 
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and chain flexibility obtain, as in polybutadiene and the silicones, dy- 
namic properties apparently accord with brittle points in implying small 
temperature coefficients of relaxation times. In fact, the temperature 
coefficient for dynamic viscosity of polybutadiene is only about 1.5 keal, 
whereas a comparable figure for polvisobutylene and polypropylene Is 
12 keal 

The frequency range in which the structural comparisons above were 
made, resides, as discussed earlier, in the zone of configurational visco- 


elasticity. That is, over-all shape changes, rather than just nearest 


neighbor interactions, are predominant even at these comparatively 


short average chain lengths. Now, other recent studies of polvisobutyl- 
ene liquids, at 5 to 100 eps frequency, exhibit no rigidity at 25°C and 
above, although they become non-Newtonian rapidly as temperature is 


PABLI \ I] Shea Dinan ( ol Poly sobut jle ne i] 


u. 1660 1,” 39.6 Poises 


Maxwel 


266 : x 10 39.2 , 10° 
160] ~ 10 38.4 ; 10° 
25300 54x 10 19.9 ( 10 
11390 5 xX 10! 19.9 10 
53060 5 & 10 Is. ] 26x 10 


reduced.) “Phe questions are, where does the configurational elasticity 
drop out, as frequency is reduced at 25°C; and does it seem reasonable 
that this dispersion correlates with a shift in frequencies at lower tem- 
peratures. Partial answers are given by very recent studies of I. L. 
Hopkins of Bell Telephone Laboratories. He has equipped the tuning 
fork vibrator described earlier” with two parallel vanes filled in between 
with a film of polymer liquid. Pure shear properties can be derived from 
the response of this svstem. Table VIE lists a few typical figures obtained 
on polvisobutylene polymer A”. These indeed show that the kilocyele 
relaxation zone (some new data by MeSkimin’s torsional pulse method 
are given for it) extends smoothly down to where dynamic and steady 
stress viscosities are equal. Seemingly there are no new “extra long 
time’ relaxation mechanisms; probably the slow relaxation times some- 
times indicated for high molecular weight rubbers are just displacements 
of this configurational relaxation to long times because of high molecu- 
lar weight and internal viscosity. 


By contrast to the conclusions associated with the data of Table VII, 
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some observations at low Trequencies On ISOVIscous properties of polyiso- 
butvlenes A” and C indicate nearly identical retardation times. Thus A” 
at 25°C and C at 61°C have n 39.6 poises. The 9 values at 266 and 
1600 cps are also about 39 polses, w at 266 eps is 3800 dynes per em 
for both liquids, and at 1600 eps is from 3.5 X 10’ to 48 & 10° dynes 
per em. 

In the final section, mechanical waves have been used to explore 
dilute polymer solutions, to see how isolated molecules behave, free of 


interaction with each other. 


DILUTE POLYMER SOLUTIONS 
Physical Principle Sin Measurements 


Precise information on dynamics of solutions approaching infinite 
dilution (and thus complete separation ol the polymer chains) is desired 
here. Again these must be shear dynamics; bulk rigidity of ordinary 
liquids is so high that a few polymer molecules added cause little effect. 
Dilution is emphasized because even at 1 per cent by volume, high 
polymer molecule coils frequently interact, especially in“ good” solvents. 
Thus, several workers have detected shear rigidity In polymer solutions, 
in one case for polymethyl methacrylate of average molecular weight 


320,000, at 1 per cent concentration in o-dichlorobenzene.”’ Very low 


frequencies used (~10 cycles) there and in an earlier study” suggest, 


however, that even here, appreciable entangling of the molecules created 
a temporary network such as studied by Ferry.’' ~ Such was certainly 
present in the 5 to IS per cent solutions of cellulose acetate in dioxane 
measured in one of the earliest observations of shear rigidity in polymer 
solutions.” 

Accordingly, since strictly linear, and hence non-interacting, mechanics 
are sought for the macromolecules in dilute solution, careful evaluation 
of experiments Is essential. Since already it appears that important 
over-all (quasi-configurational) relaxations occur for, say, polvisobuty1 
ene in the kiloevele range, and it is suspected that not all of the inter 
actions involved are between chains, the torsional crystal techniques 
are attractive. The absolute Viscosity ol these solutions is very low, so 
the ammonium dihydrogen phosphate crystal whose piezoelectric qual 
ities are appropriate for polymer liquids in the circuits previously 


noted” ~~ Is advantageously replaced by quartz 
Detailed electromechanical behaviour of such crystals in the pure 
liquids evclohexane and benzene is of first concern. The electric field 


applied to electrodes on the suspended ervstal produces mechanical 
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torsion generating pure shear waves. These waves may be modified by 
the environment around the erystal (vacuum, gas, liquid, solid) and 
react back. Thereby a mechanical resistance, Ry, and a mechanical 
reactance, Ny, are imposed on the electrical properties of the crystal 


element in the circuit. This conneetion Comes out as: 
AR, AR Vu 
Af . kK X Ww, 


where AR, is the increase in measured electrical resistance of the crystal 
element in the medium compared to in vacuum (or practically in dry 
air or nitrogen). The decrease in resonant frequency of the crystal ele- 
ment under these conditions is Af. Thus A, and A» are electromechanical 
constants, Which fundamentally may be calculated from the dimensions 
and piezoelectric constants of the cry stals.”” Now, in simple, Newtonian 
liquids, 


Ruy = xX ye w*\ ainp 


Thus, by carefully measuring AR, (or Af) on a liquid of accurately 
known density p and viscosity 7, at a given frequency f, and a given 
temperature, the constants Ay and Ay may be evaluated without as- 
sumptions and approximations of deriving them. Their constancy will 
then reflect the electromechanical stability of the system. Their be- 
haviour under various conditions will be illustrated below. One further 
point is that when a liquid or solution does exhibit shear rigidity, or, in 
other words, if the single large molecules in a dilute solution are able 
to store energy, then Ry > Ny. Hence, in this case, the observed 
quantities ARe , and especially Af require particular precision. 

In this regard, typical magnitudes of change of fy between dry air 
and pure evclohexane, at various temperatures, appear in’ Fig. 10. 
Questions often arise as to the arbitrariness of suspension of the radi- 
ating crystal, by the fine supporting and lead wires. The effects with 
the plain wires, in the solid curves of Fig. 10, are somewhat, but not 
radically, changed when a metal bead is put on, heavily loading vibra- 
tions in the wires, as shown by the dashed curves. In Fig. 11, a some- 
What larger influence of the loaded support wires is shown for the Rx 
values, but both curves, by their smoothness and shape over a tem- 
perature range Where the thermal expansion and other elastic constants 
of the metal support wires are quite different from those of the quartz 
eryvstal, affirm reliability of mounting and electromechanical coupling. 

Fig. 10 shows, even for an 80-ke erystal, that Af for an organic liquid 
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Fig. 10—Temperature variation of resonant frequency before and after add 
ing weights to mounting wires 


(or dilute polymer solution) is bothersomely small. An excellent oscillator 
at 20 ke can hardly be expected to drift less than +2 cycles, but at 20 
ke the Af like that between the sets of curves on Fig. 10 might be only 
10 cveles, so 20 to 35 per cent error could come in. Hence, a different 


? AP - 93 49 
scheme for measurement of f, than that in earlier systems Was 


evolved. The tenth harmonic of the (say 80 ke) resonant frequency was 
beat against the 79th harmonic of a controlled standard 10-ke frequency. 
An interpolation oscillator accurately readable to 1 cycle then supplies 
the many hundred (roughly 1000) difference between these two high 
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Fig 1] Temperature variation of resistance at resonance before and after 
adding weights to mounting wires, 
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ter adding weights to mounting wires 


Temperature variation of erystal 


harmonics In this Wily, and In about 350 sec a balance can be conven- 


iently achieved and the required ten-fold gain in accuracy attained, 
By these means, and with best literature values of viscosity and 
density 


which were checked in the laboratory at several temperatures ) 


for purified solvents, curves for AK, and Ky. were obtained as exhibited 
in Fig. 12 for 8O ke 


Behaviour of A, at different frequencies over a tem- 
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perature range is shown in Fig. 13, and of Ae, in Fig. 14. Fig. 14 brings 
out the significant point that in the present arrangement, where the 
oscillating crystal is immersed in the liquid studied, the dielectric prop 
erties of the liquid are important. Apparently the dielectric losses even 
of these purified hydrocarbons are different enough so that A» at 80 ke 
is quite different for benzene and cyclohexane. (Dielectric studies have 
previously indicated difficulty in preparing benzene having theoretically 
expected loss.) It is also possible that slight differences in wetting the 
crystal cause Ay to vary with the liquid used. 

The kK, and Kk, values determined for all the various conditions above 
were then used under these conditions for measurements on the poly mer 
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Fig. 14 Temperature variation of crystal constant A. over a frequency range 
with benzene and cyclohexane as standard fluids 
solutions in the kiloevele range. In the megacycle range, the balanced 
shear wave reflectance technique” gave satisfactory results over certain 
concentration zones which could fairly well be extrapolated to high 
dilutions. Thus, over the whole spectrum, there seems to be no doubt 
about the reality of the effects described below. That is, their magnitude 


far exceeds experimental uncertainty, as demonstrated in this section 


POLYISOBUTYLENE SOLUTIONS; DYNAMICS OF SEPARATE CHAINS 


Solutions of polvisobut vlene ol M. La x 10 from about O.1 to 


1.0 wt. per cent concentration in evcelohexane vield R vu and w Curves 


as shown in Fig. 15. The points coincide for the pure solvent, “as they 
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CONCENTRATION IN GRAMS (PER 100 ML. OF SOLUTION) 


Fig. | Klectromechanical interaction of solutions of polvisobutylene (M 
1.18 & 10%) in evelohexane with erystal vibrating torsionally at 20 ke 


should for a liquid having only viscosity. But, apparently as soon as any 
polymer chaims are added, the curves diverge. A stiffness coming from 
separate chain molecules is being displayed.” Qualitatively, theoretical 
expectations ol Kuhn and others seem justified, at least that there 
isa relaxation mechanism for isolated chains. 

The usual question of how best to express the dynamical results arises 
The procedure of earlier sections for polymer solids and liquuids will be 
followed. In general, a frequency dependent modified Warwell element 
as sketched on Fig. 16 will be used. However, a frequency-independent 
analy sis has also been carried out for one sumple svstem, and, from this, 
basic mechanical constants of single “average”? molecules are obtained, 
if it is reasonable to relate the mechanical models for the liquid con- 
tinuum to the discrete chains dissolved in it. 

hig. 17 shows typical results from the simple scheme of Fig. 16, where 
the pure solvent viscosity, ny, has been considered to be in parallel 
with a Maxwell element. The total shear rigidity of the solution (at a 
given concentration) is represented by un. The y ISCOSILY of the polymer 


molecule coils in solution with the solvent streaming through them is 


x < 


Wp) WT). 





Relations for ealeulation of shear stiffmess and viscosity of dilute 


SOLUTLOnS 
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taken to be Ne Thus, the steady flow viscositv, n 


/ 


Na tT Ne 1B 
= n; 
NA Na 
under steady flow, or, alternatively, approximately a “dynamic intrinsic 
Viscosity” 


"8 
NA 


can be written for any given frequency. 

The curves in Fig. 17 are frequency dependent, however, although it 
turns out that ny, is only slightly so. Nevertheless, the considerable rise 
of ny above the pure solvent \ ISCOSILY, as the concentration is Increased, 


Indicates other mechanisms are being lumped into 94. As usual, some 














GRAMS (PER 100M 

Fig. 17-—Rigidity and viscosities of polyisobutylene (M, = 1.18 & 10%) in ey 
clohexane, at 25°C and 20 ke 
extensive distribution of relaxation times is probably responsible. How- 
ever, from the chemical point of view, it is best to see if some principal 
mechanisms related to known structures can be identified. If so, they 
could be associated with new ideas about the details of polymer intrinsic 
viscosities, as well as the form of isolated molecules. " 

First, the frequency dependence of the us, of the model of Fig. 16 is 
as shown on Fig. IS. Striking regions of dispersion appear, although 
more pots are needed to define the 10 evele Zone Actually, many sets 
of data have been obtained in the 10° evele zone. Recently, an immersed 
quartz tuning fork has given the approximate value shown for 2300 


eveles. The experiments of Fig. IS were on a polvisobutvlene having 


M, = 3.9 X 10°, dissolved in evclohexane. Values of ny and ny were, 


of course, also obtained. The results were then analyze d fora svstem ol 
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three Maxwell elements in parallel with again, as in Fig. 16, a solvent 
Viscosity, this time called m, (truly absolute solvent viscosity) in parallel 
with them. The us, curve of Fig. 18, running through the observed 
points could then be calculated, by some special trial methods, with 
Which Messrs. H. ‘IT. O'Neil and O. J. Zobel of Bell Telephone Labora 
tories kindly helped 

This analysis, identifving three principal relaxation regions for the 
motions of polvisobutvlene chains in evclohexane, gave for a | per cent 
solution (taken as linear part of concentration curve and hence equiva 


lent to high dilution 


Principal rigidities S00 dynes em 


9 


Ls 3,190 dynes em 
M4 84,000 dynes em 
Principal viscosities 0.0082 poise (pure cyclohexane 

n 0.255 pose 

3 0.006 poise 

"4 0.004 Poise 
Principal relaxation frequencies f. 550 cycles 

fs 8.45 X 10' cycles 

fy = 3.52 X 10° cycles 

Tentatively, these mechanisms may be schematically described as on 

Fig. 19. Here, the polymer coil, subjected to shear waves in dilute solu 
tion, exhibits rigidities Me, bw and M4, all shown on different seales. be IS 
the configurational elasticity because of actual changes in root mean 
square separation of chain ends, as from R to R’. It is retarded by vis 
cous drag through the solvent, nm. , which is presumably the main source 
of characteristically high 7, of chain polymer solutions. The relaxation 


frequency for this mechanism is low—a few hundred eveles. It may 


come in significantly in work on more concentrated solutions at low 
4s 7 . . 
where chain entanglement is nevertheless the dom 


frequencies, 


inant factor. 

uw; is When segments of the same chain in the molecular coil temporarily 
entangle with each other. Striking evidence has recently been given by 
Fox and Flory” that because of mutual interference, the theoretical 
random flight configuration of a chain gives very much too small 
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molecular coil volume, V,. This suggests that thermal agitation tending 
toward a smaller J ‘ and excluded volume or repulsions loremng a large 
one, Will cause collisons or entanglements which might last long enough 
to give a van der Waals cross-bond as denoted by crosses on the Me 
sketch Phe actual forces in these would somewhat resemble those be- 
tween different molecules in the concentrated solutions of Ferry.” 
This mechanism has the reasonable (based on Ferry’s and others’ work 
relaxation frequeney of 845 * 10°. A small viscosity, 3 , may comprise 
frietion of slippage at the entanglement points, with both the polymer 
and associated solvent molecules 

In Fig. 19, wy is a relatively high stiffness presumed to be some average 
hindrance to rotation of one segment with respect to another. In the 
sketch, close-packed spheres representing methyl groups in polviso- 
butvlene are portraved. Their force fields overlap more in some places 
than others, in the meandering of the chain to form the molecular coil 
of course, some tail-to-tail structures may be important here; they have 
all been shown head-to-tail in the sketch). Thus, this total internal steric 
restraint on chain flexibility, with a relaxation frequency of 3.5 X 10°, 
contributes greatly to the large dispersion of rigidity in the megacycle 
range noted in Fig. 18. The related viscosity, 4, 1s again low. 

There is no doubt a considerable distribution of relaxation character- 


istics associated with each and all of these mechanisms. 


Physical Prope rhies Per Molecule 


Since the viscosities and rigidities in the dilute solutions indeed seem 
to be additive with the number of molecules present, values of these 
properties, for the hypothetical mechanisms, can be expressed per aver- 
age chain. Of course, the measured quantities are expressed as constants 
per ce of solution, but it may be useful to think of in terms of one aver- 
age chain in each ee. Then, the shear deformation of this chain could be 
denoted by a force constant. The associated viscosities remain, however, 


dependent on solvent surroundings. Thus, for the polyisobutvlene of 


M, 3.9 & 10°, in evelohexane solution, at 25°C the molecular quan- 


tities : 


17 & 107° dyne em [nm] = 1.6 & 107” poise 


"i 


6 x& 10 dyne em | 3 x 10 ' poise 


16 & 10°" dyne em [m4] 24 X 10°" poise 


In the section on polymer liquids, the high-frequeney modulus wu Was 
attributed to a nearest-neighbor glass or erystal-like interaction (since 


the actual values were indeed typical of the hardest organic solids). 
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However, in polyisobutylene (and to some degree in poly-a-methyl 
stvrene), it is especially difficult to distinguish inter-chain from tntra- 


chain crowding of methyl groups. Thus, while average center-to-center 


separation of methyls is ~4 A in adjacent chains,” it is <2.5 A within 


chains, in polyisobutvlene. This crowding is apparently strong; the ob- 
served AH,,,, is only 12.8 keal per mole instead of the 19.2 expected. 
The energy of steric hindrance thus amounts to almost half of the actual 
heat of polymerization. It is reasonable that a large part of the hardness 
of a mass of polyisobutylene chains, such as in the liquids, should there 
fore reflect the same mechanism as that for M4 (Fig. 19) in the dilute 
solutions. A rough check on this can be made. A polyisobutylene having 
considerably lower molecular weight than 3.9 X 10° and thus inter- 
mediate between the “liquid” and ‘solid’ ranges, had a Maxwell sheat 
modulus in the megacycle region (14 me) of uw = 5.3 X& 10%, at 25°C. 
The number of molecules cc, with individual [f;] given above, necessary 
to give the observed density of this polymer was multiplied by [fs], giving 
uw = 2.8 X 10° dynes, em’. Accordingly, about half of the observed high 
frequency rigidity of polyisobutylene, at 25°C, may be calculated from 
a ‘molecular constant” embodying intra-chain stiffness. 

Much more refined and detailed treatments are required to generalize 
these ‘‘molecular constants” which are after all, as shown below, de- 
pendent on using a thermodynamically “inert” solvent. However, much 
as structurally significant dipole moments can be derived from measure- 
ments in dilute solutions, it seems hopeful that macromolecular me- 
chanics can be so elucidated. Also additional structures, such as poly- 
propylene and polydimethyl siloxane compared to polyisobutylene, are 
currently being studied. 


Temperature Variation 


Some further behaviour at different temperatures and solubilities of 
separate chains in dilute solution may now be considered against this 
background of possible mechanisms. Practically, these studies will bear 
on processing and properties, lacquers, paints, and casting solutions of 
polymers, as well as on the other qualities outlined in the introduction 
Results may be conveniently discussed in terms of the modified Max 
well single element, with factors ny, nx, and ws (Fig. 16). Mostly, the 
kiloevele range, reflecting molecular coil changes, will be of interest. For 
comparison, it may be noted that at 20 ke, the polyisobutylene whose 
ux = 1061 dynes em’ in I per cent solution in evelohexane receives 889 
dynes em of this from gw , the retarded configurational mechanism; 169 
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Fig. 20> Temperature variation of rigidity of 1 per cent solution of polyiso 
butylene (M 1.18 & 10°) in cyclohexane 


dynes, em’ from the entangled segments stiffness, us ; and only 3 dynes 
em from the intra-chain stiffness, uy. Thus, chain configuration me- 
chanics, including the associated viscosities, can be well enough thought 
of in the following paragraphs, in terms of us , ng and 4. 

The exponential decrease of « with temperature familiar for polymer 
solids and liquids is much suppressed in the us vs 7 curve of Fig. 20. 
While the ws, internal rotation, mechanism for single polyisobutylene 
molecules probably has a considerable activation energy, that for the 
uw. , configurational, rigidity should be very small. Then, without re- 
tardation, the intrinsic chain modulus would rise with rising tempera- 
ture. These influences seem to combine to give the modest decline of us 
appearing in Fig. 20. If these rigidities are plotted against 1/7’, the 
temperature coefficient is 2.3 keal. This is much less than the familiar 
values for the stiffening of rubbery solids, and emphasizes that 7nter- 


chain action reigns then. 


Sol ent | ariation 


Iifeets of solvents of different (mostly positive) heats of mixing on 


state of polvisobutvlene molecules in solution have been nicely estab- 


lished by Fox and Flory.” Especially, this work has clarified principal 


factors in the intrinsic viscosity expyession 


Yr 
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Here, V. = effective volume per molecule (and hence as determined by 
chain configuration), = molecular weight, a represents change in 
linear extent of molecule because of mutual interference of segments” 
and y expresses the hydrodynamics interaction of solvent and molecular 
coil (including varving degrees of ‘straining through” the coil). ws 
Interpretation of the mechanical properties of chains in dilute solution, 
with reference to the rough concepts of Fig. 16, arouses particular in 
terest in the factor a’. For a high molecular weight polyisobutylene, 
intrinsic Viscosity theory * indicated that a the ratio for volume of 
actual coil divided by volume for ideal random flight coil was 3.81 in 
cyclohexane but only 1.42 in beazene, both at 30°C. This striking alter 
ation in equilibrium chain configuration, a variable which is not readily 
introduced into polymer liquids or solids, appears in the inherent vis 
cosity vs ¢ curves in cyclohexane, Fig. 21, and benzene, Fig. 22. The 
large difference in {yn} at 25°C, 6.00 in evelohexane vs ~1.5 in benzene, 
indeed emphasizes the different solvent powers. Likewise, the large in- 
crease of [n| with temperature in Fig. 22 accents the poor solvent qual 
ities of benzene.” Too, empirically, polymer molecules which are either 
tight coils or are actually chemically cross-linked to form microgel mole 


cules characteristically show positive slopes of inherent viscosity vs ¢ 


plots.’ Accordingly, all this evidence for large changes in the conforma 


tion of chain molecules in ‘good’’ vs “poor” solvents should show up 
in dynamics of dilute solutions. Also, technically, quite different physical 
properties are found for polymer-plasticizer compounds where compat 
ibility is high (good solvent) than where it is low (poor solvent). Here, 


6.50 
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0.250 0.500 0.750 1.000 1.250 
CONCENTRATION IN GRAMS PER 100 ML.OF SOLUTION 
Fig. 21—Inherent viscosity of polyisobutylene (M, = 3.87 X 10°) in evelo 
hexane at 25°C. 
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Fig. 22 Inherent viscosity of polyisobuty ene Vi, 3.87 & 10°) in benzene, 


il various te mper itures 


more flexible Compositions are often produced with low compatibility 
plasticizers-—indeed, sometimes with those on the verge of phase sepa- 
ration than with those with highly favorable heats of solution.” This 


would mean that the bad solvents would compress the chains so that 


they would be more easily strained than if they were in a “free chain” 


or even extended configuration. If single chain, visco-elastic stiffnesses 
are acting this way, the dynamic uv, would then actually decline as heat 


of mixing become more positive. 














0.625 0.750 1000 
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visobutvlenes in evclohexane and benzene at 25°C and 
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This seems to take place, as indicated by the lower compared to the 


> 


upper curve on Fig. 23. Here, the us of the usual modified Maxwell 


model, at 20 ke, is plotted against ¢ for the polyisobutylene of M, 


3.9 X 10°. Also, the middle curve shows ys for a polymer of about a 
third of this molecular weight; while there is a small reduction in us 
with M, in this range, it is much less than the reduction caused by 
tightening up the polymer coil. 

The us values per average molecule, |fs|, fall from IS x LO dvne 
cm in evclohexane to 7 X 10 °° in benzene. (Of course, | f,} for the inter 
mediate moleculaft weight polymer in cyclohexane is only 5 & 10 °° be 
cause sO Many more molecules are present in solution 

The temperature dependence Ol upg also becomes nearly Zero at least 
over the narrow range from 25 to 50°C, in benzene compared to evelo 
hexane. This seems to accord with the indications previously, from Fig 
20, for a lower molecular weight polymer, that different mechanisms are 
competing. These may be the configurational, with w « 7, and relaxa- 


HCreases 


tion, With w« varying in some complicated way with 7’. Thus [7] 
markedly with 7, and presumably denotes an expanding molecular coil 
tending toward the “normal” configuration in cyclohexane. At the same 
time, the relaxation processes with rising temperature tend to cause the 
decrease in uw, typical of the upper, solid, curves on Fig. 24. In engineer 
ing use, often times poorly compatible plasticizers give compounds 
Which stiffen more gradually with temperature than do ‘“‘solvent”’ plas 
ticized ones. 

For similar reasons, the dynamic molecular. coil viscosity, ny , ought 
to vary less with temperature in thermodynamically poor than in good 
solvents. This is indeed seen in Fig. 25. On the other hand, n, for the 
modified Maxwell element has been described as the solvent viscosity 
with segment hindrance and restricted rotation terms from the polyme1 
molecules tumped in with it. These latter terms are presumably littl 
affected by over-all configuration (gv. term; the gw; mechanism will be 
somewhat affected, but not the uw, , on Fig. 19). Thus, 7, should have 
comparable temperature dependence in both good and bad solvents, as 


seems to be indicated by Fig. 26. 


Microge l Molecule Solutions 


The statistical coil of linear polymer molecules may be replaced by a 
chemically fixed, cross-linked network in microgel molecules.” These 
may be made completely rigid, like Einstein spheres, or highly swell 


able. The latter are hybrids between rigid spheres and coiled chains. In 
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synthetic rubber, they confer unique flow properties, causing the excel- 
lent processibility of GR-S 60. However, dynamic tenacity, such as in 
flex crack growth, is degraded by their presence. Now presumably the 
excellent extrusion qualities of synthetic rubber composed of from 60 to 
SO per cent microgel molecules are because of their individual shear 
stiffness. Thus, if a wire coating, for instance, is extruded at high rates 
of shear, chain molecules are deformed, and store energy just as dis 
cussed in the earlier sections on liquids. After emerging from the extru- 
sion die, they relax, and cause the gross retraction, shrinkage and rough 
ness shown in the wire insulation of the upper photograph of Fig. 27. A 
polymer with about 70 per cent microgel molecules gives the smooth 
covering shown in the lower specimen of Fig. 27. Here, the shearing 
stresses of extrusion seem insufficient to distort the tiny networks of 
the microgel molecule; in any case, the covering does not roughen or 
relax. Similar effects have been found for microgel plastics. Neverthe- 
less, unlike gross or macro gelation, the whole melt can flow. 

On this basis, dilute solutions of microgel molecules ought to indicate 
high shear rigidity per molecule. The mechanism ys of Fig. 19, in which 
now the junction points are not temporary, but are primary valence 
cross-links, should be predominant. Fig. 28 shows, for a polybutadiene 


. . 67 ° ° 
microgel in cyclohexane, that ws has indeed risen, compared to equal 
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Fig 26 Temperature Variation of m4 tor ] per cent solution of poly isobut 
(My = 3.87 X 10°) in cyclohexane and benzene 
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Ietfeet of mierogel molecules in svnthetie rubber 
Vire insulation 


gel 


on smoothness ot 
Rough covering is from high-speed extrusion of GR-S 


weights of chain molecules. Further, ace ompanying the extremely high 


average molecular weight of the microgel (18.6 & 10°), the (fx| per 
average molecule is 42 & 10 


dyne em or about twenty-five times that 
4 the polyisobutvlene with M, 3.9 & 10°. Also, the temperature 
‘oefficient for a, of polybutadiene microgel is low. 

Of course, polybutadiene, as chams or as microgel molecule segments, 
has many double bonds. These will surely influence the yy, or internal 
rotation mechanism 


Further work remains to show just what is their 


elect 


in the microgel case. But, it is interesting to Compare ys, values 
for Hevea rubber chains with those for, say, polvisobutvlene, which has 
only single bonds in the chain. 
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Hevea Rubber Solutions 


The comparison of equal weight concentrations of natural rubber in 


cyclohexane with polyisobutvlene in cyclohexane is surprising: 


Hevea rubber , Ws , 1350 dynes cm’, 1 per cent 
solution (corr. ) 


Polyvisobut vlene M, 2 * ue = 1000 dynes em’, | per cent 
solution (corr, 


Both results are at 20 ke. The higher value for natural rubber may be 
because of the double bonds causing stiffening of the chain. On the other 
hand, maybe easy rotation around single bonds raises the ws part. Cer 
tainly the viscous retardation within natural rubber chains is very low, 
as noted in the section on solids. However, its interaction with, or con 
figuration, in evclohexane may be peculiar. The {fs} per average mole 


° ° oh j con " 
cule is, however, low, being 15 & 10 ~ dyne em at 25°C. 


Polystyrene Solutions 


Much work, on light seattering and othe properties, has indicated 
appreciable intra-chain stiffness for polystyrene, ‘but still much freedom 
However, this work, as well as A//,,,,, of 
17 keal compared to ~19 keal calculated for no steric hindrance, sug 


6y 
t 


compared to polyisobutyvlene. 


gests comparatively small restraints on ideal flexibility. This needs to be 
checked by a frequency analysis of dilute solution mechanics, but poly 
stvrene seems to be a reasonable example ot plastic” behaviour at room 
temperature because of interaction befween the chains. (It is recalled 
that, earlier, a-methyvl stvrene polymer was cited as plastic model show 
ing both intra- and inter-chain stiffness. Unlike in polystyrene, the 
intra-chain factor shows up in a A//,,, of 9-10 keal, a third less than 
that calculated if there were no steric hindrance.) Thus, no evidence of 
unusual stiffness appears in Fig. 29, when, indeed, the uw, values are 
considerably lower, for equal weight concentrations, than those for nat 
ural rubber. The highly milled rubber studied had M, very nearly that 
of M, = (0.234 * 10° of the polystvrene, so the [fx] per average poly 


stvrene chain, 4.5 & 10 ‘dyne em is less than a third that of the rub 


ber. No wonder that at high temperatures, where the phenyl group 
interaction between chains is much reduced, polyst vrene makes a good 
rubber. Also, in Fig. 29 are shown data for a polymer of M, = 1.2 « 10°, 


made in emulsion and having |n| 1.350 in benzene at 25°C 
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The polystyrene solutions discussed above were in benzene, a good 
solvent. Here, the situation is converse to that for polvisobutyvlene ; lor 
polystyrene, evelohexane is a poor solvent and benzene, good. Hence, if 
the previous interpretation of reduced single chain quasi-configurational 
a stiffness is general for solvents of more endothermic mixing, the 
‘plastic’ molecule polystyrene should show it in cyclohexane. This is 
indeed evident in Fig. 30, showing one of the same polystyrenes of Fig. 
29, measured at 20 ke (normalized to | per cent concentration). Also, on 
hig. 30 are shown the inherent viscosity (practically, the intrinsic vis- 


cosity, in this ease) and the absolute viscosity of the | per cent solution 





[ & Mn = 
Mn 











30 40 50 60 70 80 
FREQUENCY IN KILOCYCLES PER SECOND 
inge of shear stiffness, ug, with frequency, for 1 per cent solutions 
hh bye HZeTVWu 


under steady flow, »,. These are all plotted against temperature down 
to phase separation, at about 26° to 27°C. 

The marked positive slope of the fnn,-/e curve denotes the large con- 
traction in molecular coil volume preceding phase separation or insolu- 
bilitv. Phe absolute viscosity, n,, however, rises with declining tem- 
perature because it is dominated by solvent VISCOSILY, but when the 
polymer phase comes out, 9, abruptly falls off. 

Phe uw,» values are consistent with this steady flow behaviour, except 
that the rzse of ip it the turbidity point seems to be because a laver of 
swollen polymer-rich phase forms on the torsional erystal surface. This 
condition is seen in Fig. 30 to coincide nicely with the abrupt changes 
in steady flow viscosity 


Phe slight maximum in the we curve at about 35°C may not be real. 
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It does come near the point of minimum interaction for the whole sys 
tem. In any case, as discussed before, the average temperature coefficient 
Ol uy in the Pool solvent is vers low compared to the good solvent. The 
Values of wy are roughly 3 to } those in benzene 


GENERAL THEORY OF SINGLE CHAIN MECHANICS; KUHN AND KIRKWOOD 


As noted before, much of the present understanding Ol stress-strain 
properties of polymer chains, in dilute solutions, liquids or solids, has 
come from W. Kuhn’s long interest in this subject. Many supplementary 
contributions have been stimulated by Kuhn’s work, and new points of 
view have been introduced by others. For instance, recently new and 
different proposals have been made about the flow birefringence and 
non-Newtonian Viscosity of solutions of deformable spheres " These 
ideas could be tested on suitable microgel solutions. 

Recently, moreover, an especially significant general theory ot Visco 
elastic behaviour of polymer in solution has been constructed by Kirk 
wood.” It explicitly considers the hydrodynamic conditions leading to 
the rigidity now observed for high-frequency shear waves. It formulates 
definitely the configurational changes of isolated chains in solution when 
strained in shear. As this theory is advanced to forms where simpler 
calculations can be made, it may answer many of the questions raised 
by the new experiments on single chain properties 
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CONCLUSIONS 


To leave some impression of the elemental chemical structures which 
move around when wood, rubber, plastics, textiles and finishes are used 
mechanically—that has been the aim of this study. Polymer VIscosities 
have been found in a variety of ‘solids’; rigidities have been demon- 
strated for very fluid “liquids” and solutions. Studies of these solid and 
liquid extremes have given some chemical reality to the classical spring 
and dashpot models. 

Kxistence of compressional viscosity. has been shown for polymer 
liquids and solids. It may comprise a new quality for investigation of 
polymet structure. At present, too little is known of its origin to inter- 
pret further the effeets of intense ultrasonic irradiation of polymer solu- 
tions. Experiments of Schmid and co-workers’ early indicated degrada- 
tion of molecular weight of polystyrene, so irradiated, but whether this 
is chemical, from local heating in the solvent, or actual physical coupling 
with the wave field, is still unsettled. However, these workers also con- 
sidered a compressional stiffness of the polymer molecules in the solu- 
tions,” and showed that if there was coupling, it Was not inertial (by 
dissolving polystvrene in solvents of exactly the same density, no redue- 
tion in effect was observed \ point of general interest arises here; 
impact fractures of plastics presumably actually fracture some primary 
valence bonds. This is certainly true for many thermoset materials, and 
probably for chain compounds. Hence, if the detailed mechanism of how 
compressional Waves move and perhaps rupture polymer segments were 
known, information on the baffling problems of ultimate strength would 
be gained. The observations above on dependence of \ and X’ on molec- 
ular weight and structure provide only the barest start on this but a 
new goal is in view. Too, basic questions of how rapidly molecules being 
formed in a polymerization equilibrate in temperature with thei sur- 
roundings are elucidated by compressional wave propagation Constants. 
For instance, absolute rate measurements on velocity of chain growth 
cannot be said to be isothermal if they seem to be faster than the ther- 
mal relaxation times which the ultrasonic measurements indicate can be 
~10 to 10 Sec, 

Likewise, more thorough understanding of velocity and dispersion of 


compressional waves in polymer solutions would clear up anomalies in 


velocity measurements for a wide variety of polymers,” some of which 


have been tentatively attributed to chain branching. 
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The Reliability of Telephone Traffic Load 


Measurements by Switch Counts 
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a . : ‘ 

The switch count method of tele phone trafic measurement ts subject to sam- 

pling errors, The nature of these errors is discussed and formulas are dk 
7 } 

rived which de scribe the fl lent of the C1TOrs unde 7 normal / ¢ neounte red 


traffic conditions. 


INTRODUCTION 


Of prime importance to the telephone traffic engineer is the deter 
mination of the busy season busy hour load carried by groups of trunks 
or other circuits of a telephone switching system. Three direct methods 
of measuring such loads are found in the field today These are: 


a. Peg Count and Holding Time Method 


The number of calls carried by the circuit group during the observa- 
tion period is counted. This number multiplied by the average holding 
time per call (in hundreds of seconds) and divided by the length of the 
observation period (in hours) gives an estimate of the group load in 
units of hundred-call-seconds per hour (CCS). The major drawback to 
this peg count method is that it requires a separate determination of the 
average holding time per call for the group under observation. R. I. 
Wilkinson! has analyzed the sources of errors of holding time measure- 
ments. In addition, correlation between load and holding time introduces 
an error which has not been studied. 


b. Switch Count Method 


At fixed intervals the circuit group is scanned and the number of busy 
circuits is counted. The total number of busy conditions counted divided 
by the number of scans is, then, an estimate of the load on the group in 
units of average simultaneous calls or erlangs*. This estimate is generally 
converted to CCS (1 erlang 36 CCS) by traffic engineers since the 

*The name “erlang” for average simultaneous call was adopted at a plenary 
meeting of the CCIF at Montreux in October, 1946 
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load entries of most traffic tables are in terms of CCS. For theoretical 
studies the erlang is a more convenient unit and will be retained here. 


c. Continuous Method 


The busy condition of each circuit is represented by a fixed increment 
of electrical current through an ampere-hour meter. The instantaneous 
current is then analogous to the calls simultaneously present so that the 
meter, Which integrates the current, may be calibrated to indicate hun- 
dred-call-seconds or erlang-hours directly. Although this method is po- 
tentially the most accurate, practical difficulties have limited its use, 


In addition to these direct methods, there are several methods of in- 
direct load measurement which, relying more heavily on traffic theory, 
make use of partial load indications, such as duration of group busy or 
the number of calls finding the group busy. Such measurements are less 
reliable than the direct measurements particularly when applied to un- 
derloaded groups. 

This paper is concerned with the reliability of switch count load meas- 
urements since this method appears to have prospeets of considerably 
wider adoption in the future. Main emphasis will be placed, both quali- 
tatively and by the application of error formulas, on the relative effects 
of various measurement and traffic parameters on the accuracy of switch 
count measurements. Where long derivations of formulas are required 
they are deferred to the Appendix. 


SOURCKES OF ERROR 


\s has been described, switch count measurements yield the average 


number of calls found present when a group of circuits is scanned at 


fixed intervals during an observation period. Usually only that period 


of the day during which the load is greatest is of interest to the traffic 
engineer. Because the load during such periods also fluctuates from day 
to day, measurements of the loads for several days must be averaged to 
provide a useful load estimate. 

here are two main sources of error, therefore, in switch count esti- 
mates of telephone traffic loads: 

lL. Hach individual count of busy circuits is separated from the next 
by a time interval during which changes in load are not detected. Con- 
sequently, the load indicated by measurement may differ appreciably 
from the actual load carried. This difference can be decreased by de- 


creasing the interval between scans. 





RELIABILITY OF TRAFFIC MEASUREMENTS 


2. Even if the load carried during a measurement period were known 
very accurately, it is still only a sample of the many loads that might be 
offered by the same source of traffic under statistically identical con 
ditions. Therefore, the average of several load readings may be expected 
to be somewhat in error as an estimate of the true average of the traffic 
source. The latter will be referred to as the source load to distinguish it 
from the carried load. 

Mechanical and human errors are likely to be present as well but, since 
they are not inherent in the switch count method, they will be neglected 
here. 


SWITCH COUNT ERROR 


As shown in the Appendix, for periods of observation which are rela- 
tively long with respect to average holding time made on trafhie with 
certain assumed characteristics, the average error of switch counts in 
estimating traffic load carried in the same period is zero. The coefficient 
of variation of the error, which is the standard deviation of the error 


expressed in per cent of the traffic load carried, is given by: 


100 4 | retnh(5) 2 | 


r 1+ ¢ : i 
etnh hyperbolic cotangent of 
2 | ” 2 


re= Tt > 20 
where r ratio of scan interval to holding time 
average holding time 
carried Joad in erlangs 
number of switch counts 
length of observation period 
N number of observation periods 


and where the following assumptions are made: 

a. Calls originate individually and collectively at random.t 

b. Holding times are exponentially distributed. 

¢. Congestion loss from the group is negligible. 

* T have recently learned that these carried load formulas have been published 
by Conny Palm in Tekniska Medelanden fran Kungl. Telegrafstyrelsen, AL. nr. 7-9 


tSee T. C. Frv. Probability and Its Engineering Uses. D. van Nostrand Co 
Inc., New York, p. 216, for a definition of this condition 
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As shown in the Appendix, this formula simplifies, when r 


WO / | 
v.32 

c/T 4 ba NTt 
where ¢ 7 rate of scan in eveles per time unit. From equation (2) it 
is apparent that if the sean interval is of the order of a holding time, the 
error of an estimate of trafhie carried is inversely proportional to the rate 
of scan and inversely proportional to the square root of average load, 
holding time and hours of observation. For example, take the case where 
switch counts are made during the busy hour, five minutes apart on a 
trunk group carrying calls with an average holding time of 3 minutes 
ISO CCS). What is the error in the 


nnd an average load of 5 erlangs 
estimated lond carried i the readings for ten davs are averaged ? (As 


sume conditions (a), (b) and (¢) are met.) We have 
\ 10 observation periods 
1 how 
1 20 hour 
5 erlangs 
12 scans pe observation period 
20) average holding times per observation period 


rom equation since Tt 20 and? T ‘ct ey, 
LOO 
12 VW 6-5-10-1-1 20 


If, as proposed in the Appendix, it is assumed that the error has a 
normal distribution, there is 90 per cent assurance that observed values 


will fall within 1L.6O4V,, or inthe example within 3.52 per cent, of the true 
e*. Note that this error limit would be halved if the rate of sean 


avera 
oubled or if four times as many hours of observation were taken. 


were d 
The coefficient of variation of the switch count error for constant values 


is plotted on Fig. | for one observation period 
For loads other than one erlang the coefficient of 
by ya’N. Thus in the example we have, 


of 7) tas a tunetion of 
ol a one erlang load 
Variation is found by dividing 
Using thre dotted curve, 
Lo 
\ +10 

i sufficient number of observations are taken so that a 

neglected in making an estimate of the universe 
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or more accurately using the solid curve, 


16 


VY 910 


The error of USINY equation (2) is seen to be ne gligible for most purposes 
even when 7’ ¢ is less than 20. The probability of an observation o¢ 
curing Within a given number of standard deviations is widely published 


for the normal curve. A few values are given below: 


Fig. 2 isa plot tor 40 observations of measured load vs carried load 
Kach observation was made for a half hour period on a panel line finder 
group with switch counts made at the start and middle of the period 
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1—Accuracy of switch count estimate of load actually carried 
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é 1 16 18 20 22 24 26 28 
ARRIED LOAD IN ERLANGS ESTIMATED BY 60 SWITCH COUNTS 
hig. 2. Accuracy of switch count estimate of true average load. 


Phis is compared with the average of switch counts made every 30 sec- 
onds which has a relatively negligible error. The average holding time 
per call for the group was 176 seconds. The accuracy of only two 
counts is surprisingly good and the observations are seen to le satisfae- 


torily het 


ween the 2o limits. 


ERROR OF TRAFFIC IN A GIVEN 


PERIOD AS AN ESTIMATE OF THI 
\D 


SOURCE 


Phe average traffic carried in two different periods but generated by 


the same traffic source is subject to statistical variation. As a result, any 


measurement of load, even if measurement errors are eliminated, is only 


i sumple of the wide range of traffic loads that might have been gen- 
erated by the same source of traffie under identical circumstances. 
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J. Riordan has shown? that the standard deviation of the average 
traffic load for any one period is given by 


Pat (| , 
rr Mt 


where a the average source load 
average holding time per call 


length of observation period 


(Assumptions are as before with an additional one that all periods are 


in statistical equilibrium ) 
When ¢ 7 <1 this reduces to the form also given by F. W. Rabe 


. Yat 
| T 
or expressed in a per cent of the average 
2t 


100 { r 
a 


When N periods of length 7 are observed the coefficient of variation is 
reduced further to: 


100 4 np 


In the example of the previous section, 


COMBINATION OF ERRORS 


Evidently if switch count readings ure used to estimate the ave rage 
which may be expected in other periods, the two errors described above 
should both be taken into account. The errors are probably correlated 
but this correlation is weak and at present no method of allowing for it 
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is evident. Such a refinement would probably change the equation for 
stundard deviation only slightly from that derived for the independent 
case; therefore independence will be assumed. The standard deviation 
of the sum of two independent variables is the square root of the sum of 


thre Squares ol the component standard deviations: 


ag \ OC; 


{ } ; ta’ Yta 
| ret ) ; —— i siecaes 
| 2 INT" NG 
a’ 


\ssuming Is approximately Unity, that is, that carried load is aupproxi- 
a 


mately equal to source load, 


In the example VIVeH, 
vs LOOT; 


‘There 1s, then, a) per cent assurance that the source average is within 
L.ob & 4.96 8.1 per cent of the observed average. Note that doubling 
the switch count rate (which halves the switch count error) reduces the 
total error only to 7.6 per cent (about 6.7 per cent improvement), while 
doubling the number of hours of observations reduces the error to 5.9 
per cent (about 30 per cent improvement). Plots of the coefficient of 
Variation of a one hour observation of a one erlang load versus scan rate 
for various average holding times are given in Fig. 3 for a wide range of 
holding times. The coeflicient of variation of error in estimating other 
loads may be found from Fig. 3 by dividing the unit load coefficient 
by VYaNT. In the example, the unit load coefficient is found, by entering 
Fig. 3 with ¢ 3 minutes and rate of sean c/T 12, 1 sean cycles per 
hour, to be 35.0 per cent. Dividing by y/5-1-10 gives a coefficient of 
Variation of 4.96 per cent as before. It is evident from Fig. 3 that in- 
creasing scan rates is not a universal Way to improve the aceuracy of 


SOUTCE load estimates. 


(HOICK OF SCAN RATES 


What then governs the choice of sean rate? Evidently increasing the 


rate mereases the aceuracy ol carried load estimates to any point de- 
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sired. This is far from true if source load is being estimated. If the cost 
of making a sean is constant, increasing the number of observation 
periods and decreasing the sean rate will improve accuracy of source 
load estimates without changing measurement costs. The number of 
hours available for measuring, of course, limits this procedure, while 
the increase in accuracy becomes negligible as r becomes large. On the 
other hand, if the cost of each observation is only slightly affected by the 
cost of making additional seans, a high scan rate might be justified. 

[ni applying the above relationships to traffic measurements, the 
usual question raised by the traffic engineer will be either how many 
hours of data need he take to be reasonably sure of his estimate or, 
conversely, how sure is he of an estimate based on available data 
Assuming as before that the error distribution is normal, the per cent 
plus or minus error limits within which a proportion, 2, , of the estimates 
will fall is given by zV, ; the value of 2 corresponding to any selected 
P, may be found from tables of the normal probability distribution. 
“Reasonably sure” is often taken to mean that there is 90 per cent 
assurance that the error does not exceed 5 per cent. When P, is 0.90, 
z is 1.64, so that under this condition 1.64V, 0.05, or V, 0.0305 
Given sean rate and holding time, V, is proportional to 1/./aNT' accord 


ing to equation (9) or Figure 3. When V, is held constant, aNT is con- 


stant so that the plot of log V7 against log a is linear, as shown in 


Figs. 4 and 5. The number of hours needed to meet any chosen reliability 





A» 


switch counts. for usage measurement 
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hig. 4--Hours of measurement required for 90 per cent assurance that error in 
estimating source load does not exceed plus or minus 5 per cent when measuring 
interval is 120 seconds 


requirements may then be read directly from such graphs. In the second 
type of question, z, V7, scan rate and holding time are fixed so that 
-V, is proportional to 1) ya. Plotting log zV, against log Wa again 
gives a linear plot as shown on Fig. 6. 

In the numerical example above, the limits of error corresponding to 
90 per cent assurance may be read from Fig. 6 which is plotted for the 
appropriate assurance, average holding time and scan interval. Reading 
the error limits at the point where the 10 hours measured line crosses 
ISO CCS (5 erlangs) gives + 8.1 per cent as before. Fig. 5 may be 
entered to find the total number of hours required to reduce this error to 
5 per cent. Reading at the point where the 180 second holding time line 
crosses ISO CCS gives 26 hours. 


QUALITATIVE EXTENSION OF THERORETICAL APPROACH 


The original traffic assumptions made in deriving the theoretical re- 
sults above are: 


a. Calls originate collectively and individually at random. 





RELIABILITY OF TRAFFIC MEASI MENTS 


LOAD IN ERLANGS 
0.2 0.40.6 1.0 2 4 6 810 
he i 





. + 


+ + 


Mw 
° 
< 





IN HOURS 


@ 
° 





URING INTERVAL 
300 SECONDS 


REQUIRED OBSERVATION 

















LOAD IN HUNODRED-CALL-SECONDS PER HOUR 





4 6810 20 40 60 100 200 400 1000 ‘4000 10,000 


Fig. 5—Hours of measurement required for 90 per cent assurance that error in 
estimating source load does not exceed plus or minus 5 per cent when measuring 


interval is 300 seconds 
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b. Holding times are exponentially distributed. 

c. Congestion loss from the group is negligible. 

d. Observation periods are in statistical equilibrium. 

How do departures from these assumptions affect the rehability of 


usage measurements 7 
a Holding Time Distribution 


experience i application ol delay and loss formulas has shown that 
theories based on exponential holding times are often applicable to other 
holding time distribution cases which have a wide range. However, for a 
constant holding time distribution special theories often are called for. 
The average and standard deviation of switch count estimates of carried 
load when holding time is constant, are given in part 2 of the Appendix. 
It is shown there that for estimates of carried load, 


Ve 
1 minimum V, 


(hE) 


, r t 
maximum } 100 —-A/ —, (r = 4F, 4, #, ete.) 
. 2V aN ae 


Since constant holding times found in practice are often very short, the 
case of r > 1 is the most likely to be met. For all values of r greater than 
one, the error given by formula (1) for exponential holding times is some- 
what greater than the error given by formula (10) for constant holding 
times, so use of formula (1) for the constant holding time case is con- 
servative. For values of r less than 1, the error is an oscillating function 
of r. The coefficient of variation varies from zero to 23 per cent above 
that for exponential holding times. Where r may not be accurately known 
the formula for exponential holding times again seems appropriate. 

ln making estimates of the source load when the holding time is con- 
stant, if r > 1, each scan is uncorrelated with any other, since no call 
ean be counted twice, and may be considered a random sample of traffie. 
There area total of Ne scans which have anaverage sean of a andstandard 


deviation V a. The average error in estimating a is, therefore: 


0 


with coefficient of variation 


\, = 100 4 + 100 { “opt 
ane ( 
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Equation (12) may also be derived with the procedure used for equation 


r 


(9) using o, = o; o,. For values of r large enough to make etnh( 4 


= 1 equation (12) is approached by equation (9). For smaller values 
of r (but with r still greater than 1), V, for constant holding times is 
less than V, for exponential holding times. When r 1, there is no 
carried load error. For values of r less than 1, the coefficient of varia 


tion of error in estimating source load average will vary from 


t t (1 r ) 
a 0 Ay 7 t 
NV aNT V aNt 2 
depending on the exact value of r. It is interesting to note that V, for? 


= (0.5 is the same as for r ZS. 


b. Loss 


The effeet of loss in the group depends upon the disposition of the lost 
calls. In general, accureey in measuring carried load increases with in 
creased loss because under these circumstances fewer load changes occul 
between scans. This is evident in the extreme case of a group which is 
100 per cent loaded; a single switch count gives a correct reading for any 
length period. Obviously load readings at 100 per cent occupancy are 
not very useful in estimating offered loads since the amount of lost load 
cannot even be guessed at. However, in the cases of lost calls held 
(Poisson) or cleared (Erlang B), the offered load may be estimated from 
the carried load (less and less accurately as occupancy increases) and in 
the case of lost calls delayed the offered and carried loads are likely to be 
the same even at high occupancies. With high loss, therefore, estimates 
of source load are subject to errors not considered in deri\ ing equation 
(99): however, switch count error in estimating carried load will be 


materially less than predicted by equation (1) 


ce. Random Call Origination 


On trunk groups which are alternate routes, calls may no longer be 
considered as originating at random. The resultant grouping of call orig 
inations will tend to decrease the RECUTACY of switch count measurements 
in estimating carried load; however, there is a corresponding decrease in 
accuracy in estimating the source load from the carried load so that ac- 


curacy in estimating carried load may be less worthwhile. 
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d. Statistical Equilibrium 


Statistical equilibrium may be thought of as the absence of trends in 
subseriber calling rates or holding times with the passage of time. The 
effect of trends on switch count accuracy In measuring carried load is 
very small except where the changes in traffic level are frequent and 
abrupt with respect to the sean frequency. Such traffic behavior is rare. 

Trends within the busy hour complicate the problem of estimating the 
average source load. However, it can be shown that if the trends are small 
in the order of 10 per cent to 20 per cent) little error is introduced by 
assuming that no trend exists. Large trends (in the order of 100° per 
cent), however, may indicate that the traffic source is so unstable that 
more hours of traffic data should be taken in order to insure that the 
sumple Is representative. 

Trends from day to day do not affect the source load estimates in the 


sume Way as within hour trends. The source loads are seldom exactly the 


sume on any two days although in most offices a load pattern is repeated 


from week to week. The traffic engineer may be interested in the average 
ource load of either a typical week day in the busy season or, some- 
times, of the average of the two highest days in the week. As long as the 
source load of each particular day remains close to the average for that 
day of the week, the general average for several different days of the 
week, will be known with about the same accuracy as if they had all 
come from a common source. If, however, there is no stable pattern in 
the source load, a third error in estimating the average is generated. 
There is some difficulty in determining whether or not variations in load, 
as indicated by measurements, are due to sampling variations or to an 
unstable source. Quality control methods might be used to detect in- 
stability but gathering and processing sufficient data for such an analysis 
might prove uneconomical. In general, if a traffie engineer feels that 
his source load is unstable he will need more hours of data than indicated 
by formula (9) to meet a given criterion of reliability. 


CONCLUSIONS 


\ theoretical approach to the problem of the accuracy of switch count 
measurements in estimating carried load and average source load has 
been explored. [It is believed that the assumptions made are satisfied 
sufficiently often in practice to enable fairly wide application of the re- 
sults of this exploration to traffie measurements. However, it should be 
kept in mind that where the assumptions are clearly not valid, special 
allowances will need to be made. In any case, the confidence placed in 





RELIABILITY OF TRAFFIC MEASUREMENTS 


usage measurements by a trafhie engineer is a function of his experience 
and judgment. It is hoped that the results of this study will add to the 
knowledge essential to sound traffie engineering. 


APPENDIX 
DERIVATION OF SWITCH COUNT ERROR IN ESTIMATING CARRIED LOAD 


WITH EXPONENTIAL HOLDING TIMES 


This derivation is based on a similar derivation by R. I. Wilkinson!, 
However, since load rather than holding time is of interest here, the em- 
phasis has been somewhat shifted. 

Assume that switch count measurements are being taken on traffic 
with: 

a. Calls originated individually and collectively at random 

b. Exponentially distributed holding times 

c¢. Negligible loss 
Let 7 = interval between seans 

average holding time 


traffic carried, in erlangs 


length of observation period 


= number of holding times in a sean interval 


number scans in observation period 
i 


yy 


_ number of holding times in observation period 


number of observation periods. 


Consider that the observation period begins with the first sean and 
ends 7 time units after the last sean. It is desired to find the error in esti 
mating the true load carried by averaging the number of circuits found 
busy on each scan. Following Wilkinson’s method we will first estimate 
the error of the switch count method in measuring the contribution of a 
single call to total usage and then modify it to take account of n calls. 
Calls of two types must be considered, those originating outside the in- 
terval and extending into it, Type I, and those originating within the 
interval, Type II. Both types may be subdivided depending on whether 
or not they extend bevond the end of the observation period. These are 
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indicated in hig. 7 Only that part of a eall which falls within the obser- 
Vation period contributes to the usuge of that period. Kirst the error 


made by switch counts in mensuring this contribution will be derived. 


1 I) pe / 


Consider a call which is already in progress at the start of the obser- 
vation period. Its duration beyond that point, according to theory, will 
be exponentially distributed about an average of ¢. 

If this duration, ¢, is between O and 7, the call will be counted once (a 
meusured contribution of 7 erlang hours) and a positive error of xc l l 
will be made. The same error will be made if / 27 — x so that the call 

STOP 


~ 





¢ 


CGraphieal indication of the t tvpes of calls with their two sub- 


is counted twice and so forth. Summing all the ways of making an error 
xr, We have: 


-f(ci — 2 (1) 


where f(7 xr) is the probability of f 
‘de 


Calls lasting hevond co neither start nor end in the observation period 


so that their contribution is measured without error. For these: 
P.O P(t 


| here hore? 
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Letting 


The moment 


Type Il 


Calls of Type IT may have either positive or negative errors given by: 


Pe<o(x) dz — |fil—x) + fila 


where f, (72 probability that a Type I] call has length ne rand 


ends before the end of the observation period. 
da 


probability that a Type II call starts n r before 
the end of the observation period and ends after 
the end of the observation period. 


| nt&-z 
1 
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Mquat ion (7) becomes: 


and noting that 


(4) 
' 2) ly 
cl—e” 


The moment generating function of this pair of equations is the sum 


of then separate m.g.f.’s: 


Ps ye" dy + ‘ge (y) dy ev" 


Neglecting terms of order higher than a’, 


Vila 


t },’ rk) =  (retnt(S) 2) (re — rk + 20’ | 


Now the number of Type I calls present in an observation is a vari- 


(11) 


able with average “a” and a Poisson distribution. Similarly the num- 

ber of “Pype TT calls is a variable, independent of the number of Type I 

calls, with an average of ‘a or “are” and a Poisson distribution. Aec- 
t 
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cording to the laws governing the compounding of variables the moment 
generating function of the sum of n variables y, when n is also variable 
with generating function G(t), is GOMla)) where Mla) is the moment 
generating function of 1. 
The generating function of a Poisson variable with average “a” 
‘so that 
G( Mila dats 
(12) 
G(Mr(a)) é ih is 
These independent variables may be added by multiplying their mo- 
ment generating functions to give the m.g.f. of the total measurement 
error of the carried load 
Al lai) an. Wterritemrion terete (13) 
From (6), (11) and (13) the following parameters are found: 


0 (14) 


j = 


a, = are | (retnt(5) 


If, now, re is sufficiently large 


oy, = V are | retnh(5) es (15) 


It is more convenient to deal with the standard deviation expressed 


as per cent of the carried erlang load, the coefficient of variation. This 
is done by multiplying both sides of equation (15) by ¢ to convert the 
time dimension from holding times to hours, dividing by 7 to convert 
from erlang-hours to erlangs, dividing by a’ to convert to proportion of 
carried load, and multiplying by 100 to convert to per cent. Assuming 


[n()- 2) 


100 { rm 


When NV observations are made this reduces further to 


LOO | wr |mn(5) | 


a: | 
,, IS approximately —,: 
» ~ @ 


Now etnh(r) 
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and forr < 2 


‘Therefore 


r r 
retnh ( =2+— 
2 6 


\ 


100 cae int : (r< 2) (16 
V aNTb c Ti ba’ NT t sae 


The error in carried load may be considered as the sum of a large num- 
ber of independent errors. Its distribution may, therefore, be expected 


to approach the normal distribution. Comparison of the third and fourth 


moments of the normal distribution with those of the error distribution 
which may be obtained from equation (13)) show good agreement for 


values of a’ greater than 1. 
DERIVATION OF SWITCH COUNT ERROR IN ESTIMATING CARRIED LOAD WITH 
CONSTANT HOLDING TIMES 


Wilkinson has shown! that, for constant holding time, switch count 
error in measuring the holding time of one call has an average 


and standard deviation 


where 7 
af 
negative error 
positive error 
Divide the problem into two parts: 


1. Forr> 1] 
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Min Or 0 for 
Max. Oz 


for 


Expressing this error in terms of carried load and proceeding as 
Part I of the Appendix 


r t 
Max. V, 100 { ; 
2'¥ aNT 


100 ; | 
c/T \ ta’ NTt 


Equation (20) compares favorably with the exponential holding time 


coefficient of variation of error of 


LOO | 
c/T V ba NT t 
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Network Representation of Transcendental 
Impedance Functions 


BY M. kK. ZINN 


Manuseript received November 5, 1951) 


The pur pose of the paper vs to show that the admittance or impedance of 
certain continuous structures, such as, for example, a finite le ngth of lrans- 
mission line of any sort, or resonant cavity, can be re presente d exactly at 
all. frequencies by a network comprising lumps of constant resistance R, 
inductance L, conductance G and capacitance C. The network will contain 
an infinite number of branches, in general, although a finite number may 
he used if it ts desired to represent only certain modes. 

The procedure is based upon a proposition known to students of function 
theory as “Mittag-Lefiler’s theorem,” which amounts, roughly, to an ex- 
tension of rational functions to apply to transcendental functions of the 
type encounte red in the theory of continuous structures. 


Several illustrative examples of the network synthesis are given. 


GENERAL 


Students of network theory are familiar with the fact that the im- 
pedance at a pair of terminals in a linear network comprising a finite 
number of resistors, inductors and capacitors, connected in any manner, 
is a rational function of the frequency having, in general, the fractional 
form of one polynomial divided by another. They are also familiar with 
the partial fraction rule whereby the function can be broken up into a 
series of elementary fractions, each of which exhibits one of the poles of 
the original function. This form is sometimes useful in the problem of 
network synthesis, where the impedance function is given and the ob- 
ject is to find a network having this impedance. 

The purpose ol the present paper is to show how a similar procedure 
can be carried out for certain transcendental impedance functions per- 
taining to structures having distributed constants, such as, for example, 
a resonant cavity or a piece of transmission line. The method employs a 
well-known proposition of function theory, which is usually referred to 
as Mittag-Leffler’s theorem. This theorem provides a tool for breaking 
up a transcendental meromorphic function into an infinite series of 
simple fractions in much the same way as the partial fraction rule is used 
to break up a rational meromorphic function. The series representation 
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NETWORKS FOR IMPEDANCE FUNCTIONS 


provides a means of determining a network of resistors, inductors and 
capacitors that will have an impedance equal to the specitied transcen- 
dental impedance function. This process will be referred to as obtaining 
a “network representation” of the function. If the given function is the 
impedance of some continuous (i.e., non-lumped) electric structure, the 
result will be an equivalent network for the structure. For other pur 
poses, such as, possibly, analogue methods of computing, the given fune- 
tion may not arise from any electrical structure. In either case, the net- 
work representations to be derived are possible only if the funetion 
satisfies certain restrictions, Which are stated in the section immediately 
following. 

The discussion is confined to transcendental impedance functions be 
cause of the technological interest in the electromagnetic structures with 
which they are associated and because they have not received as much 
attention as rational functions in the literature dealing with network 
svnthesis. The problem with which this paper is concerned can then be 
stated as follows: given, a transcendental impedance function satisfving 
certain conditions: to determine a network comprising elements of 
constant resistance, Inductance and capacitance whose driving-point im- 
pedance function, at a pair of terminals, will equal the given function at 
all frequencies, real and complex (except at the poles). 

For illustration of the procedure, three examples are given. The first 
is the impedance of a short-circuited or open-circuited transmission line 
in which the distributed primary constants, R, L,G and C are assumed 
to be invariable with frequency. The second and third examples are the 
impedances of resonant cavities driven in two different modes. In these 
examples the variation of resistance with frequency, due to “skin-effeet ,”’ 
is taken into account. 


IMPEDANCE FUNCTIONS 


The funetions under discussion will be referred to as “impedance 
functions” with the understanding that the term is meant to include 
“admittance functions” as well. By reason of the duality principle that 
rulis through all eleetric circuit theory, “any general proposition devel- 
oped for one must apply to the other. The functional designation, /'(p), 
will be used to denote either an impedance or an admittance function 


When a distinction is necessary, the impedance will be designated by 


Z(p) and the admittance by ¥(p). The independent complex variable p 


is the generalized radian frequency. For sustained sinusoidal currents 
and voltages, p lw 2rif where f is the real frequency.) 


For the applications contemplated, F(p) is a transcendental mero- 
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morphic function, which term implies that the function is given by the 
ratio of two entire functions, one or both of which is transcendental, and 
that the singularities of the function are ordinary poles, except for the 
point at infinity, which is an essentially singular point. In order to realize 
the particular network developments to be given, it will be supposed 
that the function satisfies the further restrictions given below: 

(1) All the poles lie in the left half of the p-plane with none on the 
Imaginary axis, 

(2) F(p) F(p). (The superbar denotes the complex conjugate of 
the unbarred symbol.) 

(3) Real part |F(2w)| > O for all real values of a. 

These three conditions are necessary to insure that the function is the 
impedance of a possible linear, passive electric circuit structure. Inter- 
preted physically in terms of this possible equivalent structure, the first 
condition specifies that the structure shall be stable; that is, every natu- 
ral mode of oscillation dies away exponentially. The second condition 
specifies that the natural oscillations are real functions of time. The 
third condition specifies that if a sinusoidal current flows at the driving- 
point terminals of the equivalent structure, the average real power de- 
livered to it will be positive. Since these three conditions, or their equiva- 
lents, are frequently mentioned in discussions of network theory, it is 
assumed that they are understood without more detailed explanation. 

In addition to the above restrictions on the form of the impedance 
function, the following two conditions, while not necessary, will be im- 
posed to limit the scope of the discussion: 

(4) All the poles of F(p) are simple. 

(5) F(p) = O(1), exactly, as | p| — © everywhere except at the poles. 

Condition (4), while limiting the scope of the exposition required, does 
not restrict the application of the results in any important way, because 
most impedance functions for which a network representation may be 
required have only simple poles. 

Condition (5) implies that as p increases along any straight line drawn 
through the origin and not passing through any pole of F(p), the modu- 
lus of F(p) either approaches a limit or oscillates between finite limits. 
The physical implication of this condition is that the response of the 


network as a function of time to a suddenly applied cause begins with a 


discontinuity of the same degree as that of the cause. For example, the 
current response of the network to an applied step of voltage begins with 
a finite discontinuity. This behavior is a characteristic of continuous 
(non-lumped) electromagnetic structures, which furnish the principal 
application of the network developments to be described. 
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MITTAG-LEFFLER’S THEOREM? 
Let the poles of the given function F(p) be p;, pe, ps ++ , where 


O<|m| < | pel < | Ds 


and let the residues at the poles be A; , Az, A; «++ , respectively. Suppose 
that it is possible to draw a sequence of closed contours, (, , such that 
C,, encloses py, Po, ++ P, but no other poles and such that the minimum 
distance of C,, from the origin tends to infinity with n. Suppose also that 
F(p) satisfies conditions (2), (4) and (5) above. Then Mittag-Leffler’s 
theorem gives the following series development for F(p): 


F(p) = F(0) + Limit ( An 
Nv P 


NN? n= = Pn 


The notation here used employs the convention that 
Pn Dn and A 


since, by virtue of condition (2), the poles occur in conjugate complex 
pairs. The value, n = 0, then allows for a pole on the negative real 
axis. 

Given any suitable function, the procedure is to determine its value 
for p = O and the location of its poles. The residues are next determined 
by 


A, = Limit:(p — p,)F(p). 


P—Pn 


Then the Mittag-Leffler expansion can be written down at once. 


NETWORK REPRESENTATION 


In the series (1) the terms occur in pairs with conjugate complex poles 
and residues. The object is to obtain a network representation of each 
such pair of terms. If F(p) is taken as an admittance, the branches rep- 
resenting the pairs of terms will all be connected in parallel; if F(p) is 
taken as an impedance, they will all be connected in series. 

Methods for obtaining a network representation for a rational func- 
tion, such as the one comprising a pair of terms in the series (1), are well 
known. It is only necessary to describe certain procedures of particular 
application to the present problem. Brune’ has stated that the necessary 
and sufficient condition for a network representation of a rational fune- 
tion of p to be realizable is that it be a ‘positive real function,” that is, 
a function that is real for real values of p and whose real part is positive, 
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or zero, When the real part of p is positive, or zero. In view of conditions 
(1) and (2) above, only one test” need be applied to each pair of terms 
of the series (1): the sum of a pair of terms will be a positive real fune- 
tion if, and only if, the real part of their sum is greater than, or equal to, 
zero for all purely imaginary values of p. 
The general term pair for which a network representation is sought is 
A, y y 
p — P,(p) — P,(0) (2) 
p— Pp, 
evidently two cases can be distinguished at the outset, depending upon 
whether 7?,(0) is positive or negative. If P,(0) is positive, the network 
branch, in order to be realizable, should be designed to represent P,,(p). 
The left-over negative term, —/7,,(0), then can be absorbed in the posi- 
tive first term, F'(O), of the series (1); more will be said of this later. If, 
on the other hand, ?,,(0) is negative, the network branch should repre- 
sent the whole term, P,,(p) P,.(O). This procedure insures that the 
real part of the branch impedance will be positive, or zero, at zero and 
infinite frequencies. To guarantee that the resistance is positive at all 
other frequencies requires further tests now to be specified. 
Let the real and imaginary coefficients of the poles and residues of 


the n- term be 
1B n 
a, — 1b, 


(W ith this notation, a, ana 6, are always positive ; a, and b, can be either 


positive Ol negative.) Then (dropping the subscripts) 


(aa b3) + 2ap 


P(p : 
a’ + p- + Zap r Dp 


fm/2..\) 2( aa bB)(ae + Bp) + 2w 
R\P(w)} = 


T aw WwW” 


: oe F ; 3, 2 ‘ i 
2 (ala 3a bB 3aaB" + bB’)w” — 2(aa — bB)w 


R{P Ue P(O)} 

{a> + 37) [ (a? + 87)? + 2(a? — B*)w* + w*| 
The necessary and sufficient conditions” for the real part of a rational 
function of p to be positive, or zero, for purely imaginary values of p are 
that the funetion be positive for p—> +7: and have no imaginary roots 
of odd multiplicity. When this test is applied to the functions P(p) and 
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> 


P(p) P(O), as given by (3), the following conditions are obtained: 
P(p) will bea positive real function if, and only vf, 
aa b3 > 0: 1.e P(O) > O 
and 
aa + bp > 0 
P(p) P(O) will be a positive real function if, and only vf, 
la bB < 0: m i P(O) < O 
and 
aa 3a bp 3aaa" T hg < 0 


If all terms of the series satisfy one or the other of these conditions, 
network branches can be devised to represent all the terms and all the 
R, LL, G, C elements of the branches will be positive 

In case all the terms are of the type where P,,(0) is positive, so that 
the network branches are made to represent Fal P), the left-over constant 
terms can be collected and added to the first term, F(0), of the series. 
This collection of terms then must be represented bv a final branch of 


pure resistance, or conductance, of value, 


F(0) — > P,(0) 


If the sum of the variable terms approaches zero for p— +7, the 


final constant term supplies the high frequency resistance of the fune 
tion F(p) and since this must be positive, if condition (3) is satisfied, the 
final resistive element will be positive. If the series converges non-uni 
formly, the sum of the variable terms can have a value other than zero 
as p—> +7 in spite of the fact that every term approaches zero indi- 
vidually. In that case (see example 1) all or part of the high frequency 
resistance may be supplied by the sum of the variable terms 

In case all the terms are of the type where P,,(0) is negative, so that 
the network branches are made to represent the sum, P,,(p P.(0), 
of the variable and constant terms and the series is uniformly conver- 
gent, all the high frequency resistance is provided by the branches rep- 
resenting these terms. The first term, F(0) then supplies the de re- 
sistance, Which is positive by condition (3). Non-uniform convergence 
can modify this division of high- and low-frequency resistance, however. 

Cases can arise in which the series contains terms of both types. In 


such a case the de resistance, or high Irequency resistance, OF both, ol 
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the given function might be less than the sum of the variable terms for 
these frequencies, with the result that the final resistance branch would 
he negative for either the series or parallel type of network development. 

To make the procedure as concrete as possible, particular forms of 
networks are described in the section following with explicit formulas for 


computing their elements 


NETWORK FORMULAS 


Simple forms of network branches are shown in Figs. 1 and 2. Those 
of Fig. 1. referred to as branches of “the first kind” are suitable for con- 
nection in parallel where the given function F(p) is an admittance, ¥(p), 
while networks of ‘‘the second kind,” shown in Fig. 2, are suitable for 
connection in series to represent an impedance, F(p) = Z(p). The net- 
works of Figs. la and 2a apply where the value P,,(0) of the general term 
is positive, while Figs. 1b and 2b apply where ?,,(0) is negative. Figs. 3 


and 4 illustrate, respectively, networks of the types of Figs. la and 2a 


Ln 
-- _ YU -- 
Rn os — 


or Sn 
) (bd) 


Fig. 1 General branches of the first kind 





Fig. la Fig. 1b 


(use where F(p) Y(p) (use where F(p) = Y(p) 
and Y,(0) > 0) and Y,(0) < 0) 


, Bila, + 8.)'(a. + bd.) 


N 
M(ai, + B.) 


Y(Q) 
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connected to form the completed network with the final non-reactive 
branch, G, or R, , in place 

Formulas for the network elements are obtained by equating the poles 
and residues of the network impedance function to the given poles and 
residues of the general term of the series. Since both poles and residues 
occur in conjugate complex pairs, and since equality of real and imagi- 
nary parts is involved, there are four equations, which are necessary and 
sufficient to determine the four constants, R, L, G@, C, of the network. 
The formulas that are obtained by solving these equations are give! 
beneath Figs. 1 and 2. 

The values given for G, and FR, in each case assume that all the terms 
of the series are of the type specified for that case 


























(dD) 


2 General branches of the second kind 





Fig. 2a Fig. 2b 
(use where F(p) = Z(p) (use where F(p) = Z(p) 
and Z,(0) > 0) and Z,(0) < 0) 
Bilan + Ba) (an + b5) 
2M: 
M? 


9 9 


B2 (a, + b’) 


Anh, — baBn 
' M 
N 


Galin P ° 
Mai, +8;) 


R, = Z(0) — >> Z,(0) R, = Z(0) 


where V/ = an(B, - an) + ZanBrbn 


N = —a,a, 4. 3a. baBn a 3dnanB, - b,.3*, 
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In the case of the parallel-type networks (Figs. la and Ib), p, 

a, + 28, is a pole of the admittance, ¥(p), and A, a, + tb, is the 
corresponding residue. In the case of the series tvpe network, the same 
symbols represent a pole and residue of the impedance, Z(p). 

Phe networks specified by Figs. 2a and 2b are duals of the networks 
of Figs. la and Ib, respectively, and are obtained from the latter merely 
by replacing L, by C,, R, by G,, and vice versa. 

The formulas are intended to apply to complex poles. They can be 
applied to real poles by taking 6, and 8, equal to zero and doubling the 


residue, a, , but this procedure is unnecessary, because the network rep- 





T 


4 
eo? 








Network of the first kind (branches la 


Cy 























Ro Le 
Fig. 4--Network of the second kind (branches 2a 
resentation of a real pole can be found readily enough by Inspection of 
the impedance terms involved. (See Example 1.) 
Phe above discussion is intended to sketch a general picture of the 
procedure, Individual cases may involve considerable detail that can be 


understood more readily by reference to the next section. 
APPLICATIONS 


Krample fa: \ transmission line with its far terminals short-circuited 
affords a simple illustration of the equivalent network theory. Let it be 
assumed that the parameters, R, 1, Gand C of the line are constants. In 
the more advanced examples to follow, the variation of these parame- 
ters with frequency for a particular kind of line will be taken into con- 


sidet ation 
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The impedance of the short-circuited line (Fig. 5) is 
Z Zo tanh 1-0 


where Zo 18 the characteristic impedance and is the total propagation 


constant of the line. We have 


Z (= a4 
(, + pe 
IT \((R + pL ; pe 


R,L,G and C being given for the total length of line 
To obtain il deve lo yment in terms ol network brsre hue s ol thre kind 
} 


shown in Fig. 1, we consider the admittance function, 
) y, coth T 1-3 
where } 1 Zand } |Z, . Our first task is to find the poles of this 


funetion and the residues Since the complex frequency variable p occurs 


L,G,¢ 








Fig. 5 Short-circuited transmission line 


under square roots in both Zo and Ty it might be suspected, offhand, that 


the singularities of the function are branch points rather than poles 


Such is hot the case, however. There are ho branch points and all the 


poles are simple. 
The singularities ol im nre to he found among the Zeros ol tanh ¥ 


which occur at 
, LI, 
To determine them, we solve 
r 
and find these roots: 


pP 


where 
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For n 0, the above would give 


sR G 
' L’ ¢ 


sut if we let [ — 0, so that tanh T — TI, we find that only the point, 
RL, is a singularity of Y; the other point, —G/C, is a regular point. 
Therefore Y has only one real singularity. 
To find the nature of the singularities of Y, we next caleulate 


; Dp — Dn 
Limit = A, 
win E- tanh =| ; 


and find that at each p, the limit exists and has the value 


-. tb, ( 1-8) 


' 


. Ow, ae er 
where I'(p, ; I'(p), evaluated at p = p,. The fact that this limit 
ap 
exists shows that all the singularities are simple poles. The values of A, 
are then the residues at these poles. 
When we now apply formulas (6) to determine the elements in the 
general branch of the equivalent network of Fig. la, we obtain, for 


i] > 1. 


| rn Gn G 


‘S4 
Ligh LC’ C, bg (1-8) 


Phe network then comprises an infinite number of such branches in 
parallel. Each branch has the same elements FP, and L, , equal, respec- 
tively, to half the total resistance and inductance of the transmission 
line, but the elements G, and €,, decrease from one branch to the next in 
Inverse proportion to the squares of the integers. 

Phe Q of the n™ branch, which can be regarded as the Q of the asso- 


mated resonance of the short-cireuited line, 


(1-10) 
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Thus, for small dissipation, the resonances would become sharper in di- 


rect proportion to the frequency (if the parameters PR, L, G, C, were 
invariable with frequency, as assumed). 

The above described branches of the equivalent network account only 
for the complex poles (n > 1) of the admittance function. Two more 
branches remain to be calculated. One is for the real pole (n 0), which 

, , Lg : . 
occurs at po R/L, with residue, Ao 7: The required branch for 
this pole is 


Ag 
P Po 


The other is the final conductance branch, which is evaleulated as follows: 


1-12) 


G, = Y(0) 4 pa An : V - coth VGR 


n 20 P i 


: | 
-27¢ > 
nao wen? + GR 
so that, for this example, the conductance branch vanishes. The network 
is drawn in Fig. 6. 
A series type of network, as shown in Fig. 7, can be determined by 











c : oe 
Fig. 6—Network of the first kind equivalent to the short-circuited line of 
Fig. 5 
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Ll 


) 
/ \ / 
Fig. 7 Network of the second kind equivalent to the short-circuited line of 
Fig. 5. 
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similar means. Since, however, it is a dual of the parallel network of Fig. 
9 for the open-circuited line, next to be discussed, it can be drawn im- 
mediately, without further calculation, once the latter has been found. 

Lrample 1b: We now calculate a network for the same line with its 
far terminals open (Fig. 8). To obtain a network of the first kind, with 
branches in parallel, we deal with the admittance function, 


r yy) tanh T (1-14 ) 


Phe singularities of Yo are found among the zeros of coth P, which occur 
il 


-15) 


Phe points p R/L and —G)C are both regular points this time. 
G Cisa zero ot ¥.) The singularities are simple poles, as before, with 


residues, 


\ . 1-16) 
mg ~ 4 v7 .\e~ ) 
Zi (pr) V(p 
its betore. 
The network branches for the complex poles are therefore obtained 


merely by putting mw + ; inh place of the nm in all formulas of the short- 


eireuit network. There is no branch corresponding to the branch Rk + pL 


of the other network and the conductance branch is again found to be 
zero, The complete parallel network is drawn in Fig. 9 and the series net- 
work, in Fig. 10 

It will be observed that the series network of Fig. 10 is the dual of the 


6 


_— 


a 





y 4 





Fig. 8 Open-circuited transmission line. 











dequivalent to the open-eireuited line of Fig. 8 
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parallel network of big 6 and the series network of Fig. 7 is the dual of 
the parallel network of Fig. 9. These dual relationships are of course a 
result of the fact that the impedance of an open circuited line is the dual 
of the impedance of the same line when short-circuited 

Kxample 2: Short-cireuited Concentric Line (or Toroidal Cavity with 
I} Radial). Vhe preceding example considered a fictitious transmission 
line of invariable parameters, R, EB. G. €, having a perfect short circuit 
at one end. The present example has to do essentially with the same 
problem but considers it from a more practical point of view. The vari- 
ition of Rand L with frequency is taken into account and the impedance 
of the “short-circuit” is no longer neglected. 

Let the line be the piece of coaxial cable plugged at both ends with 
conducting material as illustrated in Fig. 11. Considered from an alter- 


native point of view, our line is now a toroidal cavity oscillating in the 


wf} 























477e 47re 


Network of the second kind equivalent to the open-circuited line of 


mode where the electric force Fis directed radially and the magnetic 
force H lies in planes at right angles to the axis. If we assume the cay ity 
to be excited, or “driven,” from one end,* the impedance that is effective 
in defining the selective characteristic of the cavity with respect to fre 
quency is the total impedance at that end, that is, the sum of the im 
pedance bs " viewed into the cay ity, and the impedance, Le ‘ of the ad 
jacent end plug. Therefore, we have to deal with the impedance, 


Z Z, + Zo. 2-1 


By “impedance” is here meant the same thing that one considers in look- 
ing at the problem from the point of view of transmission line theory, 
namely, the complex ratio, for exponential oscillations, of the voltage 
between the inside and outside evlindrical surfaces to the total current 

* For determining the ‘‘natural frequencies”’ of oscillation of the cavity, it is 
immaterial at what point along it the impedance is taken; the total impedance 
at every point has the same roots The impedance is, nevertheless, not the Sime 


at all points so that the behavior of the cavity, when driven, will depend to some 
extent on the driving pom\ 
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flowing axially in the inner conductor at the same point. The zeros of Z 


define the natural frequencies of oscillation of the cavity and their asso- 


ciated damping constants, or Q’s. Our task is to develop an equivalent 
network for this Z. 
We have 

















10 9 
3677’ 
), FOR COPPER IN M.K.S. UNITS 


FOR AIR IN M.K.S. UNITS 


FUNDAMENTAL MODE ) 


Toroidal CAVITY, E radial, 
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where 


R pl UI T,( oa) 
tT 4 


2ra l(ea 
G + pe 


y 
log 


a 


h, a,b cavity length, inner radius, outer radius, as shown in Fig. 11, 


(m)" 
n 
q 


o { pug) : 


all measured in meters 


wu, gare permeability, conductivity of the conducting material of the walls 
(for copper: Ir(10 ‘),g 9.8(10°) in M.KAS. units 

Ho, Jo, €o are permeability, conductis itv, dielectric constant of the dielec- 
tric material occupying the cavity (for air: yy br (10 *), go (), 


é) = (10°) 367 in M.K.S. units), p = generalized frequency vari 


able. 


I,(z), /\(z) are Bessel functions of the first kind for imaginary argument 
and of order 0, 1. 

Ko(z), Ay(z) are Bessel functions of the second kind for Imaginary argu 
ment and of order 0, 1. 

Except for ignored small deviations of the field around the corners of 
the cavity, the above formulas are exact. To arrive at results that are 
sufficiently compact to be useful, we make these approximations, at the 
start: 

| tien b 
lr a 


where 


1207 ohms 
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Krom th 


Having in mind microwave applications, where the moduli of the : 


ments of the Bessel functions are >3000, we take 
Iy(z K 
/ 


Ly b 
I log (2-12 
2r a 


R + pl. 


Also, we have in mind only air dielectric and assume any loss therein to 
be negligible: that is, we assume G Q) 


All further approximations that are made are either 
A + 2A)" =1+A 


where, for an an space enclosed by copper walls, and lor frequencies on 
the order of 80,000 megacyeles, A is on the order of 10 *. For cavities 
made of other materials, the results obtained may not be sufficiently ac- 
curate and the problem would have to be reviewed from the start. In 
particular, the results do not hold for a cavity having walls of magnetic 
material, because we assume here that the permeability of the metal walls 
is the same as that of Mr, 12... 4 wv 

Po obtain an equivalent network of the first kind, we deal with the ad- 


mittunee, which is, from (2-2 


where [1 1A 
Phe poles of Y are then the zeros of | pe, Which are obtained by 


successive approximations. We first make a close estimate of the zeros by 


assuming that the impedance of the short-circuiting plugs is zero; that 
‘ussume, Z 0, whence p |. To obtain this estimate, we have 


Pah log | 


a 





NETWORKS FOR IMPEDANCE FUNCTIONS 


3(10°) meters per second. The approximate solution is 


iaeiilh a (: ix) 


Por j and ’ (ponmg) 


where 


Next we improve our estimate of the zeros by the well-known method in 


volving the derivative of the function, | pe es with respect: top, 


evaluated at p,,. This now takes account of the actual impedance of the 


end-plugs, The values of the zeros, so obtained, are 


P, a p DP, , 1D, 


( l | ) 
" \2d6, hé 
Idd hé 


where 6,* is the real part OF do, That is, 


where 


) 


6 Wind é 


where 


hi 


As an incidental matter of interest, the above gives the Q of the cay ity 


at any resonance, hamely 


(2-16) 


For example, the dimensions, a oem., b 1.0 em., h oem, pro- 
vide a cavity that resonates at about 30,000 megacycles. Then the Q’s at 


the first three resonances would be as follows: 


Wiiy 
) 
) e 


or 
30,000 « 10) 1250 
60,000 « 10) 6OL0 
90,000 &K 10 7360 
* For any frequency, 6 oug/2)'* issometimes referred toas the ‘‘skin depth’”’ 


because it is the de pth of metal at which the current density falls to 1/e times its 
value at the surface of the metal 
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The importance of including the effect of the end-plugs in determining 
Qis shown by the fact that, if they were assumed to have zero impedance, 
(at the first resonance would be 12,120 instead of 4250. 


‘To determine the residues at the poles, we write 


I, (| p)(1 - oF " F'(p) 


2(1 err (r(p) 


and then the residue at a simple pole p,, is 
F'(p,) 
_? A (Z0e5) 
(; (p, ) 
This limit is found to exist, showing that the poles are, in fact, simple. 
The value found for the residue, A, , 1s 


A, a, + tb,, . A, n — Uy 


H Win ( | | | ) 
™m 2d5, 2h, 


FT yc ( l l ) 
wh 2d6, 2hé, 


When formulas (6) are applied to determine the elements of the tuned 
branches of the equivalent network of the first kind, the results are, for 
the n*" branch, 


Kk Th ( l l ) 
2wWo, 2d6é, 2hd,, 


2° 
Won | T 
( dé hé, 


Win 


2hé 


( 1. 8 ) 
"\ dd, | 2hé, 
In terms of the R, Land C of the piece of coaxial line, the elements of the 
AL i’ oa ) 
, : ~ 2d6, 2hé, 
a 
3... 8 ) 
26, 2h, 
3 3 
it ae) 
( 2d, t 2hé, 


nm" branch are as follows: 
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where 


h 


MY Nog 
-7T a 
2re 
™ h 
a 
The network is shown in Fig. 11. 

It will be found that a “leakage” element, G,, , appears In the equiva- 
lent network, although the air dielectric in the cavity Was assumed to 
have no leakage (G 0). This element arises from the end-plugs and is 
necessary to account for the dissipation in them 

To obtain a network exactly equivalent to the eavity at all Trequen 
cies, we should add a branch corresponding to n 0, as was done in 
example 1. This branch would make the equivalence hold down to and 
including zero frequency. But, inasmuch as the approximations that have 
been made hold only for the high trequencies, where the resonances oc 
cur, it would be inconsistent to add this branch. What has been arrived 
at, then, is a partial network representation that gives a close approxi- 


mation to the impedance of the cavity at high frequencies, only 


Example 3: Toroidal Cavity with FE Axial. For further illustration, we 


consider another mode of oscillation of the short-circuited concentric 
transmission line investigated in the previous example. This time it is 
assumed that the radial electric force vanishes while the axial electric 
force between the end-plugs exists. The magnetic force is directed in 
circles concentric with the eylindrical central conductor, as before. This 
situation is illustrated in Fig. 12, which is the same as Fig. 11, except for 
the new disposition of the /-vector. 

For the new mode of oscillation, where the wave is a eviindrical one 
propagated back and forth between the inner and outer conducting 
evlinders, the oscillatory space Is naturally thought of as a “toroidal 
cavity,”’ while, in the previous example, where the wave was propagated 
axially back and forth between the terminal dises, the space Was called 
a “concentric line.’ Actually, the cavity itself has the same geometric 
form in the two cases. A practical distinction may exist, however, in that 
the axial mode of oscillation could be more easily excited in a cavity 
whose axial length is large compared to its radius, while the cylindrical 
mode would arise more easily in a flat ‘“‘pillbox’’ cavity whose radius is 
large compared to its axial dimension. 

The approach to the problem will be that of transmission line theory, 


as before. This time, the “line” comprises two circular dises between 
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which the eylindrical wave is propagated. The series impedance and 
shunt admittance of such a line are functions of the radius and so will be 
designated Z(r) and ¥(r), respectively. Their values are given below 
») ‘ } 
r ~~ T l i 
Zr) n Wit 


2nr 
9» 
-Z) 


V(r) 
h 


These formulas take into account the losses in the flat walls but assume 


the conductance of the air between them to be negligible. Losses in the 
inner and outer “short-cireuiting” evlinders will be taken into account by 


the boundary conditions 





2(b-a)C 7, _ - ) 
mene \ 2hdé,  2(b-a)d, 
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12. Toroidal eavity, F axial. 
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If V is the voltage between the tlat faces of the cay Itvata radius r and 


/ the total current in the lower face at this radius, we have 


dV 


dr 1Z(r) 


dr 


By differentiating, 


) . IWE 
(in + lwuh) 


Which is a squared propagation Constant, inde pe ndent of r, and 
| dZ I 
Z dr r 
Therefore, 
dV iladv 
dr’ r dr 


0 


is the differential equation for the voltage. The usual solution of this 
equation is a linear combination of /)(yr) and AKo(yr) but since, in this 
case, the arguments will be almost purely imaginary, it is more con 
venient to employ the pair of functions, Jo(—?yr) and No( —7yr) 
The solution for the voltage between the upper and lower surfaces at 
radius ris 
Vo Ad tyr) + BNo( — tyr (3-6) 
and, from this, the total radial current in the lower surface, at that 
radius, is 
1 dV : P ’ 
I(r) y lr) ASi(—tyr) + BNy(—?tyr)] 


where 


The impedance at the inner radius a, looking outward, is the 


‘ V(a) : AJo(—tya) + BNo(—iya) 
Z\la 1Zy(a) : 
l(a AJ\(—iya) + BN,(—iya) 
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The total impedance at a (inward + outward) for which we require an 
equivalent network Is 
Z Z\(a T La 
where Z, is the impedance of the central plug to axial current, viz., 


nh I)(oa) 


Lo (3-9) 


2ra I,(ca) 


To evaluate the constants A and B, the following boundary conditions 


are imposed at radi a and b: 
ata: V Via), a given voltage 
at b: V [(b)Z, 


where Z, is the impedance of the other “short eireult,”’ comprising the 
outer cylindrical wall. It is given by 
Z, nh Ko(ob) (3-10) 
2rb Ky(ab) 
xcept for ignored small deviations of the field around the corners 
of the cavity, the above expressions are exact. The process of finding 
the singularities of Z by successive approximations results in expressions 
that are too long to write down here. To obtain results sufficiently com- 
pact for engineering use, we resort to the following asymptotic approxi- 
mations for the Bessel functions: 


J (2) 


sin (2 3/4) 


[o(z) Ko(2) 


: : ~ | 
I,(2) Ky (2) 


Also, with an error on the order of 10 ‘ 


hes 
Z(r) ~ —" Ki (r) = 1/Ho(r) 


2nr 
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These substitutions result in the following asymptotic formula for the 


total impedance Z at radius a 


i) » 
“7 n ‘ 

cos kx + 1 (1 + :) sin kr 
No ’ No 


a 
cos kx + ” sin ka 
No 


Z = Kya ) 


bh ; 
where k = - — landz = —v1ya. 
a 


To find an equivalent network of the first kind to represent Z, we deal 
with the admittance, Y 1/Z. It is instructive and saves much work 
to put Y in the form of exponential functions, with the substitution 
1 — 

n — 0 


which is the reflection coefficient at both inside and outside cylindrical 


p= 


surfaces of the cavity. By this means we obtain 


y « His) OO = (3-13) 


2(1 pe ed) 


This is now identical in form to the formula (2-13) of example 2, where 


the /’-vector was radially, instead of axially, directed. In fact, since 


and 


‘ y 
he 


comparison with the similar formulas of example 2 shows that all the 
results of that example can be made to apply to the present one merely 


by changing the dimensional parameters as follows: 


Example 2 Example 3 
(E radial) (i axial) 


h goes Into bh—a 


2ab log (b/a) 
d= goes into 
pee Ee * 


The first result of interest is the value of Q), which is 
] 


Qn = hi, —— Df 


“ 
14 b—a 
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where, as before, 6, is the “skin depth” equal to the real part of o, . That 


IS, 


{ 
To gain an idea of numerical magnitudes, consider the same cavity used 
in example 2. The dimensions are, as before, h oem, b — a Oem. 
Kor the square cross-section chosen, the first resonance again occurs at 
30,000 megacycles, very nearly, and we can make the following direct 
comparison of the Q’s for the two modes of oscillation : 


FE radial 
30.000 *« 10E 1250 1370 


60 000 * 10° 6010 6180 
90 000 *« 10° 7360 7560 


Due to the asymptotic approximations used, the results for example 3 
are not as accurate as those for example 2; the two sets of results show 
only that the Q of the cavity is substantially the same for the two differ- 
ent modes of oscillation. 


The poles of Y are given by 


Py ’ ia py P, 


P 
ad 


An = An — thn 


| 
2(b - 1) by 


Hy | 
2hd 2(b a)é, 


Applying formulas (6) gives the following values for the n 


th 


branch of 
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the network of the first kind: 


Ky 


m7 FE 


in all of which wo, mnv (hb 3(10° 


» meters 


per second. 


The results can be put in the same form as those obtained for the other 
cavity mode, dealt with in example 2, by employing the “primary con 
stants” of the evlindrical transmission line, viz.: 


Gia 


In terms of these constants, the elements of the n 


branch of the « quiva- 
lent network of the first kind are 


(} a) is' 


») 


a AR 1) 


Cle) (, 
Wi 


»,¢ (a) < 
’ Wir 
(, — ] : 
T-i-0 Phd 


The network is shown in Fig. 12 


As in the preceding example, a leakage element arises, in spite of the 
fact that we assumed initially that go of the air in 


i 


the cavity is zero 
This element accounts for the losses in the inner and outet evlindrical 
walls. 
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\ number of people with whom the above material has been discussed 


have given helpful comments and criticisms. | wish to acknowledge my 
debt in this respect to HE. Nyquist, 8S. A. Schelkunoff, R. M. Foster, 8S. O. 
Rice, J. Riordan and W. EH. Wise 
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Universal Equalizer Chart. D. A. Auspera’. Electronics, 24, pp. 132, 
134, Nov., 1951. 

Modification of familiar Smith chart consolidates on one time-saving plot all 
positive-value solutions to the two general equations for series, shunt, 
and bridged-T audio equalizers. 


Limits on the Energy of the Antiferromagnetic Ground State. P. W. 
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Post-War Achievements of Bell Laboratories, [. O. ki. Buckiey. Bell 
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Filamentary Growths on Metal Surfaces—* Whiskers”. K.G. Compron , 
A. Menpizza', and S. M. Arnoup'. Corrosion, 7, pp. 327-334, Oct., 
1951. (Monograph 1885). 

Filamentary growths have been found on metal surfaces of some of the parts 
used in telephone communications equipment, particularly on parts shielded 
from free circulation of air. The growths are of the same character as those known 
as ‘‘whiskers,”’ which developed between the leaves of cadmium plated variable 
air condensers and caused considerable trouble in military equipment during the 
early part of World War If. An investigation has been under way in an attempt 
to determine the mechanism of growth of the whiskers, found not only on 
cadmium plated parts but also on other metals. This paper summarizes the 
findings to date as revealed by the study of approximately one thousand test 
specimens of different metals, solid and plated, exposed under various environ- 
mental conditions. The study is being extended in the light of the findings which 
have developed during the course of the work 


An Unattended Broad-Band Microwave Repeater for the TD-2 Radio 
Relay System. R. W. Frus’ and K. D. Smirn’. Elec. Eng., 70, pp. 976 
O81, Nov., 1951. (Similar article in The Bell System Technical Journal, 
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October, 1951, Part Il, entitled The DT-2 Microwave Radio System 
reprinted as Monograph 1921 
To meet the stringent requirements of the 4,000-mile transcontinental micro- 
system, a number of new developments had to be included in the 
the repeater stations. The circuits of these uuattended stations, and 
ie maintained, are the subject of this article. 


The Bell System’s Part in Defending the Nation. F. R. Karrr.’. Bell 
Tel. Mag., 30, pp. 141-152, Autumn, 1951. 

Quickly ind accurately checks performance Of private or Common-carriel 
p-mn or f-m mobile telephone transmitters, such as those used in 30 to 44 -me 
highway and 152 to 175-me urban service. Measures r-f power output, audio 
sensitivity, signal-to-noise ratio and harmonic distortion and gives speech 
intelligibility check in few minutes. 


VWobili 1 ransmilter Te sling Set. (i. a KENT ‘ electronics, 24, pp. 106- 
109, Nov., 1951 


1 Neu Hlectrolysis Switch for Unde rground Lead Sheath Cable Drainage 
Systems. V. B. Ptke . Corrosion, 7, bp. 1, Oct., I95L. 


1 fligh Te mperature Stage for the Polarizing Microscope. H;. A: Woop’. 
lim. Mineral., 36, pp. (68-772, Sept.—Oct., 1951. 


| Pre cise Swee p Fre Q(uency Mi thod of Vector Impedance Measure ment, 


ID. A. Atspera'. Proc. I.R.E., 39, pp. 1398-1400, Nov., 1951. (Mono- 


graph IQ] ] 
The impedance of a two-terminal network is defined completely by the inser- 


tion loss and phase shift it produces when inserted between known sending and 


receiving impedances, Recent advances in precise wide-band phase and trans- 
mission measuring circuits have permitted practical use of this principle. Reac- 
tive and resistive impedance components are read directly from a simple graphi- 
rt in which frequency is not a parameter. The basic principle described 
es attractive possibilities in many cases of impedance measurements 


resent methods are in ilequate 


, . ’ . l 
lectron-Vibration Interactions and Supe rconductivity. J. BARDEEN . 
Revs. Modern Phys., 23, pp. 261-270, July, 1951. (Monograph 1912). 


The ('opper Oxide Rectifie r. W. H. BRATTAIN . Revs. Modern Phys., 
23, pp. 208-212, July, 1951. 

It is shown that the conductivity in the ohmic part of the cuprous oxide layer 

in be explained with the usual band picture of semiconductors only by assum- 
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ing the presence of some donor-type impurities in addition to the usual acceptor 
type. The energy difference between the acceptors and the filled band is 0.5 
electron volt, and the total number of impurity atoms is about 10" to LO! pet 
cm}, the number of donors being less than but of the same order as the number of 
acceptors. One finds that the density of ion charge in the rectifving laver is of 
the same order of magnitude as the difference between the donors and ac eptors 
found from the conductivity. The field at the copper-cuprous oxide interface Is 
about 2 x 10¢ volts/em; the height of the potential at the surface as compared 
with the oxide interior is about 0.5 volt: and the thickness of the space charge 
layer about 5.0 x 10> em. The diffusion equation for flow of current through this 
space charge region can be integrated to give the current in terms of the field 
at the interface and the applied potential across the space charge laver. Two 
currents are involved, one from the semi-conductor to the metal (/,) and one 
from the metal to the semiconductor (/,,) which is similar to a thermionic emis- 
sion current into the semiconductor. The net current is, of course, / I i 

One can get this “emission” current (/,,) by dividing the true current by the 
fuctor | - Xp | eVi/kT), where V, is the applied potential This emission current 
depends on the absolute temperature and on the field at the copper-cuprous 
oxide interface. At high fields the logarithm of the current is proportional to the 
square root of the field, and at low fields the current decreases more rapidly 
indicating a patchy surface having small areas of low potential maximum from 
Which all the emission comes when the field is large. 


Effect of Packaging on Corrosion of Zine Plated Equipment. IK. G. 
Compron , S. M. Arnotp, and A. Menpizza’. Corrosion, 7, pp. 365 


372, Nov., 1951. 


Physics as a Setence and an Art. WK. WK. Darrow. Phys. Today, 4, pp 
6-11, Nov., 1951. (Monograph 1914 

The last of six invited papers presented on October 25th during the symposium 
on “physics today” which keynoted the 20th Anniversary Meeting of the 
American Institute of Physics in Chicago 


Ionization by Electron Impact in CO, Nx, NO, and Oo. H. D. Hagsrrum . 
Revs. Modern Phys., 23, pp. 185-203, July, 1951. (Monograph 1916). 

Ionization by electron impact in diatomic gases has been studied in this work 
with a mass spectrometer designed to measure m/e, appearance potential, and 
initial kinetic energy lol each ion observed. Results have been obtained for the 
gases CO, N. NO, and O, with some confirmatory work in H.. Discussion is 
included of the nature and identification of dissociative ionization processes and 
of the retarding potential and appearance potential measurements. Values of 
important quantities such as the dissociation energies of CO, N. , and NO; the 
sublimation energy of C; the electron affinity of O; and the excitation energy of 
©” are determined again by electron impact in this work 
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Equivale nt Te m pe rature of an Electron Beam. M. E. HINES ., Letter 
to the Editor., J. Appl. Phys., 22, pp. 1385-1386, Nov., 1951. 


Bell System Cable Sheath Problems and Designs. F. W. Horn! and 
R. B. Ramsey’., Klee. Eng., 70, pp. 1070-1075, Dec., 1951. (Monograph 
1917). 


Engineering Planning. HU. 8. Osporne. J. Eng. Educ., 42, pp. 121 
125, Nov., 1951. 


leceptance Inspection of Purchased Material. J. . PALMER’ and E, G, 
1). Parerson . Ind. Quality Control, 8, pp. 23-27, Nov., 1951. 


A Nolte on the Partial Differential Equations Describing Steady Current 
Flow in Intrinsie Semiconductors. R. C. Pri’. Letter to the Editor. 
J. Appl. Phys., 22, pp. 1388-13889, Nov., 1951. 


General Theory of Symmetric Biconical Antennas. S. A. SCHELKUNOFF. 
J. Appl. Phys., 22, pp. 1330-1332, Nov., 1951. (Monograph 1922). 

This paper presents the input admittance of a symmetric biconical antenna 
of an arbitrary angle as the limit of a certain sequence of functions. The first 
term of this sequence approaches the exact expressions for the input admittance 
as the cone angle approaches either zero or 90°. For this reason our conjecture is 
that this term represents a good first approximation for all angles. 


Artificial Dielectrics for Microwaves. W. M. Suarpuess’. Proc. I. R. E., 
39, pp. 1389-1393, Nov., 1951. (Monograph 1923). 

This paper presents a procedure for measuring the dielectric properties of 
metal-loaded artificial dielectrics in the microwave region by the use of the 
short-circuited line method. Formulas, based on transmission-line theory, are 
included and serve as guides in predicting the approximate dielectric properties 
of certain loading configurations. 


B. Tk wks 
W. E. Co. 





Contributors to this Issue 


W. ©. Baker, B.S., Washington College, Maryland, 1935; Ph.D., 
Princeton University, 1938; Bell Telephone Laboratories, 1939—. Dr. 
Baker has carried on investigations of the molecular structure and physi- 
cal properties of polymers, particularly the fundamental constitution of 
synthetic rubbers and plastics. Harvard Fellowship, 1936-37 and Proctor 
Fellowship, 1938-39. Member of American Chemical Society, American 
Physical Society, and American Society for Testing Materials. 


J. H. Hetss, B.S. in Ch.E., Newark College of Engineering, 1942; Bell 
Telephone Laboratories, 1934-. Mr. Heiss has devoted his time to study- 
ing experimental wire coating procedures and the test methods involved, 
the experimental production of high polymers (plastics) and the exami- 
nation of their physical properties, and the properties of high polymers 
in solution. Member of American Chemical Society. 


W.S. Haywarop, Jr., A.B., Harvard University, 1943; 5.M., Harvard 
Graduate School of Engineering, 1947; Aircraft Radio Laboratory, June 
December 1948; U.S. Navy, Aviation Electronics Officer, 1944-46; Bell 
Telephone Laboratories, 1947—-. Mr. Hayward has taught telephone 
switching circuit design at the Laboratories and is currently making 
probability studies of telephone traffic. 


Brockway McMI Luan, B.S., Massachusetts Institute of Technology, 
1936; Ph.D., Massachusetts Institute of Technology, 1939; Instructor 
of Mathematics, Massachusetts Institute of Technology, 1936-39; 
Proctor Fellow and Henry B. Fine Instructor in Mathematics, Princeton 
University, 1939-42; U.S.N.R., 1942-46, studying exterior ballistics of 
guns and rockets; Los Alamos Laboratory, Spring 1946; Bell Telephone 
Laboratories, 1946-. Dr. McMillan has been engaged in mathematical 
research and consultation work. Member of American Mathematical 
Society, Institute of Mathematical Statistics, and A.A.A.S. 


R. E. Sraruier, B.B.E., College of the City of New York, 1947; 
M.E.E., Polytechnic Institute of Brooklyn, 1948; U.S. Army 1943-46, 
Communications Officer; Instructor in Electrical Engineering, Polytechnic 
Institute of Brooklyn, Fall, 1950; Bell Telephone Laboratories, 1948 
After completing the Communications Development Training Program, 


109 





HO rHE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1952 


with rotational assignments in the transmission, switching, and apparatus 


departments, Mr. Staehler became a member of a switching group 


concerned with both local and toll signaling. He is working at present on 
a new voice-frequency toll signaling development. Member of Tau Beta 
Pi and Eta Kappa Nu. 


M. W. Zinn, BS. in E.E., Purdue University, 1918; U.S. Army, 
19IS 19. Amer. Tel. and Tel., 1919-1934; Bell Telephone Laboratories, 
1034. As a transmission engineer, Mr. Zinn has been concerned with 
land and submarine loaded cables, telephone instruments, buried cable 
ind submarine cable with amplifiers and special problems. Member of 
\. 1. BE. EB. and Tau Beta Pi 











